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Inelastic light scattering is an intensively used tool in the study of electronic properties of solids. 
Triggered by the discovery of high temperature superconductivity in the cuprates and by new de- 
velopments in instrumentation, light scattering both in the visible (Raman effect) and the X-ray 
part of the electromagnetic spectrum has become a method complementary to optical (infrared) 
spectroscopy while providing additional and relevant information. The main purpose of the re- 
view is to position Raman scattering with regard to single-particle methods like angle-resolved 
photoemission spectroscopy (ARPES), and other transport and thermodynamic measurements 
in correlated materials. Particular focus will be placed on photon polarizations and the role of 
symmetry to elucidate the dynamics of electrons in different regions of the Brillouin zone. This 
advantage over conventional transport (usually measuring averaged properties) indeed provides 
new insights into anisotropic and complex many-body behavior of electrons in various systems. 
We review recent developments in the theory of electronic Raman scattering in correlated systems 
and experimental results in paradigmatic materials such as the A15 superconductors, magnetic 
and paramagnetic insulators, compounds with competing orders, as well as the cuprates with 
high superconducting transition temperatures. We present an overview of the manifestations of 
complexity in the Raman response due to the impact of correlations and developing competing 
orders. In a variety of materials we discuss which observations may be understood and summarize 
important open questions that pave the way to a detailed understanding of correlated electron 
systems. 
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I. INTRODUCTION 
A. Overview 

Raman scattering is a photon-in photon-out process 
with energy transferee! to a target material. Most of the 
light is elastically scattered from the sample, a fraction 
is color-shifted and collected at the detector. 

Light couples to electronic charge in solids and can 
scatter inelastically from many types of excitations in 
a sample, as shown schematically for YBa2Cu306.5 in 
Figure 1. Optical phonons produce sharp peaks at well 
known positions and orientations of the incoming and 
outgoing photon polarizations, while a large broad fea- 
ture centered at much higher energies due to two-magnon 
scattering occurs in compounds with antiferromagnetic 
correlations. This review article largely places empha- 
sis on the electronic Raman scattering continuum upon 
which the phonons and magnons are superimposed. 

Light scatters off of electrons by creating variations 
of electronic charge density in the illuminated region of 
a sample. By observing the frequency shift and polar- 
ization change of the outgoing photon compared to the 
incoming photon, the properties of charge density relax- 
ation is measured. However, measuring the Raman effect 
of photons scattering from electrons is a difficult proposal 
to carry forward. It is made difficult precisely because of 
the coupling of the real photon vector potential and the 
exchange of virtual photons which mediate the Coulomb 
forces between electrons. In simple metals, variations of 
the charge density will be largely screened by the mo- 
bile electrons, and the system of electrons responds col- 
lectively at a characteristic plasma frequency of several 
electron volts. In semiconductors or well-developed band 
insulators, the creation of charge density fluctuations oc- 
curs only via the population of excited states across a 
band gap - again on the scale of several electron volts. 
Therefore it is difficult to investigate the behavior of elec- 
trons at low energies. In fact, hardly any measurements 
of electronic Raman scattering in simple metals exist pre- 
cisely for this reason, and focus on semiconductors is usu- 
ally placed on plasma excitations. Since the dynamics 
of electrons lying near the Fermi surface govern the be- 
havior of transport in most systems, this would give the 
impression that Raman would have but little to offer in 
simple metals and insulators. 

Yet the Raman effect is extremely well suited to study 
electrons in systems with non-trivial electron dynamics. 
First well-studied in the context of breaking Cooper pairs 
in superconductors in the period from 1980 to 1990, the 
field of Raman scattering from electronic excitations has 
grown tremendously over the past few decades to study 
the evolution of electron correlations in a variety of sys- 
tems in which many-body interactions arc essential to 
the physics of new materials and their potential device 
applications. 

Raman spectroscopy has become an indispensable tool 
in the arsenal for understanding many-body physics. One 



of the most celebrated achievements of electronic Raman 
scattering has been the ability to focus on the nature 
of electron dynamics in different regions of the Brillouin 
zone. This distinguishes Raman scattering from most 
other transport and thermodynamic measurements, al- 
lowing the study of the development of correlations in 
projected regions of the Brillouin zone. By simply align- 
ing the polarization orientations of the incoming and 
outgoing photons, charge excitations can be selectively 
mapped and analyzed using group-theoretical symmetry 
arguments. The search for conventional as well as ex- 
otic excitations in strongly correlated matter has been 
greatly enhanced. Raman spectroscopy has provided new 
and valuable insights into unconventional superconduc- 
tivity and collective modes, excitations in charge, spin 
and/or orbitally ordered systems as well as the compe- 
tition between the various ordered phases. In addition, 
new insights into electron dynamics of metal-insulator 
transitions, quantum phase transitions, and the concomi- 
tant quantum critical behavior could be obtained. The 
purpose of this article is to review the essence of these 
new developments in a "snapshot" of the current state of 
investigation. 

The overall agenda of the paper is to provide a vehicle 
to sort through the extensive literature, learn about the 
outstanding problems, and become aware of the level of 
consensus. In order to present a detailed picture of the 
current status of electronic light scattering, other types of 
excitations, such as phonons and magnons, arc mainly ig- 
nored. There have been many reviews on inelastic light 
scattering from phonons and magnons. The reader is 
referred to earlier reviews by Klein (1982b) for a fun- 
damental treatment of scattering from phonons, while 
studies of phonons in high temperature superconduc- 
tors are summarized by Thomsen (1989) and Sherman 
et al. (2003). Recently Lemmens et al. (2003) and Gozar 
et al. (2005b) reviewed magnetic light scattering in low- 
dimensional quantum spin systems and cuprates. Due to 
space limitations we cannot give adequate commentary 
on these exciting and developing fields. 

The outline of our review is as follows. After a brief his- 
torical summary, the fundamental experimental aspects 
and theoretical developments for electronic Raman scat- 
tering are presented in the first part of the article. A gen- 
eral treatise on the theory of electronic Raman scattering 
is given in Section II with a view toward the formalism for 
both weak and strong correlations. Results from model- 
specific calculations can be found in Section II. D. Read- 
ers who are more interested in summaries of experimental 
work may want to skim these sections and skip to Section 
III, where a review of Raman scattering measurements in 
a variety of correlated materials is given with a view to- 
ward common features manifest from strong correlations. 
The presentation is generally organized in systems with 
increasing complexity of correlations and competing or- 
ders. 

In this framework, the canon of work on the high tem- 
perature superconductors in the last part of our review is 
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FIG. 1 Characteristic Raman scattering spectrum taken on YBa2Cu306.5 showing light scattering from phonons (blue lines), 
magnons (red) , and electrons (green) . Courtesy of Matthias Opel. 



presented in Section IV. This detailed part of the review 
is organized in conceptual issues of correlations, super- 
conductivity, normal state properties, and the propensity 
toward charge and spin ordering in various families of the 
cuprates. In all subsections in this part, data on a variety 
of cuprate materials are summarized. 

The review closes with a general discussion of open 
questions for both experimental and theoretical develop- 
ments in Raman scattering, and points out new directions 
in which our understanding of electronic correlations may 
be further enhanced. 

B. Historical Review 

Inelastic scattering of light was discovered indepen- 
dently in organic liquids by Raman and Krishnan (1928) 
and in quartz by Landsberg and Mandelstam (1928) who 
properly explained the observed effect: The energy of 
the incoming photon is split between the scattered one 
and an elementary excitation in the solid. Shortly there- 
after in 1930 C. V. Raman was awarded the Nobel prize, 
and his name was associated with the effect (Fabelinskii, 
1998; Ginzburg, 1998; Pleijel, 1930). Although the phe- 
nomenological description by Smekal (1923) in terms of 
a periodically modulated polarizability qualitatively cap- 
tures the relevant physics including the selection rules, 
the effect is genuinely quantum mechanical as described 
first and ahead of the experimental discovery by Kramers 
and Heisenberg (1925) in context of the dispersion in di- 
electrics. 

Soon after the observation of light scattering from vi- 
brational excitations, Verkin and Lazarev (1948) and 
Khaikin and Bykov (1956) attempted to use the new 
technique for studying electronic excitations. They 
picked one of the most ambitious subjects, i.e., light 
scattering from superconducting gap excitations in con- 
ventional metals. It is not at all surprising that they 



could not succeed. In a seminal paper Abrikosov and 
Fal'kovskii (1961) not only calculated the Raman re- 
sponse of a typical elemental superconductor but also 
demonstrated that the sensitivity in the early experi- 
ments was by approximately 5 or 6 orders of magni- 
tude too low. In 1980 light was finally scattered suc- 
cessfully from superconducting electrons in 2H — NbSe2 
(Sooryakumar and Klein, 1980). Balseiro and Falicov 
(1980) and Littlewood and Varma (1981, 1982) argued 
that the superconducting excitations in this system be- 
come Raman active mainly via their coupling to a charge- 
density wave mode (CDW). After the observation of 
gap excitations in the A15 compounds NbaSn and VsSi 
(Dierker et at, 1983; Hackl et ai, 1982, 1983; Klein, 
1982a) it was clear that light can be scattered directly 
by Cooper pairs (Dierker et ai, 1983; Klein and Dierker, 
1984). Tiitto and Zawadowski (1992) demonstrated that 
both types of coupling contribute. 

Collective excitations of normal electrons were first ob- 
served in semiconductors (Mooradian and Wright, 1966) 
following theoretical studies by Pines (1963), Platzman 
and Tzoar (1964). As a function of doping the plasmon 
peak moves across the phonon energies leading to strong 
electron-lattice interactions. In heavily doped silicon, 
with the plasma energy well beyond the vibration spec- 
trum, the evolution of the phonon line shape (Cerdeira 
et ai, 1973; Fano, 1961) clearly demonstrated the exis- 
tence of an electron continuum. In 1977, fluctuations 
of electrons between pockets of the Fermi surface of sil- 
icon were observed by Chandrasekhar et al. (1977) and 
explained subsequently by Ipatova et al. (1981). In mag- 
netic fields, transitions between Landau levels were found 
(Worlock et al., 1981)^. Strong phonon renormalization 



The subject has been reviewed in detail by Abstreiter et al. 
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effects also oceur in metallic alloys with A15 structure 
(Schicktanz et at, 1980; Wipf et at, 1978). The origin 
of the broad continuum, which interacts with phonons 
and is redistributed below the superconducting transi- 
tion (Klein and Dierker, 1984), is certainly electronic but 
as of today is still not fully understood. 

The full power of the method became apparent after 
the discovery of copper-oxygen compounds (Bednorz and 
Miiller, 1986) with superconducting transition tempera- 
tures above 100 K. It turned out that, in contrast to 
infrared spectroscopy, momentum dependent transport 
properties can be measured with Raman spectroscopy, 
since different regions of the Brillouin zone can be pro- 
jected out independently by appropriately selecting the 
polarizations of the incident and scattered photons (Dc- 
vereaux et al, 1994a). New theoretical ideas were not 
only applied to the superconducting but also to the nor- 
mal state. The spectra extend over energy ranges as large 
as electron Volts (eV) (Bozovic et al., 1987; Cooper et al., 
1988a, b; Kirillov et al., 1988) and are similar to those in 
the A15s or in rare earth elements (Klein et al., 1991). 
Both elastic (Zawadowski and Cardona, 1990) and in- 
elastic (Itai. 1992; Kostur, 1991; Virosztck and Ruvalds, 
1992) relaxation of electrons indeed produces light scat- 
tering over such a broad range of energies. It soon be- 
came clear that a continuum extending over an cV cannot 
originate from elastic scattering, but only from inelastic 
processes or interband transitions. However, it has not 
been straightforward to pin down the types of interac- 
tions. 

At very low energies, spin- (Yoon et al., 2000) and 

charge-ordering fluctuations were reported in manganites 
and, respectively, in ladder compounds (Blumberg et al., 
2002) and in the cuprates (Venturini et al, 2002b). The 
response should not be confused with that of an ordered 
spin/charge-density-wave (SDW/CDW) state (Benfatto 
et al., 2000; Klein, 1982c; Zeyher and Greco, 2002). Since 
a characteristic energy decreases rather than increases 
upon cooling (Blumberg et al, 2002; Caprara et al., 2002, 
2005; Venturini et al, 2002b; Yoon et al, 2000), which is 
in clear contrast to typical order parameter behavior. 

Along with the early studies of charge excitations, 
Flcury and coworkers observed Raman scattering from 
spin waves in antiferromagnetically ordered FeF2, MnF2, 
and K2NiF4 (Fleury et al., 1966, 1967, 1970). Elliot and 
Loudon (1963) and Flcury and Loudon (1968) presented 
a detailed theoretical description which allowed them to 
semi-quantitatively understand the spectral shape and 
the cross section. With the discovery of the cuprates by 
Bednorz and Miiller (1986) also this field experienced a 
renaissance in particular for the study at low doping close 
to the antiferromagnetic Neel state (Gozar et al., 2004, 
2005b; Lyons et al., 1988; Sugai et al., 1988; Sulewski 
et al., 1991). Raman scattering is probably the most 



precise method for determining the exchange coupling J 
though the theoretical understanding is still incomplete. 
In this context, spin-Peierls systems (Loosdrccht et al., 
1996) and ladder compounds (Abrashev et al, 1997) 
shifted very much into the focus of interest (Dagotto, 
1999). 

Very recently, light scattering from "orbitons", i.e. 
from a propagating reorientation of orbitals, has been 
proposed to explain new modes in the the Raman spec- 
tra (Saitoh et al., 2001). However, there is no agreement 
yet on whether or not orbitons can be observed indepen- 
dent of phonons or other excitations (Choi et al., 2005; 
Griininger et al, 2002; Kriiger et al, 2004). 



C. What Can One Learn from Electronic Raman 
Scattering? 

We give now a qualitative introduction into the rela- 
tionship between Raman spectroscopy and other experi- 
mental techniques. We begin by drawing the distinction 
between one-particle and many-particle properties. 

Typically, electronic states in solids are characterized 
by their energy dispersions as well as the characteristic 
lifetime of an electron placed into such a state. This 
state is characterized by the single particle propagator 
or Green's function for the electron. 



G(k,a;) 



w - - S(k, w) 



(1) 



(1984) 



Here denotes the bare energy band dispersion cal- 
culated from a solvable model. S represents the elec- 
tron self-energy which encompasses all the information 
pertaining to interactions of the single electron in state 
k to all other excitations of the system. Usually the 
self energy can only be obtained via approximate meth- 
ods. Some of these approximations are quite good - 
such as electron-phonon interactions in metals (known 
as Migdal's approach (Migdal, 1958)) for example, while 
others are more difficult - such as the Coulomb interac- 
tion between other electrons. The self energy is a com- 
plex function, S = S' -|- iS", which, in general, depends 
on temperature, momentum and energy. The real part of 
the self energy determines how the energy dispersion ^k 
is renormalized by the interactions while the imaginary 
part determines the lifetime of the quasiparticle placed 
into the state k. 

The spectral function is directly related to the analyti- 
cally continued electron's Green's function for frequencies 
on the real axis via the replacement iw w + i5 

A(k,w) = -i limG"(k,a;-Fi(5) (2) 

which is measurable via modern angle-resolved photoe- 
mission (ARPES) techniques and has provided an im- 
mense amount of information in strongly correlated sys- 
tems (Campuzano et al, 2002; Damascelli et al, 2003). 
For non-interacting electrons, A(k, w) is a 5-function 
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peaked at the pole of the propagator when the fre- 
quency uj equals the bare band energy ^k- Interac- 
tions broaden the spectral function and give it non-trivial 
temperature and frequency dependences as well as non- 
trivial anisotropics in momentum-space if the interac- 
tions among electrons are anisotropic. The spectral func- 
tion describes real electrons, hence integrals over all en- 
ergies must obey sum rules, such as (i) / (Lj A{\l, uo) — 1 
and (ii) ^ dLoJ{uj)A(\<i,uj) = n(k) with the Fermi Dirac 
distribution f{oj) and the momentum distribution func- 
tion n(k). 

If the electronic interactions are weak, one usually uses 
the nomenclature of Landau and refers to dressed quasi- 
particles replacing the electron as the fundamental exci- 
tation in the solid. These interactions may be character- 
ized by the residue of the pole (usually denoted by Zk) 
and the quasiparticle effective mass m*/mb = {Zk)~^ 
with mb the bare band mass for quasiparticles lying near 
the Fermi surface. is related to the real part of the 
self energy E' which can be expanded for electrons near 
the Fermi surface as E'(k, w) ~ So(k) + u)dT,'{k,u) = 
Q)/du). According to Luttinger's theorem, the Fermi 
surface average of Sg(k) vanishes. The enhancement 
of the quasiparticle mass over the band mass can be 
written as m* /rrih = (1 — dT,' /duo). One often defines 
m*/my, = 1 + A with the dimensionless coupling con- 
stant A > (see, e.g., (Ashcroft and Mermin, 1976)). 
Zk = (1 — dH' /dw)~^ is always smaller than 1, reflecting 
the fact that even for a; = and T ~ only a fraction 
Zk of the spectral weight (coherent part) is in the pole of 
j4(k, ui) while I — Z]^ (incoherent part) is distributed over 
larger energy scales. Equivalently, .Zk < 1 is the disconti- 
nuity at ki;- of the zero-temperature momentum distribu- 
tion function n(k). If Zk approaches zero (oc 1/lnw) the 
system is referred to as a marginal Fermi liquid (Varma 
et al., 1989a), and sum rule (i) is exhausted only at en- 
ergies much larger than ^k- Thus, knowledge of the self 
energy is an important requisite to understanding many- 
body interactions. 

For this reason, single particle methods such as 
ARPES, electron tunneling and specific heat measure- 
ments have been applied extensively to study correlated 
electron systems. Very much stimulated by the discov- 
ery of superconductivity in the cuprates (Bednorz and 
Miiller, 1986) ARPES and tunneling spectroscopy have 
developed more rapidly than any other method in the 
last decade. ARPES data has given unprecedented in- 
sight into momentum-resolved single electron properties 
and their many body effects (Campuzano et al., 2002; 
Damascelli et al., 2003), while tunneling measurements 
have provided information on pairing (Mandrus, 1991; 
Renner, 1995; Zasadzinski et al., 2002) and have recently 
elucidated many issues of nanoscale inhomogeneities in 
the cuprates and their connection to superconductivity 
(Fang et al., 2006; Hanaguri et al., 2004; Hoffman et al., 
2002; Howald et al, 2003; Kivelson et al, 2003; McEl- 
roy et at, 2003, 2005; Vershinin et al., 2004). Detailed 
knowledge of phase transitions in the cuprates has been 



obtained from specific heat studies (Loram and Tallon, 
2001; Moler et al., 1994; Roulin et al., 1998). Due to 
space limitations we only have listed some of the later ref- 
erences of important experimental papers which we hope 
serve as an entry point for the reader to search backwards 
in time to follow the developments. 

Yet knowledge of the spectral function and single- 
particle excitation spectra do not yield information about 
how the electrons may transport heat, current, entropy, 
or energy. For this one needs two-particle correlation 
functions for charge or spin which can be measured by, 
e.g., ordinary and heat transport, optical spectroscopy, 
neutron and light scattering. As an example for such a 
correlation function we consider a standard expression 
for the generalized Kubo susceptibility x" bi^) weakly 
interacting, isotropic normal electrons (see, e.g., Mahan 
(2000)), 

<b{^) = TpE^k^k / —G"{k,u;)G"{k,iv + n) 



x[f{u;)-f{oj + n)]. 



(3) 



Here V is the volume, ak,6k are the bare vertices rep- 
resenting quasiparticle charge (ok = 1) or current (ok = 
jk = ek) correlation functions, and the factor 2 accounts 
for spin degeneracy. The absorptive part of the conduc- 
tivity cr' = Xj'ji^)/^ measures essentially a convolution 
of occupied and unoccupied states. For electrons weakly 
interacting with impurities the conductivity can readily 
be calculated to exhibit a Lorentzian dependence on fl 
represented by 



(t'(0) 



(To 



1 + (Or)2 



(4) 



where ctq is the dc {^l = 0) conductivity, and the relax- 
ation time T = —h{2T,")~^ controls both the width of 
the spectral function and the conductivity as a function 
of frequency.^ A very similar expression is found for light 
scattering^. Thus in the case non-interacting electrons, 
single and two-particle correlation functions give similar 
results. 

This is also true by and large if weak but essen- 
tially isotropic interactions lead to an energy dependent 
E"(a;) and, for causality, to a finite T,'{u). Gotze and 
Wolfle (1972) and, more phenomenologically, Allen and 
Mikkelsen (1977) (for a more recent reference see (Basov 
and Timusk, 2005)) discuss how this generalization mod- 
ifies the response given by Eq. (4), which is then often 



^ Note that Eq. (3) does not return the proper transport lifetime Tt 
which differs by a factor of typically 1 — cos 6 with the scattering 
angle since events with S do not contribute to the resis- 
tivity. This deficit must be taken care of by vertex corrections 
(Mahan, 2000). 

^ The case of Raman scattering is described in detail in reference 
(Opel et al, 2000) and touched upon briefly in section IV.D.l. 
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referred to as the extended Drude model. However, in- 
teracting systems require some care. For example, in 
superconductors, both the single and the two-particle 
responses yield the energy gap. Yet two-particle corre- 
lation functions also have coherence factors which can 
be crucially important to determine the gap symmetry 
in unconventional systems. Generally, collective modes 
(such as the plasmon or excitons) appear directly in two- 
particle correlation functions but only indirectly in the 
spectral function. 

Sometimes the results from single- and two-particle 
measurements can be qualitatively different, even for 
non-interacting electrons. As an illustration, we consider 
first the metal-insulator transition occurring in a system 
of otherwise non-interacting electrons in a disordered en- 
vironment (Anderson transition). Here, backscattering 
of electrons from impurities leads to destructive phase 
interference, and the electrons become localized once a 
critical concentration of impurities is in place in three 
dimensions. Thus while the conductivity is critical and 
vanishes at the metal-insulator transition, the spectral 
function, or equivalently the density of states, is uncrit- 
ical. This distinction becomes even more pronounced if 
the electron interactions are strong and anisotropic, and 
the bare vertices along with the Green's functions en- 
tering into Eq. (3) must be renormalized by the strong 
interactions. 

As a second example, in the spinless Falicov-Kimball 
model light d— electrons strongly interact with localized 
/—electrons and are characterized by the Hamiltonian 
(Falikov and Kimball, 1969) 

-/i^(c|ci -Fwj) -|-f/^c|ciWi (5) 

where c| (cj) create (destroy) a conduction electron at 
site i, Wi is a classical variable (representing the local- 
ized electron number at site i) that equals or 1, t* is 
a renormalized hopping matrix element that is nonzero 
between nearest neighbors on a hypercubic lattice in D- 
dimensions, and U is the local screened Coulomb inter- 
action between conduction and localized electrons. 
denotes a sum over sites i and nearest neighbors j. E{ 
and /i are adjusted to set the average filling of conduction 
and localized electrons. This model has been solved ex- 
actly for electrons on a hypercubic lattice in the limit of 
large coordination number (Freericks and Zlatic, 2003b). 
The system undergoes a metal-insulator transition (MIT) 
at half-filling (one electron per site) if the interaction U 
is beyond a critical value Uc- On either side of the metal- 
insulator transition, the density of states is temperature 
independent (van Dongen, 1992), while the conductivity 
has a strong temperature dependence (Pruschke et al., 
1995) showing the development of the MIT. 

In systems with strong and anisotropic interactions, 
the differences between single and two-particle properties 



are inescapable. This is been borne out in the cuprates 
by the large amount of work using optical (Homes et 
al., 2004) and thermal conductivities (Sutherland et al., 
2005), resistivities (Ando et al., 2004), nuclear magnetic 
resonance (NMR) (AUoul et al., 1989) and electron spin 
resonance (ESR) (Janossy et al., 2003). These experi- 
ments have revealed basic properties of strongly corre- 
lated systems and have emerged as key elements to char- 
acterize the complex behavior of the high- Tc cuprates. 

Yet these two-particle measurements are largely insen- 
sitive of anisotropies, as they measure Brillouin zone av- 
eraged quantities. As a result they reveal the behavior of 
quasiparticles having the highest velocities which, in the 
cuprates, are the quasiparticles near the nodal regions 
of the Brillouin zone. As far as carriers are concerned, 
the momentum dependence of neutron scattering serves 
mainly to measure spin dynamics in different regions of 
the Brillouin zone. 

In this review, we illustrate that Raman spectroscopy 
gives complementary information to all of these measure- 
ment techniques, and also may provide detailed informa- 
tion of charge and spin dynamics of electrons in different 
regions of the Brillouin zone. This is due to the polariza- 
tion selection rules. As with phonon scattering (Hayes 
and Loudon, 2005), simple group theoretic symmetry ar- 
guments can be used to focus on electron dynamics in 
different regions of the Brillouin zone. For Raman scat- 
tering (aic^k) is replaced with 7^ which, in certain limits, 
may be represented by 7k = I]^.^ e^ef^a^^^/^'.^^/j^afc^, 
with e*'^ the incident, scattered Ught polarization vec- 
tors.^ Apart from energy-independent scaling factors and 
vertices with different k dependences there is an extra 
factor 1/51 between Raman and infrared response. It has 
been shown first by Shastry and Shraiman (1990) and ex- 
plicitly demonstrated within Dynamical Mean Field The- 
ory (DMFT) by Freericks and Devereaux (2001) that, un- 
der certain restrictions, there is a simple correspondence 
between conductivity and Raman response, 

Qa'iQ) (X x"^{Q), (6) 

highlighting that electronic Raman scattering measures 
transport properties. However, even the simple form of 
the vertices given above shows that a coincidence can be 
expected only for an isotropic material. In anisotropic 
systems, light scattering can sample parts of the Fermi 
surface which are unaccessible for infrared spectroscopy. 



D. State of the Art Experimental Technique 

Over decades Raman scattering was predominantly 
used for the study of molecular and lattice vibrations 
which produce isolated and typically narrow lines in the 



* For a microscopic derivation see section II. A. 
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spectra (see Fig. 1). The lines are used as probes which 
sensitively react to changes in the environment of the vi- 
brating atoms. Similar considerations are at the heart of 
magnetic resonance techniques such as NMR and ESR. 

If light is scattered from electrons in solids, the spectra 
are usually continuous (Fig. 1). To study their evolution 
as a function of a control parameter, such as temper- 
ature, doping, magnetic field or pressure, is rather in- 
volved since the overall shape and not the position of 
well-defined lines matters. In addition, typical cross sec- 
tions per unit solid angle and energy interval are smaller 
by several orders of magnitude than those of vibrations. 
Electronic Raman scattering in a metal typically pro- 
duces one energy-shifted photon per s, meV, and sr (unit 
solid angle steradian) out of 10^^ incoming ones. The 
low efficiency is particularly demanding in studies at high 
pressure since additional losses and complications such as 
fluorescence and birefringence arise from the windows, 
which are typically diamond anvils. Although there were 
successful early experiments (Zhou et ai, 1996) the avail- 
ability of synthetic diamonds brought substantial ad- 
vances (Goncharov and Struzhkin, 2003). 

There are three inventions which finally produced the 
required sensitivity: (i) the laser as an intense light 
source providing lines of high spectral purity in a wide en- 
ergy range, (ii) as an early application of the laser, holo- 
graphically fabricated gratings without secondary images 
(ghosts) and an extremely low level of diffusely scattered 
light and, finally, (iii) the invention of charge-coupled 
devices (CCD) as a location-sensitive detector with an 
efficiency at the quantum limit and negligible dark count 
rate. 

Gratings have an extremely well defined number of 
lines per unit length (cm in the cgs system). This is 
the origin of the energy unit cm"-'^. The following con- 
versions are frequently used: 

1 meV ^ 11.604 K 

1 meV = 8.0655 cm"^ 
cm" ^ 

fee = 0.69504 

Since the CCD has a high spatial resolution down to a 
few /zm it is a superior replacement of the photographic 
plate. It facilitates recording complete spectra in a single 
exposure with energy ranges from meV up to approxi- 
mately 1 eV, depending on the desired resolution. 

The essentials of a setup for inelastic light scattering 
with polarized photons are shown schematically in Fig. 2. 
The coherent light at energy Hui from the laser (Ar+ and 
Kr+ gas lasers are still very popular) is spatially filtered. 
A prism monochromator (PMC) selects the desired fre- 
quency and suppresses incoherent photons from the laser 
medium. A combination of a A/2 retarder and a polar- 
izer (PI) facilitates the preparation of a photon flux of 
a well-defined polarization state and intensity. For ex- 
citation, the polarization inside the sample counts. The 
same holds for the selection of the proper polarization 
for the scattered photons at hujs- The best results are 




FIG. 2 Schematic drawing of the light path. The laser light 
with energy hiOi is first spatially filtered. A prism monochro- 
mator (PMC) is to suppress the plasma lines of the laser 
medium, only the coherent one at hwi passes the slit. The po- 
larization is prepared with polarizer PI and the Soleil-Babinet 
compensator (SB). The A/2 retarder in front of PI allows one 
to adjust the power. Before hitting the sample the light is once 
again spatially filtered to maintain an approximately Gaus- 
sian intensity profile in the spot. If the angle of incidence 
is not close to zero phase shift effects at the sample surface 
must be taken into account since the polarization inside the 
sample is important. High speed optics collects the scattered 
light. The polarization state is selected by a A/4 retarder and 
polarizer P2. The A/2 retarder in front of the entrance slit 
rotates the light polarization for maximal transmission of the 
spectrometer (here single stage). Except for the compensator 
most retarders and polarizers work also for light propagating 
at small angles (up to approximately ±5°) with respect to the 
optical axis. The configuration shown here is usually referred 
to as back-scattering geometry since incoming and outgoing 
photons have essentially opposite momenta in the sample. 

obtained by using a crystal polarizer (PI, e.g. of Clan- 
Thompson type) and a Soleil-Babinet compensator for 
the incoming light and an achromatic A/4 retarder and 
another crystal polarizer (P2) for the scattered light. In 
this way all states, including circularly polarized ones, 
can be prepared. The A/2 retarder in front of the en- 
trance slit of the spectrometer rotates the polarization 
into the direction of highest sensitivity. 

Since we wish to discriminate between the 10^°^"'^^ elas- 
tically scattered photons at huji and the few Raman pho- 
tons at fiujs at very small shifts hfl — \hLUi — hujsl < 
1 meV, a single monochromator is insufficient. A modern 
instrument for Raman scattering in metallic samples has 
three stages consisting of essentially independent grat- 
ing monochromators. The first two are usually subtrac- 
tively coupled and select a band from the spectrum of 
inelastically scattered photons. The third stage disperses 
the band transmitted through the two stages of the pre- 
monochromator into a spectrum which is recorded by the 
CCD. In this configuration, the dispersion is given only 
by the third stage while the first two discriminate the 
elastically scattered laser light. If all stages are coupled 
additively, the resolution is improved by a factor of 3. 
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Because of losses at the mirrors and gratings, only 15 to 
20 % of the photons entering the spectrometer arrive at 
the detector. 

Since very interesting physics is going on at energies 
even below 1 mcV (sec, e.g., section IV. D. 3) the discrim- 
ination is a cardinal point. Out of the two options only 
the premonochromator gives satisfactory results below 
10 meV. The price one has to pay is a loss of intensity 
of approximately 60 %. Alternatively, for energy shifts 
above 10 meV an interferometric notch filter can be used. 
The latter device is widely used for commercial applica- 
tions which develop rapidly since the introduction of the 
CCD. The main fields are quality control and analytics. 

At finite temperatures, T > 0, inelastically scattered 
light is found on either side of Huji. As a consequence of 
time-reversal symmetry and for phase space arguments 
the energy gain (Anti-Stokes) and loss (Stokes) spectra 
are related by the principle of detailed balance (equiva- 
lent to the fluctuation-dissipation theorem) (Landau and 
Lifshitz, 1960; Placzek, 1934) 

with A^sT(AS) s-nd hil the rate of photons per unit time 
collected on the Stokes (Anti-Stokes) side and the energy 
transferred to the system, respectively. Equation (7) can 
be used to determine the temperature of the laser spot. 

If spectra are measured in large energy ranges, the sen- 
sitivity of the instrument has to be taken into account. 
To this end, the spectral response of the whole system, 
including all optical elements between the sample and the 
entrance slit of the spectrometer, the spectrometer itself, 
and the detector, must be calibrated. This is best done 
by replacing the sample with a continuous light source 
of the same size as the laser spot with a known spectral 
emissivity. A continuous source is of crucial importance 
for including the energy dependence of the dispersion in 
addition to the bare transmission. In addition, the fre- 
quency dcpendcnic:c of the; sample's index of refraction, 
\/e = n + ik, requires attention in order to get the inter- 
nal cross section. 

The main limitations of present commercial systems 
come from geometrical aberrations of the spectrometer 
optics and from the relatively low total reflectivity of 
the large number of mirrors. It is a matter of resources 
to improve these caveats. Recently, an improved type 
of triple spectrometer with aspherical optics has been 
described by Schulz et al. (2005). CCDs and gratings are 
close to the theoretical limits. 

For many studies, light sources with continuously ad- 
justable lines in an extended energy range would be de- 
sirable. This holds particularly true for organic materi- 
als (e.g. carbon nanotubcs or proteins) which have rel- 
atively sharp resonances in the visible and the ultravi- 
olet. The synchrotron, free electron "lasers", as well as 
dye and solid-state lasers are developing rapidly and will 
gain influence on the field of Raman spectroscopy in the 



near future. The same holds for near-field techniques 
(Hartschuh et at, 2003) which are capable of improving 
spatial resolution by at least an order of magnitude below 
the diffraction limit. 



II. THEORY OF ELECTRONIC RAMAN SCATTERING 
A. Electronic Coupling to Light 

The aim of this section is to formulate the theoretical 

treatments for inelastic light scattering in general. Much 
has been done in the development of theories for Raman 
scattering, particularly in semiconductors and supercon- 
ductors. Excellent and early reviews of electronic Raman 
scattering have been given by Klein (1983) and Abstreiter 
et al. (1984) focusing on semiconductors. More recent 
reviews by Devereaux and Kampf (1997) and Sherman 
et al. (2003) have focused on theory in superconductors 
with applications towards the cupratcs. We outline the 
general formalism for treating systems with weak and 
strong correlations and return to a discussion of various 
theoretical models in connection with materials in the 
following section. 



1. General Approach 

We first consider a Hamiltonian for N electrons cou- 
pled to the electromagnetic fields (Blum, 1970; Pines and 
Nozieres, 1966): 

H = ^ ^ 2m + -f^Coulomb + -fffields, (8) 

i 

where p = — zfi.V is the momentum operator, e is 
the magnitude of the elementary charge (the electronic 
charge is (j,. = — e) and c the speed of light. A{ri) is the 
vector potential of the field at space-time point ri and 
m the electron mass, -ffcouiomb represents the Coulomb 
interaction and iJficids the free electromagnetic part. We 
use the symbol A to denote operators. We expand the 
kinetic energy to obtain 

H = ^f' + ^ [pi • A(ri) + Ain) ■ Pi 

i 

+ ^EA(rO-A(rO, (9) 

i 

with H' = Ha + //fields and Hq = (l/2m) J2^ Pf + 
-ff Coulomb- Generally we choose |a) to denote eigenstates 
of Ho with eigenvalues £„: Ho\a) = Ea\a). The eigen- 
statc is labeled by all the relevant quantum numbers for 
the state, such as combinations of band index, wavevec- 
tor, orbital and/or spin quantum numbers, for example. 
The eigenstates may be considered to be Bloch electrons 
when the electron-ion interaction is included in Hq, as 
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plane- wave states if it is neglected, or may represent Hub- 
bard states if iJcouiomb is taken to include short-range 
Hubbard- like interactions between electrons. 

The electromagnetic vector potential can be expanded 
into Fourier modes A(ri) = X]q s*'* "^'-^q- second 
quantized notation, the electromagnetic field operator 
takes the form (Mahan, 2000) 



-* 1 1 

GqAqJ 



(10) 



with V the volume and a^, Oq are the creation, anni- 
hilation operators of transversal photons with energy 
hiijq = hc\(\\ having a polarization direction denoted by 
the complex unit vector eq. 

Electronic Raman scattering measures the total cross 
section for scattering from all the electrons illuminated 
by the incident light. The differential cross section is 
determined by the probability that an incident photon 
LOi is scattered into a solid-angle interval between U, and 
U,-\'d£L and a frequency window between and LOs + dLOg. 
A general expression for the differential light scattering 
cross-section is given via the transition rate R of scatter- 
ing an incident (q^, w^, eq'') photon into a outgoing state 



= hri—R. 



(11) 



Here, tq = /mc? is the Thompson radius, and R is 
determined via Fermi's Golden Rule, 

with [3 = Z the partition function, and Mpj — 

{F\M\I) where M is the effective light scattering opera- 
tor. The sum represents a thermodynamic average over 
possible initial and over final states with k vectors in the 
solid angle element dVt of the many-electron system hav- 
ing energies Ei , Ep , respectively. Here = Wi — is 
the transferred frequency and we denote q = — the 
net momentum transfered by the photon. Multiplying 
Eq. (11) by the incident photon flux gives the number of 
scattered photons per second into the solid angle incre- 
ment dD, within the frequency range duig, while multiply- 
ing Eq. (11) by LOg/^i gives the power scattering cross 
section (Klein, 1983). ^ 



From here on we consider the case relevant to Raman 
scattering in the visible range with photon energies of 
typically 2 eV. Since the momentum transferred to the 
electrons, q ^ 1/5, with 5 the skin depth at the light 
energies (Abrikosov and Fal'kovskii, 1961), is much less 
than the relevant momentum scale of order fcp, the Fermi 
momentum in metallic systems, the limit 5 ^ is a good 
approximation in practically all cases. ^ However, finite 
q should be considered if incident light from frequency- 
doubled or synchrotron radiation is used where transi- 
tions between initial and final states at finite q can be 
probed. Then the structure of the Landau particle-hole 
continuum in weakly correlated systems or transitions 
across a finite q Mott gap in strongly correlated insula- 
tors can be studied.^ 

Mpj has contributions from either of the last three 
terms in Eq. (9): the first two terms coupling the elec- 
tron's current to a single photon and the third term cou- 
pling the electron's charge to two photons. This is shown 
in the schematic cartoon in Figures 3 and 4. Here we con- 
sider two bands - one partially filled and the other com- 
pletely filled - in which the incident photon excites an 
electron from either the partially or the completely filled 
band, shown in Figures 3 and 4, respectively. In the non- 
resonant intraband case, the photon gives up part of its 
energy to leave behind a particle-hole pair, while in the 
interband case, an intermediate state is involved, which 
decays via a particle from the partially filled band into 
the hole left behind in the filled band. The latter scat- 
tering may be resonant if the incident or emitted photon 
energy corresponds to that of the energy gap separation, 
otherwise it is non-resonant. In this simple cartoon, one 
can see that excitations lying near the Fermi surface are 
predominantly probed by non-resonant intraband scat- 
tering while excitations involving transition between dif- 





INITIAL STATE FINAL STATE 

▲ ▲ 



^ We note that Eqs. (11) and (12) describe scattering inside the 
materiah Trivial (Fresnel-formulas) and non-trivial (q^ integra- 
tion) (Abrikosov and Fal'kovskii, 1961; Fal'kovskii, 1990, 1991) 
transformations, which we do not discuss here, are required to 
fully describe the cross section outside. The integration origi- 
nates from the lack of momentum conservation perpendicular to 
the surface of an absorbing medium and can change the spectra 
qualitatively. From here on, is always the energy transferred 
to the system. 



FIG. 3 Cartoon showing light scattering via non-resonant in- 
traband scattering. 



^ The applicability of the 5 = limit is discussed in more detail at 
the beginning of section III. 

For a review of relevant work in this regard, the reader is referred 
to references (Devereaux et al, 2003a, c; Kotani and Shin, 2001; 
Platzman and Isaacs, 1998). 
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ferent bands - such as the lower and upper Hubbard band 
for example - are probed by intermediate state scattering. 
The Feynman diagrams representing these contributions 
to Mpj are shown in Figure 5. 

Referring to Eq. (9), the current coupling has odd spa- 
tial symmetry and involves single photon emission or ab- 
sorption, while the second term is even in parity and 
involves two photon scattering of emission followed by 
absorption and vice-versa. The cross-section or the tran- 
sition rate is thus determined via Fermi's Golden rule by 
the square of the matrix elements shown in Figure 5. 

The resulting Feynman diagrams of the contributions 
to the cross section are shown in Figure 6. However, not 
all of them give rise to inelastic light scattering. Some 
of these terms vanish either because they represent con- 
tributions to the renormalized photon propagator (Fig- 
ure 6 (a)), or due to parity arguments (Figures 6 (b)-(d)) 
in the limit of small q scattering. The remaining terms 
can be classified as non-resonant (Figure 6 (e)), resonant 
(Figure 6 (h)-(j) and mixed terms (Figure 6 (f)-(g)), since 
in the former case the initial and final states must share a 
large sub-set of quantum numbers, while the other terms 
can involve transitions through intermediate states well 
separated in energy and distinct from the initial and fi- 
nal states. However, we remark that the response is only 
truly resonant if the photon energies are tuned to the en- 
ergy gap between intermediate and initial or final states. 

To obtain a general expression for the matrix element 
Mpj for Raman scattering, we use second quantized 
notation for the fermions in which the single-particle 
wave function and it's conjugate are given by V'(r) — 

EaCa<Pa(r) and V^(r) = EaCl</'a(l"): ^ith Lf.Lp* the 

eigenstates of the Hamiltonian Hq. Electron states a,/3 
are created, annihilated by cl^,C0 respectively, and the 
indices refer to the quantum number associated with the 
state, such as the momenta and/or spin states. The ma- 
trix element Mp / can thus be written as 



M. 



F.I 



es'^PaA^){F\cicp\ 



I) 



a,l3 



INITIAL STATE INTER. STATE FINAI STATE 




FIG. 4 Cartoon showing light scattering via interband tran- 
sitions. 



I' a, a', 13, 13' 



clc(3> 



I) 



El — E^ -\r huji 



El — El, 



(13) 



Here |/), | F),\i^) represent the initial, final and inter- 
mediate many-electron states having energies Ejp^, re- 
spectively. The many-electron states could be labeled by 
band index and momentum as, for example, for Bloch 
electrons. They may also consist of core and valence 
electrons on selected atoms for x-ray scattering, or may 
represent states of the many-band Hubbard model for 
correlated electrons. pa,^(q) = / d^i'ip'^{r)e^'^''^ Lp/3{r) = 
{a |e*'i "'| P) is the matrix element for single-particle den- 
sity fiuctuations involving states a, (3. The momen- 
tum density matrix element is given by Pa.f3i<ii.s) = 
(a |p • ei^se^"'''^ ''! /?). The first term in the expression 
arises from the two-photon scattering term in Eq. (9) 
in first order perturbation theory. The remaining terms 
arise from the single-photon scattering term in Eq. (9) 




b) 








d) 



FIG. 5 Feynman diagrams contributing to the effective light 
scattering operator M. The top diagram represents single- 
photon absorption, the second two photon scattering, while 
photon emission (absorption) followed by absorption (emis- 
sion) is shown in the 3'"'^ (4*'^) diagram from top. The panels 
on the right are the time-reversed partners of the left dia- 
grams. 
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in second order via intermediate states v and involve dif- 
ferent time orderings of photon absorption and emission. 
The p • A coupUng does not enter to first order since the 
average of the momentum operator is zero. 



2. Importance of Light Polarization 

At this point, httle progress can be made in evaluating 
the matrix elements for Raman scattering without speci- 
fying the quantum numbers of the electronic many-body 
states. Yet, from Eq. (13), one can apply symmetry ar- 
guments to view what types of excitations can be created 
by the incident photons. In this subsection, we employ 
a general set of symmetry classifications and put specific 
emphasis on models in later subsections. 

The first term in Eq. (13) only arises if the incident and 
scattering polarization light vectors are not orthogonal, 
as electronic charge density fluctuations are created and 
destroyed along the polarization directions of the incident 
and scattered photons. Thus, for instance, this term does 
not probe electron dynamics in which the charge density 
relaxes in a direction orthogonal to the incident polariza- 




FIG. 6 Feynman diagrams contained in the cross-section, (a) 
gives a renormalized photon propagator in the solid, while (b)- 

(d) vanish due to parity and thus (a)-(d) do not contribute 
to inelastic light scattering. Of the remaining contributions, 

(e) refers to intra-band (sometimes "non- resonant" ) scatter- 
ing within a single band, (f) and (g) are referred to as "mixed" 
contributions while (h)-(j) describe transitions within a single 
or between different bands via intermediate band states. If 
the light energy is equal or close to the energy difference of 
the states involved resonance effects with a strong enhance- 
ment of the cross section occur. Therefore the contributions 
themselves are sometimes referred to as "resonant". 



tion direction. 

In the limit of small momentum transfer q 0, the 
matrix element simplifies to Pa.fiisi ^ 0) = i^^p and thus 
this term gives rise to scattering from fluctuations of the 
isotropic electronic number density. 

The remaining terms in Eq. (13) have contributions re- 
gardless of photon polarization directions. However, one 
can further classify scattering contributions by separat- 
ing the sum over intermediate states \ v) into states which 
share some quantum numbers with the initial many-body 
state, such as band index, and states which do not. Since 
the photon momenta are much smaller than the relevant 
electron momenta, contributions of the terms where the 
intermediate states include the initial states are roughly 
a factor of v-pjc smaller than the first term in Eq. (13) 
and can be neglected (Pines and Nozieres, 1966; Wolff, 
1966). However contributions where \v) includes higher 
bands cannot in general be neglected, particularly if the 
energy of the incident or scattering light lies near the en- 
ergy of a transition from the initial state Ei to an inter- 
mediate state . These terms thus give rise to "mixed" 
and "resonant" Raman scattering. 

The polarization dependence of Raman scattering can 
be generally classified using arguments of group theory. 
In essence, the charge density fluctuations brought about 
by light scattering are modulated in directions deter- 
mined by the polarizations of the incident and scattered 
photons. These density fluctuations thus have the sym- 
metry imposed on them by the way in which the light 
is oriented, and the charge density fluctuations obey the 
symmetry rules governing the scattering geometry. This 
is manifest in the dependence of the Raman matrix ele- 
ments on the initial and final fermion states. In general, 
the Raman matrix element Mpj = Mj^efe^ can be de- 
composed into basis functions of the irreducible point 
group of the crystal <I>^ (Devereaux, 1992; Hayes and 
Loudon, 2005; Klein and Dierker, 1984; Monicn and Za- 
wadowski, 1990; Shastry and Shraiman, 1991) 

MF,iiq^0) = J2M^^^, (14) 

with /i representing an irreducible representation of the 
point group of the crystal. Which set of contributes 
to the sum is determined by the orientation of incident 
and scattered polarization directions. As an example, if 
we consider the £>4h group of the tetragonal lattice, as in 
the cuprates, the decomposition can be written as 

Mfj = io^a)(er< + efeD 
+ ^0^(2)(efef) 

+ ^OB.,,(efe^-feref) 
+ ^OA..,«e^-eref) 
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FIG. 7 Schematic weighting of the Ught scattering transition 
for polarization orientations transforming as Big and -B29 for 
a Z)4h crystal. High symmetry points are indicated. Here a 
typical Fermi surface for optimaUy doped cuprates is repre- 
sented by the solid line, and the orientations of the incident 
and scattered polarization light vectors are shown with re- 
spect the to copper-oxygen bond directions. 



(15) 



with Op the corresponding projected operators and ef^ 
the light polarizations. This classification demonstrates 
that there is no mixing of representations for (7 = 0, i.e. 
the correlation functions read R ~ (O^O^/) = R^5^^^' 
and there are 6 independent correlation functions, each 
selected by combinations of polarization orientations. 

Following Shastry and Shraiman (1990), we list in Ta- 
ble I some common experimentally used polarization ge- 
ometries in relation to the elements of the transition rate 
R selected. One observes that a complete characteriza- 
tion of M can be made from a subset of the polarization 
orientations listed in the table. However, additional po- 
larization orientations can be useful to calibrate data and 
compare symmetry decompositions from different combi- 
nations of orientations. For geometries with polarizations 
in the (a — 6) plane in D4/J crystals, the irreducible repre- 
sentations cannot be accessed individually and must be 
separated by proper subtraction procedures. As a min- 
imum set, four independent configurations are required, 
while additional polarizations may be used for consis- 
tency checks. 



In addition, we have listed in Table I the representative 
basis functions <I>p (k) taken from the complete set of Bril- 
louin zone (BZ) harmonics for D^h space group (Allen, 
1976). This directly points out the connection between 
polarizations and the coupling of light to electrons. By 
virtue of the k-dependence of the light scattering transi- 
tion rate, excitations on certain regions of the BZ can be 
correspondingly projected out by orienting the incident 
and scattered light polarization vectors. Thus, Raman 
is one of the few spectroscopic multi-particle probes (the 
other being inelastic X-ray scattering) able to examine 
charge excitations in different regions of the BZ. 

For example, as demonstrated in Figure 7, for crossed 
polarizations transforming as -Big, light couples to charge 
excitations along the BZ axes {kx or fcy = 0), while for 
B2g, excitations along the BZ diagonals (fc^ = ±fcj,) 



are projected accordingly. Operators like OA^g can- 
not be accessed independently by linear polarizations 
alone. Only sums including circular polarizations al- 
low the isolation of components. Light scattering 
in this orientation can be coupled to chiral excitations. 
These important symmetry classifications have been ex- 
tremely useful to point out anisotropic electron dynamics 
in correlated insulators (Devereaux et al., 2003a, c; Shas- 
try and Shraiman, 1990), superconductors (Devereaux 
et ai, 1994a), disordered (Devereaux, 1992; Zawadowski 
and Cardona, 1990) and correlated metals (Einzel and 
Manske, 2004; Freericks and Devereaux, 2001; Freericks 
et ai, 2001). More recently they have been related to 
sum rules referring to BZ-projected potential energies in 
correlated systems (Freericks et ai, 2005). In the remain- 
ing part of this review, such symmetry classifications will 
be featured prominently. 



B. Formalism: Single-Particle Excitations and Weak 
Correlations 

We now consider specific cases where simplifications 
can be made to Eq. (13). First we assume that the inter- 
mediate many-particle states only differ from the initial 
and final states by single electron excitations. This is ex- 
act in the limit of non-interacting electrons, yet ignores 
the effects of many-body correlations, and specifically 
the role of Coulomb interactions on the reorganization 
of the initial state into the intermediate states by creat- 
ing many-particle excitations. We now focus on weakly 
interacting systems, where single-particle excitations are 
relatively well defined, and discuss strongly correlated 
systems in later sections. 



1. Particle-hole Excitations 

Eq. (13) can be simplified by replacing E,y in the de- 
nominators by Ei — EjB' + Ep and Ei — Ea' + Ea in the 
first and second terms, respectively, and use the closure 
relation — 1- Commutator algebra eliminates 

the four fermion matrix element, 



M;^,, = ^7„,^(F|4c^ |/), 

a,l3 

where 

7a,/3 = Pa,/3(q)ei • + — / 



(16) 



m yf\Efj - Efj' + hwi 



Efj — Epi — hujs 



(17) 



Specifying to states a, (3 indexed by momentum quantum 
numbers (such as Bloch electrons), from Eq. (11), the 
Raman response simplifies to a correlation function S of 
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Geometry 
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Basis functions <I'^i(k) 
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TABLE I Elements of the transition rate R for experimentally useful configurations of polarization orientations (given in Porto 
notation) along with the symmetry projections for the D4h point group relevant for the cupratcs. Here wo use notations in 
which X and y point in directions along the Cu-0 bonds in tetragonal cupratcs, while x and y arc directions rotated by 45°. L 
and R denote left and right circularly polarized light, respectively. In our convention left circular light has positive liclicity. (In 
a right-handed system the polarization rotates from x to y while the wavefront travels into positive z direction by A/4.) Note 
that in back-scattering configuration (see Fig. 2) with ei,s pinned to the coordinate system of the crystal axes the representation 
for incoming and outgoing photons with circular polarizations change sign in order to maintain the proper helicity. 



an effective charge density p, 

^-^nrl'^sici,in^n + io). (is) 

Here the Raman effective density-density correlation 
function is 



5(q,in)=^ 



,-0Ei 



dTe^"-(7|T,p(q,r)p(-q,0)|7), 



Tr is the complex time r ordering operator and 

^(q) = X]7(k,q)4+q,^Ck,a• 
k,c^ 



(19) 



(20) 



The scattering amplitude 7 is determined from the Ra- 
man matrix elements and incident/scattered light polar- 
ization vectors as 



7(k,q)=^7„,^(k,q)efef, 

a,/3 



(21) 



with 

7a,/3(k,q) = (5c«,/3 + ^Z] 



(k + q|pf |k.)(k. Ipflk) 



(k + q I j// I k,.)(k,. I j/J I k) 

-Ek-l-q — 7Sk„ — fl^s 



■ (22) 



Here, p"^ = pag±iqi,, r rpj^g effective dynamical density- 
density correlation function or Raman response S can be 
written in terms of a dynamical effective density suscep- 
tibiUty X via the fluctuation-dissipation theorem, 

Sid, n) = -7T-'{i + n{n, T)}x"(q, n), 



with n(fi, T) the Bose-Einstein distribution and 



X(q,n) = ([p(q),p(-q)])o. 



(23) 



Thus, for non-interacting electrons, the Raman response 
is given as a two-particle effective density correlation 
function and can be calculated easily using, e.g., diagram- 
matic techniques or via the kinetic equation (Devereaux 
and Einzel, 1995a). This reduces to evaluating the bubble 
diagram depicted in Figure 6 with vertices 7 depending 
upon the incident and scattered photon frequencies.^ 

The vertex 7 depends on polarization, but does not 
depend sensitively on q for g << fcp- This can be made 
more obvious if we consider the sum over intermediate 
states k^ in Eq. (22). The sum over intermediate states 
includes both the band index of the states created from 
the initial state (i.e. the conduction band) as well as 
intermediate states separated from the conduction band. 
The matrix elements of the former are proportional to the 
momentum transferred by the photon, which in the limit 
q « kp are smaller by a factor of (vp/c)^ than the other 
terms, with up (c) the Fermi (light) velocity, and can be 
neglected (Pines and Nozieres, 1966). For the remaining 
sum over the intermediate states separated from the con- 
duction band, we assume that Uis «\ Ei^^ — I and 



* We remark that technically these expression must be modified if 
a resonant condition is satisfied. In that case one needs to ex- 
pand to higher terms in the vector potential to capture resonance 
effects as perturbation theory breaks down. Yet for low energy 
Raman scattering, in most cases this is not crucially important 
since in real materials the intermediate states reached via a direct 
resonance are quite broadened by interactions and the resonant 
terms are not orders of magnitude larger than the non-resonant 
terms. Thus this treatment may be of more general utility. 
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recover the widely used effective mass approximation^ 

1 d^Ei, 



7a,/3(k,g ^ 0) = 



(24) 



The symmetry classifications listed in Table I thus can 
connect excitations created by certain polarization orien- 
tations to properties of the band structure. Yet it must 
be kept in mind that this connection can only be made 
in the limit of small oJi « . 



2. Im(l/e) and Sum Rules 

We consider first the case when the incident polariza- 
tion is parallel to the scattered polarization and the band 
Ek is isotropic and parabolic, as for the electron gas or 
lightly-filled isotropic band metal. The scattering ampli- 
tude 7 is then independent of k and the effective density 
p(q) is simply proportional to the pure charge density 
p. Using the definition of the complex dielectric function 
(Mahan, 2000; Pines and Nozieres, 1966) 



(25) 



is obtained with Vq the bare Coulomb interaction, and 
Xsc is the screened or irreducible polarizability. It is de- 
termined from the full polarizability x via 



Xsc(q,^) = 



X(q,f^) 



1 - i'qx(q, ^) ' 



(26) 



but is most easily identified diagrammatically as all 
contributions to the polarizability which are irreducible 
with respect to the interaction. The dynamical density- 
density correlation function (or structure factor) follows 



as 



5(q, O) = [1 -I- n{n)] Im 



TTVa 



1 



e(q,r!) 



(27) 



The Raman response S is proportional to S with a con- 
stant of proportionality determined by the ratio of effec- 
tive to free electron mass. 

Inelastic light scattering occurs via the creation of 
charge fluctuations inside the unit cell which are cou- 
pled via the Coulomb interaction to charge fluctuations in 
other unit cells. These intercell excitations are therefore 
well screened by the Coiilomb interaction and reduce the 
scattering cross section at small q. In particular, for small 
q, xiq, ^) ^ NpVpq"^ /Q'^ , and the response is governed by 
the plasma frequency. The Raman response thus obeys 
the longitudinal sum rule resulting from particle-number 
conservation (Pines and Nozieres, 1966), 



dfl il Im 



1 



e(q,n). 



2"P^ = ^;^' ^''^ 



where Qpi is the plasma frequency and A'^ is the number 
of electrons of mass m and charge — e in the system. Thus 
the only contribution for q ^ comes from exciting the 
plasmon being the only charge excitation available for 
light scattering at small g in a free-electron gas. 



3. Intra- vs. Inter-cell Charge Fluctuations 

In the more general case of light scattering in solids, 
the Raman response may have other contributions com- 
ing from intra-cell charge fluctuations, provided the band 
structure is non-parabolic, as pointed out by Platz- 
man (1965) and Wolff (1968). Then light can create 
anisotropic charge fluctuations which are zero on aver- 
age inside the unit cell and thus are not screened via 
the long-range Coulomb interaction, as pointed out by 
Abrikosov and Genkin (1973). The general expression 
for the screened Raman response function x^'i-y can be 
written as (Devereaux and Einzel, 1995a; Monien and 
Zawadowski, 1990) 



x: 



7.7 



X7,7 ■ 



X7,lXl,7 I X7ilXl,7 ^, 

T o XsC) 



Xi,i 



(29) 



^ This is derived in Appendix E of Ashcroft and Mermin (1976). 



where Xsc = Xi,i(l ~ I'qXi.i)"^- This is an exact expres- 
sion, where the subscript 7 denotes the effective Raman 
density and 1 denotes the pure charge density, obtained 
when the usually momentum dependent vertex 7 is re- 
placed by a constant. The respective x-s describe the 
full density-density, density-Raman density, and Raman 
density-Raman density susceptibilities which are again 
each irreducible with respect to the interaction. 

The first term in Eq. (29) is the bare response for a 
neutral system, and the other terms represent the back- 
flow needed to enforce particle number conservation of 
charge density fluctuations and gauge invarianee. These 
terms arc important for light scattering conflgurations 
which transform according to the symmetry of the lat- 
tice, such as Aig in _D4h crystals. In particular, if we 
consider scattering from pure charge density fluctuations 
where 7 is a constant independent of momentum, which 
is an Aig representation, the first two terms in Eq. (29) 
cancel and Xsc and Eq. (27) are recovered. This is in- 
escapable for g = 0, since then the scattering operator 
is given in terms of the total density of electrons, which 
commutes with the bare Hamiltonian and therefore can- 
not give inelastic scattering channels to light. On the 
other hand, if we consider the scattering vertex 7 to de- 
pend on wavevector, the backfiow terms are not capable 
of completely canceling the bare Raman response. 

Momentmn dependence of the vertex 7 is quite gen- 
eral for electrons in solids. In particular, for crossed 
light polarizations projecting out representations of lower 
symmetry than that of the lattice x-y.i is identically zero 
by symmetry for g = and the backfiow terms make 
no correction to the Raman cross-section. This occurs 
for Big , B2g , and Eg scattering geometries in D^i^ sys- 
tems such as the cuprates, for example. While there is 
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no conservation law for light scattering from the exci- 
tations created by crossed polarizations (Kosztin, 1991), 
there are sum rules which relate the Raman intensity to 
model-dependent potential energies projected in different 
regions of the BZ (Freericks et ai, 2005). 

C. Formalism: Strong Correlations 

1. General Approach to Treating Correlations 

Section II. B outlined a general approach to inelastic 
light scattering when the intermediate states differ from 
the initial and final states only by individual single elec- 
tron energies. This holds in the limit of weakly interact- 
ing electrons. In this section we show that the formal- 
ism is also valid in the Hcisenbcrg limit of the Hubbard 
model including the manifold of zero and 1 doubly oc- 
cupied sites. These are two limits in which either the 
kinetic energy or the potential energy of the electrons is 
dominant. However, the more general case of interest to 
most systems is tackling the problem when both kinetic 
and potential energies are roughly equal. The interac- 
tions are sufficient to give broad spectral functions as 
measured by ARPES, where the incoherent part of the 
spectral function is manifest from the strong many-body 
interactions mixing individual electron states. 

In this situation, one must resort back to Eq. (13), and 
correlation functions involving two, three and four par- 
ticles are needed, as depicted in Figure 6. Yet usually 
one is not interested in treating many-body correlations 
over various bands and puts focus on a few bands close to 
the Fermi level having strong correlations. For example 
in the cuprates, one usually takes downfolded Hamilto- 
nians involving only Cu 3dx2_y2 and O 2px,y orbitals, 
with short-range Coulomb interactions for two electrons 
in on-site or neighboring orbital states. The downfolding 
procedure, such as that described by Lowdin (1951), re- 
moves all other bands (effectively moving them infinitely 
far away in energy from the focus bands), such as the 
apical oxygen, Cu 4s and other Cu d orbitals, and treats 
the electrons in the bands of interest as having renormal- 
ized energy dispersion. In the cuprates the downfolding 
results in either a few bands from local density approxi- 
mation (LDA) approaches in one case or cell perturbation 
theory in the strongly interacting case. Thus we consider 
a downfolded tight-binding Hamiltonian 

Ho = J2 *i'oci<TCj,a + Hint, (30) 

where tij are the effective hopping integrals resulting 
from the downfolding procedure, and Hj^t describes the 
relevant interactions. For example, we may consider a 
square lattice of electrons with strong on-site repulsion U 
much greater than the electron hopping t in the Hubbard 
model, 

H = -tY^ 4,<TCj> + t^Xni,|ni,4, (31) 



where denotes a sum over nearest neighbors. 

The electronic eigenstates fall into two bands for large 
U at 1/2 filling - e.g., the occupied lower and unoccupied 
upper Hubbard bands - separated by an energy U for 
double occupancies. Away from 1 /2 filling, quasiparticles 
develop. The microscopic Hamiltonians can be viewed 
as families of models related to the Hubbard model, 
such as the Falicov-Kimball model, Anderson model, or 
Anderson-Fano model, without loss of generality. 

In essence, the sums over intermediate states in 
Eq. (13) are separated into groups of bands lying far away 
in energy from the initial and final states (i.e. the bands 
projected out) and bands lying nearby the initial and 
final states (considered bands with correlations). The 
former grouping of intermediate states are considered as 
in Section II. B to be approximated by the effective mass 
contribution, Eq. (24). The remaining terms involve ma- 
trix elements of the current operator between the remain- 
ing band of interest^°. 

The interaction of light with these downfolded elec- 
trons can be treated via the Peierls construction, in which 
the creation and annihilation operators develop a phase, 

Q,. ^c,,„e-^^/-»^-<^^. (32) 

The resulting scattering Hamiltonian obtained by ex- 
panding in powers of A reads 

2 

Hint = • A + ^ Yl ^«7a/3^/3, (33) 

where 

= E ^^-('^ + - (34) 

k 

is a component of the current operator j and 

^"'^^'^^ = 5 afc^k^"^^"" + "^/^^'"^^ ~ "1/^^ ^^^^ 

is the stress tensor operator. Both operators are thus 
formed from the energy dispersion of the downfolded 
band structure. The matrix element can be written in 
compact form: 

Mf,/(q)-E<e^^"'''(q)' 

a,f3 

M«/3(q) = (F|7,,^(q)|7) (36) 



We remark that this is not exact as it inaccurately treats resonant 
scattering processes occurring within the conduction band. Since 
these processes are usually taken to be damped, even though this 
approach is strictly speaking not exact, it provides a good start- 
ing point for considering Raman scattering in correlated single- 
band systems, and as such, has been widely used (Shastry and 
Shraiman, 1991). 
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with the sum over intermediate states v of the Hamilto- 
nian Eq. (30). The Raman cross section can be separated 
into non-resonant, mixed, and resonant contributions: 

R{n) = RN{n) + iiM(O) + (37) 

where the nonresonant contribution is 

RN{n) = ^ ^^P(^/^-^^) Y^i^ s{Ef -Ej- Ml), 

(38) 



the mixed contribution is 

exp{-f}Ei) 
Z 

I,F,v 



-As) -At) 

■iF.u-'u.I 



Rm{^) = 



Ai) . . 



^^'^ \ - Ei - hvi E^- Ei + hws 



~Ai) -As) 

■ll,v-lv.F 



-As)-A^) 

JI,vJv,F 



\ Ely — El — fkui Ely — El + hojs 
X SiEp -E,y-Ml), 

and the resonant contribution is 
^ exp(-/3£j) 

I,F,uv' 




Rum 



:-w -As) 

JI,vJv,F 



As)A^ 

■II,vlv,F 



E,y ~ El — hjji E,y - El + huJs 



-As) -Ai) 

JF.u'Jv' J 



-S) -As) 

Jf,v'Jv'I 



\ E,ji — Ej — hcoi Eiyi — Ej + HlOs I 
X 5{Ef - El - Ml). (40) 

We have introduced the symbols 

af3 a 

(41) 

and denote Ok,a as the matrix element {k\0\X). Some of 
the contributions are depicted diagrammatically in Fig- 
ures 8, 9, and 10 for the non-resonant, the mixed, and 
the resonant terms, respectively. In the limit D ^ oo 
these are the main diagrams to consider. Additional di- 
agrams involving multi-particle vertex renormalizations 
generally contribute for finite dimensions (not shown). 

In general, the matrix elements that enter into 
Eqs. (38-40) are not easy to calculate for an interacting 
system, so the summations are problematic to evaluate. 



In particular, one needs to evaluate the irreducible stress 
and current vertices, depicted in Figures 8-10 by the 
hatched symbols. Contributions to these vertex dressings 
include many-particle renormalizations. A particularly 
complicated one is shown in Figure 11 which represents 
4-particle vertex corrections. Moreover, anal}d;ic contin- 
uation must be performed to obtain the Raman response 
i?(ri) on the real axis from the imaginary axis. While this 
is relatively straightforward for the non-resonant case, 
mixed and resonant cases are problematic because of 
the complicated dependences on each of the frequencies 
LUi^s,Q which have to be analytically continued. While 
the continuation has been worked out recently for the 
general case, evaluating these diagrams for general inter- 
actions has proved elusive. 

The overall complexity of the problem limits the evalu- 
ation of the light scattering cross section to generic inter- 
acting systems. Only recently, these diagrams have been 
evaluated exactly in a DMFT treatment of the Falicov- 
Kimball model (Shvaika et al., 2004, 2005). 




FIG. 8 Fcynmari diagrams for non-rosonant Raman scatter- 
ing. The wavy and the solid lines denote photon and electron 
propagators, respectively. The cross-hatched rectangle is the 
reducible charge vertex. The symbol 7 denotes the stress- 
tensor vertex of the corresponding electron-photon interac- 
tion. Prom Shvaika et al. (2005). 




FIG. 9 Fcynmau diagrams for the mixed contributions to Ra- 
man scattering. The symbols jf and ji remind us to include 
the relevant vertex factors from the current operator in the 
electron-photon interaction. Prom Shvaika et al. (2005). 
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dimensional Hubbard model. Following Shastry and 
Shraiman (1991) we consider a system with N interacting 
electrons in which the manifold of states can be classified 
by the number n of doubly occupied sites. The well- 
known Heisenberg Hamiltonian emerges from projecting 
the Hubbard model down onto the reduced Hilbert space 
containing no double occupancies: 



-f^Heisenberg — J ^ ^ Sj • 



(42) 
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FIG. 10 Feynman diagrams for the resonant contributions to 
Raman scattering. From Shvaika et al. (2005). 

2. Correlated Insulators - Heisenberg Limit 

To emphasize the generality of Eq. (13) we now con- 
sider the large U limit for the insulating half-filled two- 



with J = At^ /U the Heisenberg exchange constant. 
Higher manifolds containing empty holes and doubly oc- 
cupied states can be labeled according to the net spin 
configuration {cr} of the N — 2n singly occupied sites, as 
well as the locations {R} of the empty rhoie and dou- 
bly occupied rdoubie sites. We denote these states as 
I n; {a}, {R}). These states are connected to each other 
via 



1; {u'y, {R}) = 4(rdoubie)ca(rhoie) |0; {u}), (43) 
IlrC^^ (r) I vac) and | vac) denotes the 



where | 0; {cr}) 
vacuum. 



Light scattering thus occurs via transitions out of the 
manifold of singly-occupied states. To leading order for 
large U, only the n = and n = 1 manifold of states 
contributes to light scattering via Eq. (13), with n = 
denoting the ground state and n = 1 the manifold of 
intermediate states having one doubly occupied and one 
empty site. The first term containing ma,/3 cannot con- 
tribute for the half-filled lattice, and thus only interband 
scattering between the upper and lower Hubbard bands 
occurs via the p • A term. The energy difference between 
these excitations is U to lowest order in i/i7, allowing us 
to write the matrix element Eq. (13) in the form 



i 



" k' 



k+q a ^ " k'+q 



FIG. 11 Feynman diagrams for a typical parquet-like renor- 
malization. This resonant diagram has a simultaneous hor- 
izontal and vertical renormalization by the two-particle re- 
ducible charge vertex. From Shvaika et al. (2005). 
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Mf,i = -t^-f } I ^'^ W e« • <5 I 1; I • 5'ji(r') I 0, {of}) 



1 



1 



U -hcui U + fiLOs 



with the current operator defined as 



ii,s(r) = it[cl{r + 5a ■ ei,s)ca{r) - 4(r)c<^(r + 5a ■ ei,^)] 



(44) 



(45) 



Here 5 is a unit vector connecting a site with its nearest neighbors. The intermediate states u represent a sum over 
spin configurations and locations of both the doubly occupied and the hole sites. Substituting Eq. (43) into Eq. 
(44) collapses the intermediate state sum, leaving four terms connecting initial and final states. Using the identity 
1/2 — 2Si • Sj = c^r + 5a)c{r)c^r)c{r + 5a) valid in the manifold of singly occupied states, one obtains the light 
scattering Hamiltonian of Elliot and Loudon (1963) and Fleury and Loudon (1968), 



^^EFL = Sr • Sr+<5a(es ' 5){ei ■ 5) 

r,S 



1 



1 



u — huji u + hujs 



(46) 



We note that the polarization dependence is crucial as 
well. For xx + yy polarizations projecting the fully 
symmetric components, the light scattering Hamiltonian 
commutes with the nearest neighbor Heisenberg Hamil- 
tonian and thus does not give inelastic scattering in the 
Aig channel. Moreover, i?2g {xy) is also identically zero. 
As a result, a large signal appears only in the Big chan- 
nel {xx — yy). These restrictions are lifted however if 
longer range spin interactions are considered (Shastry 
and Shraiman, 1991). 

The collapse of the intermediate states allowed us to 
replace the operators with projected spin operators con- 
fined to the restricted Hilbert space of the n = 0, 1 man- 
ifolds. Thus, in this limited Hilbert space, the formalism 
is similar to non-interacting electrons in that the opera- 
tors appearing in the scattering matrix may be simplified. 
If the Hilbert space is enlarged to include larger mani- 
folds, then this would no longer be the case, and thus 
including terms to higher order in t/U becomes highly 
non-trivial and is still one of the challenges to merge a 
weakly interacting picture into a strongly interacting one. 

We note that Equation (46) was derived effectively as 
an expansion in t/(C/ — JvjJi). Therefore, the scattering 
Hamiltonian is limited to cases when both the number 
of holes and double occupied sites arc restricted and off- 
resonance conditions apply, with the incident photon en- 
ergy huji far away from U. Efforts to extend the treat- 
ment to more general conditions involve understanding 
the motion of holes or doubly occupied sites in an ar- 
bitrary spin background. This has proved to be a hard 
task. 



now specify the electronic states from which light can be 
scattered and review the various theoretical treatments 
for specific models of interacting electrons. Emphasis is 
placed upon how symmetry can be used to highlight elec- 
tron dynamics on regions of the BZ, and general features 
for each model system will be presented. We first con- 
sider the case where the correlations among electrons are 
weak. 



1. Weakly-Interacting Electrons 

In this subsection, we consider electrons as having very 
well defined eigenstates labeled by energy and momen- 
tum and having sharp spectral functions. Apart from 
the form of the energy dispersion ^k, the results are 
rather general and governed largely by phase space con- 
siderations. We must, however, consider the long-range 
Coulomb interaction in order to account for charge back- 
flow and screening, and wc utilize the results derived in 
Section II. B and the general expression Eq. (29). 

The use of Xsc rather than x takes into account the 
most drastic manifestation of the long-range Coulomb 
interaction, viz. screening. For weakly interacting elec- 
trons, the Random Phase Approximation (RPA) is ac- 
ceptable, which replaces Xsc by the Lindhard function 
for non- interacting electrons. The response functions are 
determined by the Lindhard kernel, 



2 



/(a)-/(^k+q) 

'^k - ^k-Fq + W.- iS' 



(47) 



D. Electronic Charge Relaxation 

In Sections II. B and II. C we reviewed the general for- 
malism of the theory of Raman scattering for weakly 
and strongly correlated systems. In this subsection we 



for general vertices a,b appearing in Eq. (29). In a free 
electron gas, the Raman response is given by Eq. (27) 
or equivalently, the last term in Eq. (29). For large q, 
collective excitations are unimportant and light scatter- 
ing occurs via creation of particle-hole excitations in the 
Landau continuum. However, as the only phase space 
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for creating particle-hole pairs comes from finite q trans- 
ferred from the photons, the resulting response is a con- 
tinuum varying linearly with CI at small frequencies and 
extending up to a cut-off Qc — v-pq from the borders of 
the continuum (Mahan, 2000). The low-energy intensity 
is proportional to q^, and the only excitation left at q = 
is the collective plasmon. The Raman response for the 
free electron gas is shown in Figure 12. 

For electrons in a solid, however, the non-parabolicity 
of the energy dispersion results in charge fluctuations 
which are anisotropic in the small q limit and thus can 
survive screening and give more weight at low energy 
transfers. Formally, an additional contribution to the re- 
sponse is given by the first two terms in Eq. (29). Yet 
phase space restrictions still produce an inescapable cut- 
off at V,c = vpq (Wolff, 1968), and the response resembles 
that shown in Figure 12. This is also the case if scatter- 
ing occurs for the 2 dimensional electron-gas (2DEG) in 
the absence of a magnetic field (Jain and Das Sarma, 
1987; Mishchenko, 1999) or for complex Fermi surfaces 
(Ipatova et ai, 1983). An RPA treatment for resonant 
scattering has been given by Wang and Das Sarma (1999, 
2002). 

Recently, substantial progress has been made in under- 
standing the Raman response in the integer or fractional 
quantum Hall regimes of the 2DEG. Space limitations do 
not allow us to review these systems; so for brevity, we 
cite only a recent reference (Richards, 2000). 



2. Impurities 

Excitations at low energies in non-interacting elec- 
tronic systems can arise for small q via electronic scat- 
tering from impurities, where momentum contributed by 
impurity scattering can provide phase space for electron- 
hole creation which is anisotropic in the Brillouin zone. 
For example, if one considers a general electron-impurity 



interaction of the form 



KT < qVp 

KT«Ep 




FIG. 12 Raman response of an electron gas. From Platzman 
(1965). 



^imp = ^ VkM'c]^ a^k',cr, 



E 

k,k',(T 



(48) 



with an anisotropic interaction Vk,k', the gauge invari- 
ant Raman response is given via the diagrams presented 
in Figure 8. If we make a symmetry decomposition of 
the scattering amplitude 7(k) = J2l'^l^lO^) in terms 
of basis functions <i>L of the Brillouin zone, the result- 
ing response in channel L corresponding to a particu- 
lar light polarization orientation has a Drude Lorentzian 
form (Devereaux, 1992; Fal'kovskii, 1989; Zawadowski 
and Cardona, 1990): 



l + (f7r£)2' 



(49) 



with A'f the density of states at the Fermi level. 1/t£ = 
1/t -f Dq^ — \/tl is the effective scattering rate where 



1/r = l/Ti=o 



S 



Vk,k' 



(50) 



involving an integration of the Fermi surface S'k normal- 
ized to the Fermi area S. Here is the impurity con- 
centration, D = -^VpT is the diffusion constant. The 
anisotropy of the impurity scattering is characterized via 
orthonormal basis functions $l 



Vk.i 



L,L' 



$2(k)$L(k')Fj 



L,L' 



(51) 



using an intelligent basis where the interaction is diagonal 
Vl,l' = Sl.l'Vl- Then, 1/tl = 2nn,NpVL and 1/tl=o 
is the dominant contribution. 

The resulting Raman spectrum (Figure 13) grows lin- 
early with frequency fi, decays as l/fi, and has a peak 
when flrl equals 1. The width of the Lorentzian re- 
flects the rate for which charge density excitations hav- 
ing symmetry L decay into all other channels. The light 
polarizations select the type of excitation L created, and 
thus allow a way to probe the anisotropy of the impurity- 
electron interaction. The decay of the charge density 
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FIG. 13 Raman response from impurity scattering in an oth- 
erwise non-interacting system. From Zawadowski and Car- 
dona (1990). 
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fluctuations can occur via finite q through the diffusion 
term and all contributions other than Vl which relax elec- 
trons out of the state L. Put another way, Dq^ + l/r are 
"scattering out" processes, while I/tl is a "scattering in" 
process, giving an effective scattering rate 1/t£. This is 
a consequence of a gauge invariant treatment including 
charge backflow (Coulomb interaction) as well as density 
preserving scattering (impurity vertex corrections). 

In the limit of weak scattering, the response collapses 
into a delta-function, reflecting momentum conservation. 
One can note the obvious connection of the response plot- 
ted in Figure 13 to the Drude conductivity, although even 
for simple impurity scattering the two response functions 
are not related by a power of frequency as soon as the im- 
purity potential has any momentum anisotropy. Other- 
wise, for purely isotropic impurity scattering the conduc- 
tivity and the Raman response are related by a power of 
frequency - the so-called Shraiman-Shastry relation given 
in Eq. (6) (Freericks and Devereaux, 2001; Shastry and 
Shraiman, 1990). 

3. Interacting Electrons - Non-Resonant Response 

However, the most important application of light scat- 
tering is for systems where the electronic correlations are 
strong and cannot be treated in standard RPA. Thus, 
while long-range Coulomb screening is still important in 
order to maintain gauge invariance, the interactions in- 
troduce generally complex dynamics in specific regions of 
the BZ. In this case, the electron self energy S, as well as 



the vertex corrections to the light scattering amplitude 7, 
depend normally on both momentum and energy, mak- 
ing the light scattering evaluation more difficult. On the 
other hand, anisotropies of the electron dynamics can be 
explored. 

Here we start by considering non-resonant scattering, 
since this is an area in which by far most theoretical 
treatments lie, as it is simpler to evaluate than the 
mixed or resonant terms. We note that many calcula- 
tions of /m(l/e) have been performed from ab-initio ap- 
proaches in the context of inelastic X-ray scattering (see 
Gurtubay et al. (2004) and references therein for recent 
work). There, the focus is largely on the q— dependence 
of the response, and electron-electron interactions have 
been treated in various ways (Ku et a/., 2002). Yet, to 
our knowledge, no calculation exists for Raman scatter- 
ing in a simple Fermi liquid in which inelastic scattering 
processes via the Coulomb interaction are incorporated 
exactly, although recently dynamical mean field theory 
(DMFT) in correlated metals has been used (Freericks 
and Devereaux, 2001). This is because the irreducible 
charge vertex is not generally known in models with 
strong correlations, with the exception of the Falicov- 
Kimball model. Thus we focus more on the polarization 
dependence and investigate contributions to Raman scat- 
tering from non-conserved charge fluctuations. 

The general expression for the two-particle correlation 
function describing the non-resonant Raman response 
reads 



= ^,i^) = -yQY.Y. 7(k)G(k, ioj)G{h, loj + if2)r(k; iw; i^). (52) 

iu! k 

Similar expressions are obtained for x-y,! and xi,i where the vertices 7 are successively replaced by 1 to be inserted 
into Eq. (29), or may be generally represented in terms of the reducible Raman vertex as shown in Figure 8. In the 
ladder approximation, the renormalized vertex is given by a Bethe-Salpeter equation: 

r(k; iu; in) = 7(k) + y^^^Vi^- k', ioJ - iu;')G(k', iw')G(k', iu' + ifi)r(k'; iu'; ifl). (53) 

iuj' k' 



Here, V{k,u)) is the generalized electron-electron inter- 
action, and we have suppressed spin notation. If one 
neglects vertex corrections such that the theory is not 
gauge invariant, the Raman response has a particularly 
simple form given by Eq. (3). The effect of the long-range 
Coulomb interaction is treated formally in the same way 
as in Eq. (29), with the vertices replaced by the normal- 
ized vertex as a solution to the Bethe-Salpeter Eq. (53). 

Eqs. (52) and (53) have been the starting point for 

many studies of light scattering treating electron-electron 
interactions in effective models. These include systems 



which have nearly nested Fermi surface segments (Vi- 
rosztek and Ruvalds, 1991, 1992) or antiferromagnetic 
spin fluctuations (Devereaux and Kampf, 1999; Kampf 
and Brenig, 1992). Similarly, a slave boson approach 
to the t — J model (Bang, 1993), electron-phonon in- 
teractions (Itai, 1992; Kostur, 1991, 1992; Rashkecv and 
Wendin, 1993) and fluctuation-exchange (FLEX) treat- 
ments of the Hubbard model (Dahm et al., 1999) have 
been considered. While these studies involve approxi- 
mate solutions, more recently the use of DMFT has pro- 
vided exact results in the limit of strictly local correla- 
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tions in the Hubbard (Freericks et al, 2001, 2003a) and 
Falicov-Kimball models (Freericks and Devereaux, 2001). 

Two aspects of the Raman response are generally in 
the main focus: the frequency dependence of the broad 
continuum extending well past qvp, and the polarization 
dependence. We discuss first the spectral response. 

In the context of the cuprates, Varma and cowork- 
ers pointed out that a flat, nearly frequency indepen- 
dent response could be obtained if the imaginary part of 
the electron self-energy depended linearly on frequency 
(Varma et al., 1989a; Varma, 1989b). The response is 
then given in terms of the scale- invariant ratio of fku/k^T 
and approaches a constant at large frequency transfers. 
This can be understood phenomenologically by replac- 
ing l/r£ with 1/T£(0,r) ac ma.x{kBT,hn) in Eq. (49)." 
A scale-invariant response at low frequencies is a gen- 
eral consequence of systems in proximity to a quantum- 
critical point, but this scale invariance is broken outside 
the quantum critical regime. "Marginal Fermi liquid" 
behavior emerges for instance when scattering is consid- 
ered in a nested Fermi liquid (Virosztek and Ruvalds, 
1991, 1992), low Fermi-energy systems (Dahm et al., 
1999; Devereaux and Kampf, 1999), and slave-boson sys- 
tems (Bang, 1993). A broad background very similar 
to "marginal" behavior is also found for strongly coupled 
elcctron-phonon systems (Itai, 1992; Kostur, 1991, 1992). 
As a representative example, we show the response calcu- 
lated by Virosztek and Ruvalds (1992) for a nested Fermi 
liquid in Figure 14. 

Low energy electron dynamics can be extracted by 
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FIG. 14 Raman response as a function of temperature ob- 
tained by Virosztek and Ruvalds (1992) for a system with a 
nested Fermi surface. The upturn at low frequencies is to the 
(log(f2)) frequency dependence of the effective mass in this 
model. 



studying the Raman response in the limit ^ 0. Ne- 
glecting vertex corrections, the low frequency response 
roads (Devereaux and Kampf, 1999; Venturini et al., 
2002a) 



2S(:(C,T) 



(54) 



Here, Np is the density of electronic levels at the 

Fermi energy Ep, SJ^ is the imaginary part of the 
single-particle self energy related to the electron life- 
time as h/2Y.'l{Lo,T) = Ti^{oj,T), Zi,{u},T) = (1 - 
d^[^{Lu.T)/duj)^^ is the quasiparticlc residue, is the 
equilibrium Fermi distribution function, and (• • •) de- 
notes an average over the Fermi surface. Thus the inverse 
of the Raman slope 



(55) 



measures the effective "scattering rate" of the quasipar- 
ticles in a correlated metal, and can be best thought of 
as a "Raman resistivity" . 

In systems with isotropic interactions, the polarization 
dependence drops out and the slope of the low-frequency 
Raman response; is given in terms the low energy quasi- 
particlc scattering lifetime, r^(T) oc h/T{T), as an ex- 
tension of Eq. (49). Yet, in strongly correlated systems, 
the quasiparticlc residue Z and, importantly, vertex cor- 
rections, enter as well. In a correlated or a strongly disor- 
dered metal (near an Anderson transition, e.g.), however, 
a finite energy might be necessary to move an electron 
from one site to another one. Thus, in spite of a non- 
vanishing density of states at the Fermi level, as observed 
in an ARPES experiment for instance, no current can be 



1000 



Q. 

o 



c 
o 
£ 
o 



> 



100 




0.2 0.4 0.6 0.8 

Temperature 



This is only an approximation since 1/T£(f2,r) depends now on 
energy. Causality requires that the relaxation function has real 
and imaginary part, M{Q,,T) = f2A + i/r£ (Gotze and Wolfle, 
1972; Opel et al, 2000). 



FIG. 15 Inverse Raman slope (see Eq. (55)) close to a metal- 
insulator transition at a value U = \/2 in the Falicov-Kimball 
model for D = oo (Freericks and Devereaux, 2001). All ener- 
gies and temperatures are measured in terms of the hopping 
t. 
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transported and r^(T) ^ h/T{T). This is an important 
difference between single- and two-particle properties. 

Figure 15 displays the inverse Raman slope defined in 
Eq. (55), as determined via a DMFT treatment of the 
Falicov-Kimball model in the vicinity of a metal-insulator 
transition, as a function of the Coulomb repulsion U 
(Freericks and Devereaux, 2001). It provides an illus- 
trative example of how Fermi-liquid-like features evolve 
as the lifetime of putative quasiparticles increases due to 
decreased role of correlations. For small U, the corre- 
lated metal displays an inverse slope oc as a canonical 
Fermi liquid in the metallic state. A pseudogap open- 
ing in the density of states with increasing U drives the 
inverse slope into insulating behavior, increasing as the 
temperature decreases. 

As a second important application Eq. (55) can illu- 
minate the anisotropy of electron dynamics due to the 
momentum-dependent weighting factors of the polariza- 
tion orientations and self energies. The geometry of light 
scattering orientations, as given by the form factors listed 
in Table I, project out the ratio of quasiparticle residues 
and scattering rates in different regions of the BZ. As a 
consequence, the Raman spectra show polarization de- 
pendent behavior determined largely by the self energies 
and vertex corrections near the regions projected by the 
scattering vertices 7. Figure 16 plots the Big and B2g 
Raman response calculated in a spin-fermion model in 
which electron scattering is most pronounced involving 
antiferromagnetic reciprocal lattice momentum transfers 
Q = (tt, tt), leading to "hot" quasiparticles near the BZ 
axes (projected by Big form factors) and "cold" quasi- 
particles along the BZ diagonals (projected by B2g form 
factors). Therefore, the Raman response has a sharp 
quasiparticle peak for i?2g scattering at low energies due 
to the long quasiparticle lifetimes, while the response in 
Big is dominated by strong incoherent scattering leading 
to a suppression of the quasiparticle peak at low energies 
and an essentially structureless continuum. 

Lastly, we note that the non-resonant Raman re- 
sponse has also been calculated for exchange of fluctu- 
ation modes at wavevectors Q and -Q for systems near a 
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FIG. 16 Polarization-dependent Raman response in a spin- 
fermion model for a fixed temperature for three different val- 
ues of the coupling constant. From Devereaux and Kampf 
(1999). 



spin-density wave instability (Brenig and Monien, 1992; 
Kampf and Brenig, 1992; Venturini et ai, 2000) and a 
charge-density wave instability (Caprara et ai, 2005). 
Here, the Raman response is sensitive to the light po- 
larizations and has a peak centered at twice the energy 
of the fluctuating mode. 



4. Interacting Electrons - Resonant Response 

In addition to the non-resonant response, one has the 
mixed and resonant contributions to c;onsider. Typi- 
cally these diagrams are neglected in the weak correlation 
limit, as they can be summed into two-particle response 
functions as discussed in Section II. B. In the insulating 
case, only the resonant terms are kept, as the studies fo- 
cus on excitations across a charge-transfer or Hubbard 
gap. This has been calculated in systems exhibiting ID- 
Luttinger behavior (Kramer and Sassetti, 2000; Sassetti 
and Kramer, 1998; Sassetti et al., 1999; Wang et ai, 2004) 
where bosonization techniques can be applied. Yet, gen- 
erally treating all diagrams on equal footing is technically 
demanding. Only recently, an exact evaluation in DMFT 
has been performed (Shvaika et al., 2004, 2005). 

It is well known that many of the Raman signals in 
correlated metals and insulators display complicated de- 
pendences on incoming photon frequency uji. For exam- 
ple, the Big two-magnon feature at roughly 350 meV in 
the thoroughly studied insulating parent cuprates has a 
resonance for incident photon energies near 3 eV. As a re- 
action to the experimental results in the cuprates, much 
theoretical work has been devoted to Raman scattering 
in a two-dimensional Heisenberg antiferromagnet using 
the Fleury-Loudon model Eq. (46). 

In the nearest neighbor Heisenberg antiferromagnet, 
one treats the spin operators using a Dyson-Maleev rep- 
resentation of magnons with dispersion 2J. Due to the 
Big form factor, in the absence of magnon-magnon inter- 
actions at T=0, a sharp peak at 4 J appears as the top of 
the magnon dispersion is projected out, and the response 
abruptly drops to zero for large Raman shifts (Sandvik 
et al., 1998). However, since the light scattering is local- 
ized to neighboring spins, magnon-magnon interactions 
must be included, and the peak becomes more symmetric 
and shifts to ^ 3 J via breaking exchange bonds between 
local neighbors, as shown in Figure 17. 

There have been many developments on the Fleury- 
Loudon model, which has been addressed via critical 
fluctuation analysis (Halley, 1978), series expansions 
(Singh et al, 1989), lower (Morr and Chubukov, 1997; 
Parkinson, 1969) and higher (Canali and Girvin, 1992; 
Chubukov and Frenkel, 1995b) order spin wave theories, 

J studies at flnite doping (Prelovsek and Jakhc, 1998), 
exact diagonalizations of small clusters (Tohyama et al., 
2002), excitonic cluster approaches (Hanamura et al., 
2000), flnite temperature Quantum Monte Carlo meth- 
ods (Sandvik et al., 1998), and studies of bilayer effects 
(Morr et al., 1996), 2-leg spin ladders (Jurecka et al.. 
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2001) and ring exchange (Katanin and Kampf, 2003), 
giving a thorough treatment of two magnon scattering 
from spin degrees of freedom in the non-resonant regime. 
Scattering from channels other than Big, and describing 
the anisotropic hneshape of the response, have been ad- 
dressed via longer range spin exchange interactions and 
by exact diagonalizations of magnons coupled to phonons 
(Freitas and Singh, 2000), using an earlier approach of 
Lorenzana and Sawatzky (1995). Lastly, recent develop- 
ments concern scattering from orbiton degrees of freedom 
(Okamoto et ai, 2002) and scattering within a resonant 
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FIG. 17 Cartoon of the two-magnon scattering process in a 
2D Heisenberg antiferromagnet. An incident photon causes 
an electron with spin a to hop leaving a hole and creating a 
double occupancy in the intermediate state with energy oc U. 
One particle of the double with spin —a hops back to the 
hole site liberating a photon with energy ~ f/ — zja, leaving 
behind a locally disturbed antiferromagnet with z exchange 
bonds broken in the final state as indicated by the dotted 
lines. 
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FIG. 18 The dependence of (a) two-magnon Big Raman in- 
tensity (shown in the inset for uii = 8t) and (b) absorption 
spectrum on the incoming photon energy uji in a 20-site clus- 
ter of the Hubbard model with U = lOt (J = OAt). The solid 
line in (b) is obtained by performing a Lorentzian broadening 
with a width of OAt on the delta functions denoted by vertical 
bars. From Tohyama et al. (2002). 



valence bond picture (Ho et al., 2001).^^ 

These approaches fail when the laser frequency is tuned 
near an optical transition. In this regime, based on 
a spin density wave approach, Chubukov and Frenkel 
(1995a, b) have formulated a so called "triple-resonance" 
theory from which important features of the spectra can 
be derived. Using a SDW approach to the Hubbard 
model, they found that additional resonant diagrams of 
the type shown in Figure 10 (h)-(j) contribute to the 
usual Loudon-Fleury terms, and derived a resonant pro- 
file in good agreement with experiments. In addition, re- 
cent results by Tohyama et al. (2002) have been obtained 
for the Raman response in the resonant limit from both 
spin and charge degrees of freedom. In Figure 18 we show 
their results from exact diagonalization of the Hubbard 
model with a 20-site cluster (Tohyama et al., 2002). The 
two-magnon response at roughly Ml = 2. 7 J is resonantly 
enhanced when the incident photon frequency is tuned to 
the Mott gap scale U, in qualitative agreement with the 
results of Chubukov and Frenkel (1995a, b). Both ap- 
proaches predict a resonant profile for two-magnon Ra- 
man which differs from the absorption profile, as shown 
in Figure 18. Tohyama et al. (2002) have pointed out 
that the resonance energies for the absorption spectrum 
and the two-magnon response are not the same, due to 
differences in SDW coherence factors. 



5. Interacting Electrons - Full Response 

An approach treating the full fermionic degrees of free- 
dom and, simultaneously, treating non-resonant, mixed, 
and resonant scattering on equal footing, is still in its in- 
fancy. The theoretical challenge in calculating the full in- 
elastic light scattering response function is that the mixed 
diagrams involve three-particle susceptibilities and the 
resonant diagrams involve four-particle susceptibilities. 
Only in the infinite-dimensional limit, where most of the 
many-particle vertex renormalizations vanish (all three- 
particle and four-particle vertices do not contribute; only 
the two-particle vertices enter), one can imagine arriving 
at exact results. As an exception, the full Raman re- 
sponse function can be calculated in the Falicov-Kimball 
model, because the two-particle irreducible charge vertex 
is known exactly in the limit of large dimensions (Freer- 
icks and Miller, 2000; Shvaika , 2000). 

Recently, Shvaika et al. (2004, 2005) obtained the fuU 
electronic Raman response function, including contribu- 
tions from the non-resonant, mixed, and resonant pro- 
cesses within a single-band model. In general, the res- 
onance effects can create orders of magnitude enhance- 
ment over the non-resonant response, especially when the 



The reader is also referred to the review article by Lemmens et 
al. (2003) for reviews on the theoretical treatments of magnetic 
light scattering in low-dimensional quantum spin systems. 
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Frequency [t] 

FIG. 19 Raman response in the 1/2-filled Falicov-Kimball 
model {U = 2t) as a function of transfered frequency for vari- 
ous temperatures for fixed incident photon frequency uji — 2t. 
The thickest curve is T/t = 0.05, and the temperature in- 
creases to 0.2, 0.5, and 1 as the curves are made thinner. 
From Shvaika et al. (2004). 

incident photon frequency is slightly larger than the fre- 
quency of the non-resonant feature. The resulting Raman 
response is a complicated function of the correlations, the 
temperature, the incident photon energy, and the trans- 
fered energy. It was found that resonance effects are dif- 
ferent in different scattering geometries, corresponding to 
different symmetries of the charge excitations scattered 
by the light. 

Resonance effects were found as a function of both the 
incoming and the outgoing photon frequencies tUi^s- A 
double resonance - occurring when the energy denomi- 
nators of two pairs of the Greens functions, appearing 
in the bare response shown in Figure 10, approach zero 
- gives the strongest resonant enhancement of the re- 
sponse(Shvaika et ai, 2004). In addition, an interesting 
resonance effect on both the charge-transfer peak and 
the low-energy peak was found when the incident pho- 
ton frequency is of the order of the interaction strength, 
showing that in general the total response cannot be well 
described as a uniform resonance enhancement of the sep- 
arable non-resonant response. In agreement with the re- 
sults of Tohyama et al. (2002) , for an antiferromagnetic 
system this is a direct consequence of the inseparability 
of energy scales in the correlated electron problem, in 



contrast to non-interacting electrons. 

Shown in Figure 19 is the temperature and symmetry- 
dependent Raman response, including non-resonant, res- 
onant, and mixed terms in the Falikov-Kimball model. 
In the insulating phase, spectral weight is depleted for 
small energy transfers and piles up into the excitations 
at energies of order U as the temperature is lowered. The 
transfer of spectral weight from lower to higher energies 
occurs across a temperature independent so-called isos- 
bestic point. An isosbestic point also appears in studies 
of the Hubbard model (Freericks et ai, 2001, 2003a), im- 
plying that it is a generic feature of the insulating phase, 
regardless of the microscopic origin of the phase. We 
note that isosbestic behavior already appears in the non- 
resonant contributions for Big scattering. In the Aig 
and symmetries, it emerges only if resonant terms 
are included. 

The local treatment of self energies in the single-site 
DMFT approach imposes limitations on the theory of 
light scattering in correlated systems. In particular, the 
full polarization dependence of the Raman spectra would 
uncover the way in which correlations affect electron dy- 
namics in regions of the BZ, providing a two-particle 
complement to ARPES, for example. Progress here most 
likely will come from cluster dynamical mean field the- 
ory able to treat nonlocal and anisotropic interactions in 
a coarse-grained manner. 



6. Superconductivity 

As discussed in Section II.D.l, in the absence of inter- 
actions, there is no phase-space for low energy Raman 
scattering for q = momentum transfers. In the super- 
conducting state, phase space restrictions are lifted since 
light can break q = Cooper pairs if the energy of the 
light is greater than 2 A. Thus the Raman response be- 
comes non-trivial, yet easily formulated, in BCS theory. 
As a consequence, there has been an enormous amount 
of theoretical work devoted to light scattering for tem- 
peratures below Tc as an extension of the theory for non- 
interacting electrons in the normal state. We review that 
work here. 

In the superconducting state, focus has been tradi- 
tionally placed on the two-particle non-resonant response 
in BCS theory. Formally the Raman response is given 
by generalizing Eqs. (52-53) in particle-hole space using 
Pauli matrices Ti=o..3 in Nambu notation (Nambu, 1960): 



2 



^ ^ Tr [7(k)G(k, iw)f (k; iw; ^^^)G(k, ioj + in) 



(56) 



25 



where Tr denotes the trace, and 



f(k; iw; iO) = 7(k) + ^ ^ y,(k - k^ - ^wO^'^G(k^ iw^f (k^w^ i^^)G(k^ iw' + iO)fi. (57) 
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Here the bare Raman vertex of coupling to charge is 
^ = and the interaction Vi determines the channel of 
the vertex corrections. For example Vi=z corresponds to 
interactions coupling electronic charge, while Vi=o corre- 
sponds to spin interactions. 

For the case of weak correlations, the Green's functions 
appearing in Eqs. (56) and (57) are given by the BCS 
expression 



G(k, iio) 



iijjTa + e(k)f3 + A(k)-ri 
(ia;)2-£;2(k) ' 



(58) 



with £'2(k) = ^^(k) + A2(k) the quasiparticlc energies. 
In the weak coupling limit for the BCS approximation, 
Vi=3 = — y for phonon- mediated pairing. 



By far, the superconducting state has received the 
largest amount of attention from theory, starting from 
the seminal contribution of Abrikosov and Fal'kovskii 
(1961), which predated the observation of the effect by 
19 years. The main focus in the early years was to study 
the 2A features in conventional s— wave superconduc- 
tors with small and large coherence lengths (Abrikosov 
and Fal'kovskii, 1961, 1987: Abrikosov and Falkovsky, 
1988; Klein and Dicrkcr, 1984), including the effects of 
Coulomb screening (Abrikosov and Genkin, 1973) and 
examining the temperature dependence (Tilley, 1972). 
If one neglects vertex corrections, the q-dependent Ra- 
man response in a superconductor is given by a projected 
Maki-Tsuneto function (Maki and Tsuneto, 1962), 
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with the coherence factors A± (k, q) = 1 ± ^('^)?(^^q)^^(^)^'^('^+q) _ more common expression is the Raman response 
for q = which simplifies to 



Xa,6(q = 0, + i^) = ^ ak6k 
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The hill Raman response, including charge screening, is 
once again given by Eq. (29), in which the vertices a,b 
are replaced by the Raman (a, 6 = 7) and pure charge 
(a, 6 = 1) vertices (Abrikosov and Genkin, 1973). 

For the case of an isotropic gap (A(k) = A) and mo- 
mentum transfers q = 0, a threshold and a square-root 
discontinuity appears at twice the gap edge A, reflecting 
the two-particle density of states. For finite q, the singu- 
larity is cut-off due to breaking Cooper pairs with finite 
momentum, and the peak is shifted out to frequencies of 
roughly v-pq as in the normal state (Figure 20). Qualita- 
tively similar behavior is obtained for disordered s— wave 
superconductors (Devereaux, 1992) in which 1/t£ (see 
Eq. (49)) assumes the role oiv^q ■ 

Further advances in the theory for conventional 
s— wave superconductors were made for energy gaps with 



small anisotropy (Klein and Dierker, 1984), coexistence 
with charge-density wave order (Balseiro and Falicov, 
1980; Littlcwood and Varma, 1981, 1982; Tiitto and Za- 
wadowski, 1992), layered superconductors (Abrikosov, 
1991), impurities (Devereaux, 1992, 1993), and final-state 
interactions (Devereaux, 1993; Devereaux and Einzel, 
1995a; Klein and Dierker, 1984; Monien and Zawadowski, 
1990). 

We note in particular that the variation with k of the 
Raman vertices 7(k) is coupled to the k-dependence of 
the energy gap A(k) (see, e.g., Eq. (60)), leading to 
a strong polarization dependence of the spectra. For 
isotropic s-wave superconductors, the vertex does not af- 
fect the lineshape, and thus the spectrum is polarization 
independent, apart from an overall prefactor. For this 
case, a polarization dependence can be generated in BCS 
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theory by taking into account channel-dependent final- 
state interactions (Bardasis and SchriefFer, 1961) and/or 
impurity scattering. However, for the most part this only 
produces a channel dependence in the vicinity of the gap 
edge, and thus the main feature of the response is the 
uniform gap existing for all polarizations. For anisotropic 
energy gaps, the symmetry dependence of the spectra is 
a direct consequence of the k-summation (angular aver- 
aging) , which couples gap and Raman vertex and leads to 
constructive (destructive) interference if the vertex and 
the gap have the same (different) symmetry. 

Generally, in superconductors with nodes of the energy 
gap, power-laws in the low frequency and/or tempera- 
ture variation of transport and thermodynamic quanti- 
ties emerge, replacing threshold or Arrhenius behavior 
ubiquitous in isotropic superconductors. However, due 
to the averaging over the entire Fermi surface, the power- 
laws themselves do not uniquely identify the ground state 
symmetry of the order parameter, but only can give the 
topology of the gap nodes along the Fermi surface, e.g., 
whether the gap vanishes on points and/or lines. Thus, 
one cannot distinguish between different representations 
of the energy gap which have the same topology. For in- 
stance, for the case of d-w&ve tetragonal superconductors, 
there are five pure representations which have line nodes 
on the Fermi surface. Two-particle correlation functions, 
determining the density, spin or current responses, do not 
have the freedom to probe various portions of the gap or 
its phase around the Fermi surface. 

With the advent of high- Tc cuprates, a flurry of ac- 
tivity ensued on theory of Raman scattering in d— wave 
superconductors (Falkovsky, 1990; Monien and Zawad- 
owski, 1989), specifically including polarization depen- 
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dences (Devereaux et al., 1994a; Devereaux and Einzel, 
1995a), collective modes (Dahm et ai, 1998; Devereaux 
and Einzel, 1995a; Wu and Griffin, 1995a, b), impurities 
(Devereaux, 1995b; Devereaux and Kampf, 1997; Dev- 
ereaux, 2003b; Wu and Garbotte, 1998), temperature 
dependences (Branch and Garbotte, 1995; Devereaux, 
1995b; Devereaux and Kampf, 1997; Devereaux, 2003b), 
screening (Branch and Garbotte, 1996; Devereaux and 
Einzel, 1995a; Manske et ai, 1997; Strohm et ai, 1998a), 
band structure and bi-layer effects (Branch and Garbotte, 
1996; Devereaux et al., 1996; Krantz and Gardona, 1994; 
Strohm and Gardona, 1997), surface and c-axis contri- 
butions (Wu and Griffin, 1996; Wu and Garbotte, 1997), 
resonant effects (Sherman et ai, 2002), and mixed-state 
pairing (Devereaux and Einzel, 1995a; Lee and Ghoi, 
2002; Nemetschek et al., 1998). The focus was largely 
on how the symmetry selection rules could locate the po- 
sitions of the gap maxima and nodal points around the 
Fermi surface. 
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FIG. 20 Raman rosporiso for an s— wave superconductor for 
= 0.1,0.5,1.0,2.0,4.0, and 8.0 with ^ = Hvp/tvA the 
Pippard-BCS coherence length ^ and g ~ 1/5 the momentum 
transfer in a metal with skin depth S. From Klein and Dierker 
(1984). 



FIG. 21 Raman spectra for miscrccnod Aig (top), Big (mid- 
dle) and B2g (bottom) response as a function of reduced tem- 
perature t = T/Tc- Here 2 A ~ 37meV, and the higher peak 
in Aig and Big channels is a van Hove feature. From Branch 
and Garbotte (1995). 
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For (ij.2_j,2 superconductors, the interplay of polariza- 
tions and gap anisotropy can be simply drawn. Refer- 
ring to Figure 7, Big orientations project out excitations 
around the principle directions (M-points or anti-nodal 
regions) of the BZ where the superconducting gap is max- 
imal and where the van Hove singularity is located, while 
B2g orientations project the nodal regions along the di- 
agonals. As a consequence, the Raman response has a 
peak at 2Amax for Big and at slightly lower energy for 
B2g. The polarization dependence also enters the low fre- 
quency behavior. Since line nodes yield a linear depen- 
dence on energy of the density of states, the B2g response 
depends linearly on hfl in the limit f7 ^ for a gap van- 
ishing on the diagonals. For Big orientations, in contrast, 
the Raman vertex vanishes along with the energy gap at 
the same points in the BZ. This yields an additional f2^ 
contribution from the line nodes of the vertex, and the 
resulting response varies as fi^. The unscreened Aig re- 
sponse measures an overall average throughout the BZ 
and thus picks up the gap maxima as well as the linear 
density of states. This is shown quantitatively in Fig- 
ure 21. The frequency power laws also translate in low 
temperature power laws of the response in the small fre- 
quency limit (Devereaux and Einzel, 1995a). 

Disorder effects generally smear peak features at larger 
energies and change the Big exponent to 1, similar to 
the change in the low temperature NMR rate for d-wave 
superconductors (Devereaux, 1995b). Moreover, simi- 
larities between the in-plane conductivity and B2g Ra- 
man follows from the BZ weighting around the nodes, 
while the Big response is qualitatively similar to the c- 
axis conductivity due to the weighting around the anti- 
nodes (Devereaux, 2003b). For example the residual in- 
plane conductivity as T ^ is universal and given by 
a{T = 0) = ne"^ /mnAo, the slope of the B2g response 
2Np/nAo is also universal and insensitive to impurity 
effects, while the Big channel and c-axis conductivity 
are non-universal, having additional impurity dependent 
prefactors (Devereaux, 1995b; Devereaux and Kampf, 
1997; Devereaux, 2003b). The slope of the B2g response 
follows the temperature dependence of the in-plane con- 
ductivity, and both possess a peak at intermediate tem- 
peratures due to a balance of DOS and lifetime effects 
as temperatures are lowered from Tc. Yet, both the out- 
of-planc conductivity and the Big response do not show 
a peak due to the more rapid variation of the projected 
DOS coming from antinodal portions of the BZ. 

We note that for a dxy energy gap, the above discussion 
applies accordingly, with the role of Big and B2g sym- 
metry reversed. It was indeed an important development 
to show that Raman scattering is unique in determining 
two-particle electron dynamics independently in different 
regions of the BZ in the superconducting state. 

While the low frequency power-laws are insensitive to 
band structure (such as the shape of the Fermi surface), 
the polarization selection rules can select features of the 
band structure at higher energies. For example, in the 
cuprates, the van Hove singularities from (tt, 0) and re- 



lated points yields a peak at twice the quasiparticle en- 
ergy E(k) in Aig and Big channels, as shown in Figure 
(21). For multiple bands, including the case of several 
Fermi surface sheets, the responses for crossed polariza- 
tions are simply additive, yet for Aig channels due to 
backflow effects an additional interference term can be 
present if the charge fluctuations are different for the dif- 
ferent sheets (Devereaux et al, 1996; Krantz and Car- 
dona, 1994). 

Quite generally the backflow yields substantial reor- 
ganization of spectral weight around 2Aniax compared 
to the bare response (Branch and Carbotte, 1995; Dahm 
et al., 1999). This is because the term X7,i) which con- 
tributes in channels having the symmetry of the lattice, is 
peaked and large at the same position as the unscreened 
term x-yn- ^'^^ surprisingly, the reorganization depends 
delicately on the relative momentum dependences of the 
Raman vertices and energy gap (number of BZ harmonics 
for example, as shown in Figure 22), as well as on details 
of the band structure (Branch and Carbotte, 1996; Dev- 
ereaux et al., 1994a; Krantz and Cardona, 1994; Strohm 
and Cardona, 1997). 

For the cuprates, Raman vertices have been calculated 
using LDA (Strohm and Cardona, 1997), but limited 
progress has been made in including the contributions of 
substantial electronic correlations. In most other cases, 
either the effective mass approximation has been used in 
calculations or simply a symmetry classification has been 
made. While a detailed lineshape analysis can be applied 
based purely on symmetry as explained above, it must be 
kept in mind that even a comparison of overall intensities 
between different geometries can at best be qualitative. 

Yet, the continuation of these treatments from the su- 
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FIG. 22 Comparison of the screened Aig response ob- 
tained for different number of d^2_y2 gap harmonics A(k) = 
Ao{[cos(A;2,a) — cos(fcj,a)]/2 + Ai[cos(fcxa) — cos(fc3,a)]^/8-|- 
A2[cos(fc2,a) — cos(fcya)]'''/32}. Here a — d correspond to the 
set Ai,2 = (0, 0), (1, 0), (0, 1), and (1, 1), respectively, and Ao 
has been rescaled to give the same value for the maximum 
gap. Generally the large peak at 2A (as shown in Figure 21) 
is suppressed by the backflow terms. From Devereaux et al. 
(1996). 
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perconducting to the normal state is not straightforward. 
As can be seen directly from Eq. (60), the intensity van- 
ishes proportional to as T approaches Tc- Hence, to 
avoid phase space limitations in the absence of Cooper 
pairs, an additional source for electronic scattering, such 
as the one mediating the formation of Cooper pairs, must 
be included. While strong coupling extensions of Raman 
scattering in rf— wave superconductors have recently been 
presented (Dahm et al, 1999; Devcrcaux and Kampf, 
2000; Jiang and Carbotte, 1996), a merging of the nor- 
mal and the superconducting states is poorly understood. 
This would require a not yet existing microscopic descrip- 
tion of the formation of d— wave superconductivity from 
the normal state. 



7. Collective Modes 

Raman scattering has the almost unique ability to sort 
out collective modes of the two-particle response in dif- 
ferent symmetry channels, owing to the freedom to inde- 
pendently adjust the two polarization vectors. The col- 
lective mode spectrum one obtains depends upon which 
interactions are included in Eq. (57). We first discuss 
the general consequences based on gauge-invariance and 
focus on exciton-like modes. 

In order to form a fully gauge-invariant theory, the in- 
teractions responsible for superconductivity appear not 
only in G, but must also be included as vertex renormal- 
izations F. In this way, the Raman response from pure 
charge density fluctuations in the superconducting state 
yields the Goldstonc mode from the broken gauge sym- 
metry - the phase or Anderson-Bogoliubov mode (An- 
derson, 1958; Bogohubov et al., 1959; Nambu, 1960). In 
the absence of the long-range Coulomb interaction, this 
mode is a soft sound mode, yet the Coulomb interac- 
tions - inescapable for q=0 - push the sound mode up 
to the plasma frequency via the Higgs mechanism. As 
a result, particle- number conservation is satisfied in the 
superconducting state and Xsc(q = 0; f^) — 0; indepen- 
dent of whether one considers Bloch states (Abrikosov 
and Genkin, 1973; Klein and Dierker, 1984; Monien and 
Zawadowski, 1990) or Anderson exact eigenstates of the 
disordered problem (Devereaux, 1993). 

However, additional modes of excitonic origin may ap- 
pear if one considers further interactions between elec- 
trons in clean (Bardasis and Schrieffer, 1961) and dis- 
ordered (Fulde and Strassler, 1965; Maki and Tsuneto, 
1962) conventional superconductors. These excitons ap- 
pear split off from the continuum at Ml < 2A if the 
interaction occurs in higher momentum channels orthog- 
onal to the BCS condensate. 

Since Raman scattering couples to anisotropic charge 
density fluctuations with symmetry selectivity to differ- 
ent channels L, the light polarizations can be used to 
determine the exact nature of bound states. Balseiro 
and Falicov (1980) considered the formation of a phonon- 
Cooper pair bound state due to electron-phonon coupling 



though neglecting channels higher than L = 0. How- 
ever, this mode is canceled by the backflow applying 
generically to all systems. Finite L exciton formation 
in clean and disordered superconductors, and the result- 
ing appearance in Raman scattering, have been consid- 
ered explicitly by Monien and Zawadowski (1990) and 
Devereaux (1993), respectively, bringing the symmetry 
of the exciton and the polarization dependence to light. 

We show in Section III.B below that the effect of fi- 
nal state interactions can be substantial in strongly cou- 
pled conventional superconductors. This demonstrates 
the strength of the electron-phonon coupling, not only 
in general, but also specifically in channels orthogonal 
to the ground state. Interestingly, the lattice instability 
found in some of these materials has the same symme- 
try as the collective mode and the electronic states which 
apparently drive the transition (Weber, 1984). 

For d.j.2_y2 superconductors, the collective mode spec- 
tra have been investigated thoroughly by Devereaux and 
Einzcl (1995a) and others (Dahm et al, 1998; Manske 
et al, 1997, 1998; Strohm et al, 1998a; Wu and Griffin, 
1995a,b). It was shown that the Anderson-Bogoliubov 
mode appears in Aig channels and massive modes can 
appear in other channels. Since the pair state has 
only one representation in the Z)4h group, massive col- 
lective modes arise when one considers interactions in 
orthogonal channels. Recently, it has been suggested 
that the presence of collective modes may allow one 
to distinguish charge- or spin-mediated c?— wave pairing 
(Chubukov et al., 1999, 2006), highlighting the possible 
importance in the context of the cuprates. 

Generally, the collective mode spectrum can be quite 
diverse in unconventional superconductors. In principle, 
additional broken continuous symmetries can exist, such 
as 5O3 spin rotational symmetry in spin-triplet systems 
and orbital rotational symmetry in spin-singlet sys- 
tems, if the gap does not possess the full symmetry of the 
lattice. Furthermore, massive collective modes can arise 
if the energy gap is degenerate or has an admixture of dif- 
ferent representations of the point group. The massive 
modes can in principle lie below the gap edge, and can 
thus be relevant for the low frequency dynamics of cor- 
relation functions. In fact, Raman-active modes in spin- 
triplet superconductors such as Sr2Ru04 have drawn re- 
cent theoretical interest (Kee et al., 2003), although the 
experimental challenges are not negligible because of the 
low Tc and the related small energy gap in these materi- 
als. 

Though very interesting, spin-triplet pairing or spin- 
orbit effects are rare and more on the exotic side in 
superconductivity. Competing ground states, however, 
are quite common whenever correlation effects come into 
play. This is not at all confined to the cuprates, but 
occurs also in, e.g., spin (SDW) and charge (CDW) den- 
sity wave systems. Usually, density- wave formation with 
long-range order at least partially suppresses supercon- 
ductivity such as in 2H — NbSe2 (sec below). Then, ad- 
ditional modes appear as a result of the competition be- 
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tween CDW ordering and superconductivity, and collec- 
tive modes appear as one modulates either one or both 
order parameters. 

Littlewood and Varma (1981, 1982) and Browne and 
Levin (1983) considered a direct coupling between charge 
density and superconducting gap amplitudes, modulated 
for example by a CDW phonon, although this was not 
specified. They obtained an additional "gap" mode be- 
low 2A. Yet this mode was only considered in the L = 
channel, and Coulomb interactions once again remove 
this mode. Lei et al. (1985) considered an effective CDW- 
SC coupling via a phonon and realized that in anisotropic 
systems such collective modes in L ^ channels may 
appear. Finally, Tiitto and Zawadowski (1992) treated 
electron-phonon and CDW amplitude-phonon coupling 
on equal footing in finite angular momentum channels, 
showing generally that the collective modes in these chan- 
nels are unaffected by Coulomb screening. The modes 
obtained split off from the gap edge and appear as exci- 
tations below the quasiparticle spectrum, much like ex- 
citons in semiconductors. Evidence of mixed CDW-SC 
pairing may be seen in Raman experiments via the pres- 
ence or absence of these modes. This has recently been 
extended by Zeyher (2003) to MgB2, having multiple en- 
ergy gaps on different electron bands. 

The collective mode spectrum of coupled d— wave 
charge density and superconductivity was investigated 
by Zeyher and Greco (2002) along the lines developed by 
Tiitto and Zawadowski (1992) for conventional CDW and 
superconducting systems. As for s— wave CDW super- 
conductors, collective modes split off from the maximum 
of the gap edge. As an important difference in d— wave 
systems, the modes distinctly affect the various symme- 
try channels. Besides the reorganization of Aig spec- 
tral weight, additional modes alter the Big spectrum, as 
shown in Figure 23. 

Density waves instabilities need not necessarily com- 
pete with superconductivity but rather can provide an ef- 
fective coupling mechanism (Castellani et al., 1995; Perali 
et al., 1996) as long as quantum and thermal fiuctuations 
suppress long range order. Then, collective modes may 
appear as fluctuation-induced modes. The Raman scat- 
tering is usually determined from Aslamazov-Larkin fluc- 
tuation diagrams considered for the conductivity (Asla- 
mazov and Larkin, 1968). To overcome the q = phase 
space limitations, the Raman response is given by the 
exchange of two fluctuations modes at wavevectors Qc 
and — Qc, yielding generally a mode at energies of twice 
the mass of the fluctuation propagator. Once again the 
polarization dependence can select different fluctuation 
modes corresponding to different ordering wavevectors 
coupling to either charge or spin density modes. This 
was investigated for spin (Brenig and Monien, 1992) and 
charge (Caprara et al., 2005) fluctuations in the normal 
state, and for novel spin resonances in the superconduct- 
ing state of the cuprates (Chubukov et al., 1999, 2006; 
Venturini et al, 2000). 

Lastly, we remark that many other types of collective 



modes arc possible if one considers more exotic ground 
states with different symmetry classifications. For ex- 
ample, a chiral spin liquid has been investigated by 
Khveshchenko and Wiegmann (1994) in which helical ex- 
citations were conjectured to exist and are in principle 
measurable in A2g orientations which can be projected 
out via proper sums of spectra taken with both linearly 
and circularly polarized light. Other examples are modes 
induced by magnetic fields or optical modes resulting 
from Dzyaloshinskii-Moriya interactions in Heisenberg 
antiferromagnets as observed recently in lightly doped 
La2-xSrxCu04, < .x < 0.03 (Gozar et al, 2005b) and 
discussed by Silva Neto and Benfatto (2005), directly 
demonstrating the importance of spin coupling to the 
local environment. 



III. FROM WEAKLY TO STRONGLY INTERACTING 
ELECTRONS 

In this section we review experimental results in sys- 
tems other than doped semiconductors (see reviews by, 
e.g., Abstreiter et al. (1984) and Pinczuk and Abstreiter 
(1989)) and cuprates (see section IV) with a view to- 
wards signatures in the Raman spectra arising from the 
development of strong electronic correlations. We discuss 
various types of superconductors and summarize results 
on correlated metals and other strongly interacting sys- 
tems. 

The light scattering cross section in absorbing media, 
such as systems with free carriers, is generally weak since 
the interaction volume is small for the short penetration 
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FIG. 23 Electronic Raman spectra of B symmetries for three 
different doping levels including coupling of d-CDW and d- 
superconducting amplitudes. The superscript denotes the 
spectra in the absence of collective modes. From Zeyher and 
Greco (2002). 
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depth of visible light, 5 Xi = 2Trc/uJi. As a conse- 
quence, the momentum perpendicular to the surface is 
not conserved, and the transfer q is not given any more 
by the difference of the vacuum momenta of the involved 
photons kj — ks but essentially hy S = A/(47rfc) with k the 
imaginary part of the index of refraction (Abrikosov and 
Fal'kovskii, 1961; Mills et ai, 1970). Even in strongly ab- 
sorbing materials with k > 1, l/S <^ ir/a holds where a 
is the lattice constant, and the limit of small momentum 
transfer is still effective. This introduces a new energy 
scale hvpq hvp/S with vp and q being the magnitudes 
of the Fermi velocity and the momentum transfer, respec- 
tively. In all considerations, this scale must be put into 
relation to the other relevant energies, such as the elec- 
tron scattering rate T = h/r in the normal and the gap A 
in the superconducting state. These apparent academic 
considerations have major impact on both the observ- 
ability and the interpretation of electronic spectra. 



A. Elemental Metals and Semiconductors 

In addition to the small scattering volume due to the 
absorption of light by free carriers, a parabolic dispersion 
and a spherical Fermi surface reduce the cross section 
of single-electron excitations in metals and degenerate 
semiconductors strongly, since in such systems the asso- 
ciated density fluctuations are screened by the long-range 
Coulomb interaction. The few spectra wc are aware of 
have been taken on elements with a more complex band 
structure such as Nb (Klein, 1982a; Klein and Dierker, 
1984) or Dy (Klein et al, 1991). In Dy a broad con- 
tinuum similar to that in the high- Tc cuprates (Bozovic 
et al., 1987) is found. In Nb the superconducting state 
was studied. Due to the low transition temperature Tc, 
the correspondingly small energy gap A(T) and the small 
ratio with 5 the penetration depth of the light and 
^ the superconducting coherence length (for details see 
section III.B) the characteristic redistribution of scatter- 
ing intensity is very hard to observe. The peaks found 
at 1.8 K in the expected energy range close to 2A(r) are 
very weak, and no normal state spectra have been mea- 
sured for comparison (Klein, 1982a; Klein and Dierker, 
1984). 

In fact, superconductors rather than normal metals 
were the main focus in the early days of electronic Ra- 
man scattering. Only after the discovery of the cuprates 
(Bednorz and Miiller, 1986), with generally complicated 
and sometimes very surprising electronic properties, did 
studies of the normal state become increasingly attrac- 
tive (see section IV). 



B. Conventional Superconducting Compounds 

Among superconductors, intermetallic compounds like 

NbsSn or VsSi with A15 structure, can be considered 
conventional both above and below Tg. They are strictly 



3D, superconductivity is mediated by phonons leading to 
an essentially isotropic s-wave gap, and correlations are 
believed to be of minor importance. This does not mean 
they are simple. For instance, the Fermi velocity is very 
small, close to the velocity of sound, and the Fermi sur- 
face is multi-sheeted. Sufficiently perfect single crystals 
of NbsSn and VsSi undergo a structural transformation 
from a cubic to a tetragonal lattice at low temperature. 
Nevertheless, A15 compounds are paradigms of strong- 
coupling s-wave superconductors with a high density of 
electronic states at Ep. Materials like the borocarbides, 
MgB2 or 2H — NbSe2 are certainly more complex and 
correlations or multiband aspects come into play. 

1. A15 compounds 

Superconductivity induced structures close to twice 
the gap edge were found in mono-crystalline NbsSn 
(Fig. 24) and V.-jSi (Dierker et al, 1983; Hackl et al, 
1982, 1983; Klein, 1982a) two years after the discovery 
of gap modes in 2ff — NbSe2 by Sooryakumar and Klein 
(1980) (section III.C) and after an early but imsuccess- 
ful attempt in polycrystaUine NbaSn by Fraas (1970). 
For T < Tc the scattering intensity is redistributed 
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FIG. 24 Raman spectra of NbaSn. Lower curves in (a) and 
(b) arc at 40 K and upper curves are at 1.8 K. Data in (c) 
and (d) are at 1.8 K. Below Tc = 18 K the intensity at low 
energies is strongly suppressed with respect to the normal 
state. Beyond a threshold of approximately 50 cm~^ a new 
peak appears. The symmetries are Eg (a, d), T2g (b), Eg-\-Aig 
(c top). Eg (c middle), and Aig (c bottom). The A\g data 
in (c) are obtained by subtracting the middle from the upper 
curve. All solid lines except for the upper two in (c) are 
theoretical fits to a broadened Maki-Tsuneto function (see 
section II.D.6). Prom Dierker et al (1983). 
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with a suppression below and a pile-up at approximately 
2 A ~ 50 cni^^. The well-defined peak in Eg symmetry 
follows the BCS prediction for the temperature depen- 
dence of the gap up to approximately 0.85 Tc (Hackl 
et al, 1983, 1989). Somewhat unexpectedly, the peak 
frequencies of the superconductivity-induced features de- 
pend on the selected symmetry (Table II). Independent 
of minor differences in the absolute numbers stemming 
from the data analysis the Eg peaks are significantly 
below those having Aig and symmetries. At first 
glance one could think of a gap anisotropy to manifest 
itself. However, there is no support from the tunneling 
results which rather indicate the possible gap anisotropy 
to be opposite in VaSi and NbsSn and very large or from 
calorimetric studies which should track the smallest gap 
(Table II). In addition, the shapes of the Raman spectra 
are strongly symmetry dependent in that the Eg peak 
is much narrower than the others. The meaning of this 
anisotropy was a matter of intense discussion. 

The results in Aig scattering symmetry in NbsSn 
(Dierker et al, 1983; Klein, 1982a) and later in VsSi 
(Hackl and Kaiser, 1988a) demonstrate clearly that the 
structures below Tc originate in light scattering from 
Cooper pairs (Fig. 24), since there exist no Raman ac- 
tive excitations at this symmetry, such as phonons or 
other bosonic modes in the A15 structure from which 
the electrons can borrow intensity. There is not even an 
electronic continuum above Tc (see Fig. 24 (c) and Hackl 
and Kaiser (1988a)). In spite of similar band structures 
and densities of states at the Fermi level, N{Ef), (Klein 
et al., 1978) the intensities of the modes are quite dif- 
ferent in the two compounds as is the overall scattering 
cross section. For this reason, the weak Aig mode in VsSi 
escaped detection for a while (Hackl and Kaiser, 1988a). 
Since there is nothing to interact with the peak frequen- 
cies, the Aig structures should be close to the energy gap 
in the respective material. One actually observes coin- 
cidence of both the Aig and T2g Raman energies with 
those of bulk methods such as calorimetry and neutrons, 
while the Eg energies are substantially lower (Table II). 

We first note that surface sensitive methods such as 
tunneling return somewhat smaller gap energies than 
bulk methods. Optical spectroscopy results are also 
smaller most likely due to surface treatment. Strain 
or disorder can indeed reduce Tc in A15 materials since 
N{Ef) decreases rapidly (Mattheiss and Weber, 1982). 
Similar reasons might apply for the Raman data in NbsSn 
of Dierker et al. (1983) although the fits (see Figure 24) 
reveal gap values slightly below (5-10 %) the peak posi- 
tions. Spectra of cleaved surfaces, such as those of VsSi 
and NbsSn taken by Hackl and Kaiser (1988a) and Hackl 
et al. (1989), respectively, apparently give gaps closer to 
the bulk values. 

For these reasons, it seems worthwhile to look for other 
sources of the anisotropy, and we consider an interpreta- 
tion in terms of final state interactions (Bardasis and 
Schrieffer, 1961; Klein and Dierker, 1984; Zawadowski et 
al., 1972). This means that the two single electrons of 



Sample 


Raman Energy (cm ^) 

Aig Eg T2g 


Reference Data (cm ^) 


NbsSn 


52 41 50 a 
67 48 70 b 


50 tunneling e 

35-13 tunneling f 
53 tunneling g 
62 calorimetric h 
56 neutrons i 


VsSi 


40 - c 
55 42 52 d 


37 tunneling j 
40-50 tunneling k 
46 IR 1 
41 IR m 
49 calorimetric h 



TABLE II Gap energies in A15 compounds as measured by 
Raman scattering and other methods, a and c refer to re- 
sults from fits (sec Fig. 24 and section 11. D. 6), b and d arc 
peak frequencies. In two cases an anisotropy was found by 
tunneling being indicated by a range (f and k). The first 
and the second numbers are for [100] and [111] directions, 
respectively. Results of the following publications are used: 
a (Dierker et al, 1983), b (Hackl et al, 1989), c (Klein and 
Dierker, 1984), d (Hackl and Kaiser, 1988a), e (Rudman and 
Beasley, 1984), f (Hoffstein and Cohen, 1969), g (Geerk et 
al, 1984), h (Junod et al, 1983), i (Axe and Shirane, 1973), 
j (Moore et al., 1979), k (Morita et al., 1984), 1 (Tanner and 
Sievers, 1973), m (Perkovitz et al., 1976). 



a broken Cooper pair can still interact in channels or- 
thogonal to the pairing channel. The strongly coupled 
Eg phonon (Schicktanz et al., 1980, 1982; Weber, 1984; 
Wipf et al, 1978) is in fact orthogonal to the fully sym- 
metric (s-wave) pairing channel. Hence, it is capable of 
forming a bound state below the pair-breaking threshold, 
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FIG. 25 Raman spectra in Eg symmetry of VsSi. From 
Monien and Zawadowski (1990). 
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explaining both the reduced energy and the hnewidth of 
the Eg gap mode (Monien and Zawadowski, 1990). Fits 
to the results in VsSi are substantially improved by in- 
cluding the bound state (Fig. 25) in comparison to those 
neglecting it (Klein and Dierker, 1984). Additional ex- 
perimental support comes from the evolution with tem- 
perature of the spectra in VsSi and NbaSn (Hackl et at, 
1983, 1989). In either compound, the integrated spectral 
weight in Aig symmetry increases significantly because a 
new scattering channel opens up below due to the for- 
mation of Cooper pairs while staying essentially constant 
in Eg symmetry because the weight is being transferred 
from the phonon to the bound state (Fig. 26). 

In contrast to Eg symmetry, the pair-breaking features 
in symmetry are weak and essentially at the Aig po- 
sition. The question arises as to why there is no bound 
state although there exists a phonon. Clearly, the T2g 
phonon intensity is weak and the line width is small, 
reflecting the moderate coupling as opposed to Eg sym- 
metry where the complete line width and the asymmetric 
Fano shape stem from the coupling to conduction elec- 
trons (Weber, 1984; Wipf et ai, 1978). The bound state's 
energy split-off by approximately 30 % indicates that the 
very strong interaction drives the system close to an in- 
stability of the s-wave ground state. On the other hand, 
the T2g mode is only weakly coupled and the interaction 
with the conduction electrons is not strong enough to 
substantially renormalize the spectrum. 

Symmetry arguments, the unique line shape and the 
intensity transfer in Eg symmetry, as well as the com- 
parison to calorimetric results, make us believe that the 
formation of a bound state is more likely an interpreta- 
tion of the Eg results in A15 compounds than the man- 
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ifestation of a gap anisotropy or a two-gap scenario. In 
contrast, both effects may cooperate in 2D MgB2 discov- 
ered to be a superconductor just recently by Nagamatsu 
et al. (2001). 



2. MgB2 and the Borocarbides 

Electronic Raman studies on MgB2 have explored the 
superconducting energy gap and changes in phonon line- 
shapes occurring below Tc, starting with the work of 
Chen et al. (2001) and followed thereafter by Quilty et al. 
(2002, 2003). Here the symmetry dependence of the re- 
sponse allowed a direct observation of the pairing gap 
on the two-dimensional a bands and the 3D tt bands. By 
orienting the light polarizations along the c-axis of MgB2 
(perpendicular to the hexagonal planes) the a bands can- 
not be probed for having little dispersion and thus the tt 
bands are projected out, giving a value 2A7r = 29 cm~^. 
Other polarizations are able to detect a larger pairing 
gap 2A(j = 100 cm~^. zz-polarized spectra in the super- 
conducting state are shown in Figure 27 along with fits 
from the theory for disordered s— wave superconductors. 

The gap values are consistent with those from other 
techniques, yet the fit yields values of the disorder-related 
scattering rate different from those of the resistivity by 
a factor of 2 (Quilty et al, 2003). Zeyher (2003) has re- 
analyzed the fit where the direct coupling of light to the 
a band is zero and the cr gap appears as a result of a 
coupling to the Raman active E2g phonon, believed to 
be largely responsible for pairing. The xx spectrum in 
the superconducting state can be understood then as a 
superposition of a phonon line, a background, and a col- 
lective bound state due to residual interactions between 
electrons, similar to that observed in A15 compounds (see 
Figures 25 and 26 (a)). 

The superconducting energy gap has also been stud- 
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FIG. 27 Raman spectra for xx and zz polarization geometries 
in the superconducting state of MgB2 taken by Quilty et al. 
(2003). The solid lines are fits to the data using the theory 
of Devereaux (1992) for disordered s— wave superconductors 
using a single gap and two gap model for zz and xx polariza- 
tions. The XX spectrum has also been interpreted in terms of 
a collective bound state by Zeyher (2003). 
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ied in some detail in the borocarbidc superconductors 
i?Ni2B2C iR=Y, Lu) by Yang et al. (2000a,b). Sharp 2A 
peaks were observed in Aig and B2g symmetries, while 
the maximum in Big symmetry is less pronounced and 
20% higher in energy. All peaks showed a typical BCS- 
type temperature dependence and disappeared above the 
upper critical field iJc2- Due to the high surface quality 
and improved instrumentation, the residual scattering in- 
tensity below the gap edge is much smaller than, e.g., in 
the A15 compounds but finite with an approximately lin- 
ear variation with energy. 

Since a direct coupling to a Raman active mode was 
not found in the borocar bides, it is more complicated 
than in the A15 compounds or in MgB2 to sort out 
whether a gap anisotropy, multi-gap superconductivity 
or collective modes are responsible for the variations in 
energy and line shape at the different symmetries. Of 
course, the existence of bound states is not related to 
Raman-active modes; only the experimental verification 
is more indirect, e.g., via the line shape (see Figure 25). 
Similarly, the linear low-energy scattering can suggest 
either the presence of strong inelastic scattering due to 
large coupling constants A (Allen and Rainer, 1991; Yang 
et al., 2000b) or gap nodes in pairing states with lower 
symmetry, such as s -I- g-wave superconductivity (Lee and 
Choi, 2002). A quantitative analysis on the basis of a re- 
alistic band structure could possibly help clarify these, 
at present, open issues. 

While superconductivity was the dominant correlation 

in the A15 compoimds, MgB2 was more complex due to 
the interplay between 2D and 3D behavior. In the boro- 
carbides magnetic order as a second instability compet- 
ing with superconductivity comes into play (Canfield et 
al., 1998). Although the Raman studies were performed 
on non-magnetic compounds (Yang et al, 2000a, b) the 
vicinity of different types of order is characteristic for this 
and the following classes of systems. 



C. Charge Density Wave Systems 

The competition or coexistence of different ground 

states was studied intensively in layered dichalcogenides 
in the 1970s and 1980s. The interest in these charge 
density wave (CDW) systems was revived after the dis- 
covery of superconductivity in the cuprates for two rea- 
sons: in both compound classes superconductivity com- 
petes with one or more other instabilities, and, secondly, 
the dramatically improved instrumentation allowed qual- 
itatively new and unexpected insights into materials like 
2J?— NbSe2. As an example, the electronic scattering 
rate exhibits marginal (Varma et al, 1989a) rather 
than Fermi liquid like temperature and energy depen- 
dences (for a discussion and for references see, e.g., Cas- 
tro Neto (2001)). 



1. 2i?-NbSe2 

2/?— NbSe2 is a layered, though 3D, superconductor 

with an in-plane coherence length ^|| of approximately 
70 A and = 25 A (de Trey et al, 1973). The pene- 
tration depth for visible light 5 is of the order of 200 A, 
hence ^ 5. The discontinuity at 2A is expected to in- 
crease with 2A/ {hvpq) w S/£, (Klein and Dierker, 1984). 
In addition, the material can be cleaved easily, facilitat- 
ing the preparation of atomically flat surfaces from which 
diffuse scattering of laser light is minimized. These are 
favorable (though not easy!) conditions for observing gap 
structures close to the elastic line. 

In Fig. 28, the first observation of the redistribution 
of scattering intensity in the superconducting state of 
2_ff — NbSe2 with the sample immersed in superfluid He is 
reproduced. The effect is measured for two samples with 
slightly different impurity concentrations. In either case, 
the fully symmetric A and the E responses are shown. 
The peaks have slightly different energies, and are located 
at 18 and 15 cm"-*^, respectively, close to the essentially 
k independent leading edge gaps found in recent photoe- 
mission experiments (Valla et al., 2004). In the normal 
state at 9 K the new low energy modes are absent, while 
the maximum related to the CDW state is still present. 




FIG. 28 Raman spectra of 2//— NbSc2 above and below Tc- 
The lower curves of each panel are taken in the normal state 
(9 K) the upper ones at 2 K well below Tc with the sample 
immersed in superfluid He. At the A (|| — _L) and the E (_L) 
symmetry the superconducting spectra arc offset by 40 and 
by 20 counts, respectively. According to the strength of the 
CDW mode (labeled by C) Sample B and M have slightly 
diS'erent impurity concentrations. From Sooryakumar and 
Klein (1980). 



In the (incommensurate) CDW phase the symmetry representa- 
tions Aig etc. of the Dg^ point group do not apply any more. 
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FIG. 29 Raman spectra in A symmetry of 2i?— NbSe2 (sam- 
ple B) in various magnetic fields at 2 K. Prom Sooryakumar 
and Klein (1980). 



As can be seen in all panels, the CDW mode hardens 
below the superconducting transition. 

The difference between the two samples is apparently 
the impurity concentration affecting the strength of both 
the CDW and the gap mode. In fact, the CDW transi- 
tion can be suppressed by either pressure or an increas- 
ing number of defects, which may be quantified by the 
residual resistance ratio (Huntley and Frindt, 1974). In a 
systematic study of impurity effects, Sooryakumar et al. 
(1981) showed that the gap excitations go away along 
with the CDW mode while the superconducting transi- 
tion temperature is essentially unchanged. It is tempting 
to assume that the gap modes are directly coupled to 
the CDW mode and exist only along with it. This inter- 
pretation is supported by results obtained in a magnetic 
field (Fig. 29). Upon increasing the field the gap fea- 
ture in A symmetry is gradually suppressed while the 
CDW mode gains intensity leaving the energy integral 
over the Raman response x"{^) constant to within 7% 
(Sooryakumar and Klein, 1981). In E symmetry no clear 
sum rule could be found (Sooryakumar and Klein, 1981), 
and it is possible that some of the gap intensity appears 
independent of the CDW. 

Particularly the result in A symmetry (Fig. 29), where 
the gap mode gains intensity at the expense of the CDW 
mode, triggered the theoretical work to follow. The 
available data are not supportive of a sum rule in E sym- 



1* For convenience we again give the related references being dis- 
cussed thoroughly in the context of collective modes and in the 
previous paragraph: (Balseiro and Falicov, 1980; Browne and 
Levin, 1983; Klein and Dierker, 1984; Lei et al, 1985; Little- 
wood and Varma, 1981, 1982; Monien and Zawadowski, 1990; 
Tiitto and Zawadowski, 1992). We would like to draw the read- 
ers's attention also to the closely related Raman experiments in 
superfluid ''He (Greytak and Yan, 1969) and their theoretical 
description (Zawadowski et al., 1972). 



metry (sec Fig. 28 a), demonstrating similarities with the 
A15 compounds where also both electronic scattering and 
coupling to phonons was observed. In order to bring some 
light into the rather involved discussion, it is worthwhile 
to reconsider the influence of impurities (Sooryakumar 
et al., 1981). 

At first glance, the reaction to disorder points in the 
same direction as the results in magnetic fields. How- 
ever, defects not only suppress the formation of the CDW 
(Huntley and Frindt, 1974) but also, independently, re- 
duce the intensity close to 2A (Devereaux, 1992, 1993) 
while normally leaving the transition temperature Tc of a 
conventional s-wave superconductor unchanged (Ander- 
son, 1959)^^. For hvp/S < A, h/r the intensity at 2A 
is proportional to rA (Devereaux, 1992, 1993). In this 
respect Raman is just opposite to the optical conductiv- 
ity where the gap can l)e observed only if h/r is of the 
order of A or larger (Mattis and Bardeen, 1958). This 
implies that the Raman gap feature can be wiped out by 
impurities while Tc remains essentially constant; at the 
same time, though independently, the CDW transition is 
suppressed. Hence, it is possible that the gap features in 
2_ff— NbSe2 exist on their own as pair-breaking effect but 
the interaction with the CDW leads to a bound state. 

Most of the other CDW systems are not supercon- 
ducting, but show very interesting behavior around the 
transition to the charge-ordered phase. Some of them 
have been studied earlier using light scattering. Here we 
briefly discuss a recent study of the temperature-pressure 
phase diagram of the CDW state. 



2. lT-TiSe2 

In lT-TiSe2 a commensurate CDW is established be- 
low TcDW 200 K. The amplitude of the CDW couples 
to zone-boundary acoustic phonons which are folded to 
the center below Tcdw (Snow et al., 2003). Pronounced 
soft-mode behavior can be observed as a function of tem- 
perature. In the limit T — !■ two strong lines at 115 and 
75 cm~^ in Aig and Eg symmetry, respectively, dom- 
inate the low-energy spectra. By increasing the pres- 
sure the CDW state flrst stiffens along with the lattice 
then disappears rapidly in the pressure range of 5 to 
25 kbar. Above 25 kbar a quantum disordered (essen- 
tially isotropic) metallic or semi-metallic state is found 
although the Raman continuum typical for a metal is 
not reported. The quantum mechanical melting of the 
CDW order is in many ways similar to classical 2D melt- 
ing, with the appearance of crystalline and disordered 
CDW regimes, as well as an intermediate "soft" CDW 



In A15 compounds the high density of electronic states at Ep 
(partially responsible for the high Tc) depends sensitively on dis- 
order (Mattheiss and Weber, 1982). Hence, disorder reduces Tc 
fast as opposed to what one would expect from the Anderson 
theorem. 
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regime in which the CDW exhibits strong fluctuations 
and loses stiffness. Here, measurements on the develop- 
ment and polarization dependence of the electronic con- 
tinuum raises the possibility of following quantum critical 
behavior in other systems with competing orders. This 
is a promising direction for future studies. 



D. Kondo or Mixed-Valent Insulators 

All correlations discussed so far arc related to electron- 
lattice interactions in systems with screening lengths of 
the order of the interatomic spacing. With decreasing 
electronic density, new phenomena develop originating 
from the competition between kinetic and potential en- 
ergy of the conduction electrons. Well known examples 
are the Mott metal-insulator transition or Wigner crys- 
tallization. In either case, the material becomes an in- 
sulator at low temperature due to "immobilization" of 
electrons rather than an energy gap as in band insula- 
tors. 

Raman scattering on band insulators and semiconduc- 
tors has been well-documented, with a focus placed on 
high energy charge transfer excitations. Yet the devel- 
opment of the Raman response at low frequencies can 
in principle shed light on the development of electronic 
correlations with temperature and/or doping. 

Experiments by Nyhus et al. (1995a) on Kondo insulat- 
ing FeSi shown in Fig. 30 and Nyhus et al. (1995b, 1997a) 
on mixed- valent SmBe have indeed shown the transfer of 
spectral weight from low to high energies as the temper- 
ature is lowered into the insulating state. This "univer- 
sality" suggests that there is a common mechanism gov- 
erning the electronic transport in correlated insulators. 
As these materials are cooled, they all show a pilcup of 
spectral weight for moderate photon energy losses with 
a simultaneous reduction of the low-frequency spectral 
weight. This spectral weight transfer is slow at high tem- 
peratures and then rapidly increases as the temperature 
is lowered towards a putative quantum-critical point cor- 
responding to a metal-insulator transition. 

A characteristic energy appears which separates the 
region of intensity depletion from intensity enhancement 
with temperature. This characteristic frequency is es- 
sentially independent of temperature and is thiis called 
an isosbestic point in the spectra (as shown in Figure 
30). Finally, it is often observed that the range of fre- 
quency where the Raman response is reduced as T is 
lowered is an order of magnitude or more larger than the 
temperature at which the low-frequency spectral weight 
disappears. As discussed in the Section II. D. 4 these find- 
ings are consistent with the loss of low- frequency scatter- 
ing due to the thermal depletion of excited states. The 
channel-dependence has not yet been a focus of interest. 
If it were measured, important information concerning 
the evolution of the potential energy with doping could 
be inferred (Freericks et al, 2005). 
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FIG. 30 Temperature dependence of the Raman response 

measured in FeSi by Nyhus et al. (1995a). Here Ac denote 
the position of energy gap developing in the continuum at low 
temperatures, transferring spectral intensity into the peak at 
energy Aq . The sharp low energy features are phonons. 



E. Magnetic, charge and orbital ordering: Raman 
scattering in Eu-based compounds, Ruthenates, and the 
Manganites. 

There has been a great deal of interest in the relation- 
ship between diverse and exotic low-temperature phases 
of strongly correlated systems (Dagotto, 2005). In par- 
ticular, the manganites, ruthenates, Eu-oxides, and hexa- 
borides display charge ordered, paramagnetic insulating, 
and ferromagnetic metallic phases as a function of dop- 
ing, temperature, and/or temperature. Due to the com- 
plex interplay between spin, charge, and orbital degrees 
of freedom, these systems present a battleground where 
different ordered phases compete for primacy as knobs of 
the Hamiltonian are changed (Imada et al., 1998). 

The phase diagram of strongly correlated materials is 
more complex in systems having strong electron-lattice 
interactions as well as orbital ordering tendencies. Ra- 
man spectroscopy in systems such as the manganites and 
ruthenates have provided important information on the 
evolution of lattice, charge, and spin dynamics across 
phase boundaries. In many cases the transitions can be 
induced by applying pressure. 

While Raman scattering from phonons has tradition- 
ally provided important information concerning the de- 
velopment of locally or globally symmetry-broken states 
accompanied by the formation of static charge ordering, 
electronic Raman spectroscopy can be brought to bear on 
this problem as it can detect both fluctuating or static 
charge and/or spin ordering, and may reflect on the ten- 
dency toward orbital ordering as well. In addition, the 
polarization dependence can shed light on the types of ex- 
citations that are created in or near the ordered phases 
which may serve as signatures that certain interactions 
are more prominent than others. Thus Raman is a pow- 
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crfiil spectroscopic method by which the dynamics across 
quantum phase transitions can be investigated in corre- 
lated systems. 

Recent Raman studies on EuBg (Nyhus et al, 1997b), 
Eui-ajLaj-Bg, EuO (Snow et al, 2001), and Eui_2.Gda;0 
(Rho et ai, 2002), show that the metal-semiconductor 
transition in these materials is accompanied by distinct 
changes of the electronic continuum. A high-temperature 
paramagnetic semimetallic phase is well characterized by 
scattering from diffusive charge excitations which become 
less diffusive at lower temperatures when correlation ef- 
fects have not yet set in. However, the diffusive scatter- 
ing rate, when fit with Eq. (49), increases with decreasing 
temperature and scales with the magnetic susceptibility 
as the system begins to develop short-range magnetic or- 
der, typical of insulating behavior, as shown in Figures 19 
and 30. Finally, at low temperatures, a ferromagnetic 
metallic phase occurs, showing a flat continuum charac- 
teristic of a strongly correlated metal. The doping, po- 
larization, and magnetic field dependence of the spectra 
implies that the metal-semiconductor transition is pre- 
cipitated by the formation of bound magnetic polarons 
above the ferromagnetic ordering temperature. 

Concerning the ruthenates, recently Snow et al. (2002) 
and Rho et al. (2003) have studied the evolution of the 
spin and lattice dynamics through the pressure-tuned 
collapse of the antiferromagnetic Mott-like phases of 
Ca2Ru04, Ca3Ru207, and Ca2-a;Srj,Ru04 into a ferro- 
magnetic and possibly orbitally ordered metallic state at 
low temperatures. The studies have shown many char- 
acteristic features resulting from the interplay of strong 
electron-lattice and electron-electron interactions. These 
include (i) evidence of an increase of the electron-phonon 
interaction strength, (ii) an increased temperature de- 
pendence of the two-magnon energy and linewidth in 
the antiferromagnetic insulating phase, (iii) evidence of 
a charge gap development significantly below the metal- 
insulator transition (Tmi), and (iv) a hysteresis asso- 
ciated with the structural phase change. The latter two 
effects are indicative of a first-order metal-insulator tran- 
sition and a coexistence of metallic and insulating compo- 
nents for T < Tmi. The measurements have not yet been 
extended to probe the unconventional superconducting 
state at low temperatures. 

Raman measurements on cubic and layered manganites 
have been used to explore the interplay of spin, charge 
and orbital degrees of freedom. Yamamoto et al. (2000) 
and Romero et al. (2001) observe a suppression of the 
low-energy continuum at Big symmetry upon entering 
the charge- and orbital-ordered state. The interpreta- 
tion is not yet settled. Possible candidates arc spin den- 
sity or dynamical charge-orbital fluctuations (Yamamoto 
et al., 2000) or a collective CDW excitation Romero et al. 
(2001). The controversy can probably not be solved with- 
out a quantitative theoretical description. 

Bjornsson et al. (2000) have measured cubic 
Lai_j;Srj;Mn03 and via phonon lincshapc analysis have 
shown strong electron-phonon interactions involving lo- 



cal lattice distortions in the high temperature paramag- 
netic state which gradually vanish below the ferromag- 
netic transition. A broad hump in the electronic spectra 
develops at low temperatures in the metallic state around 
400 cm~^ for x=0.2, and weakens in intensity and shifts 
to higher energies for x=0.5. Although a polaronic peak 
would weaken with increasing carrier density, the shift to- 
wards higher energies with doping was interpreted rather 
as a low-energy plasma excitations within the ferromag- 
netic metallic phase. Since recent ARPES studies on 
the same compound (Mannella et al., 2004) and bi-layer 
manganite (Mannella et al., 2005) have revealed coexis- 
tence of quasiparticle and polaron features in the metallic 
phase and do not show evidence for low energy plasma 
excitations, more work is needed to clarify this issue. 

A comprehensive study via reflectance and Ra- 
man measurements on Pro.7Pbo.2iCao.o9Mn03, 
Lao.64Pbo.36Mn03, Lao.66Pbo.23Cao.iiMn03, and 
Pro.63Sro.37Mn03, by Yoon et al. (1998), have shown 
that the electronic continuum displays a change from 
diffusive polaronic peaks at high temperature to a 
flat featureless continuum, similar to that observed 
in the cuprates, in the low temperature ferromagnetic 
phase. A broad polaronic peak around 1200 cm~^ shifts 
to lower energies with increases doping indicative of 
weakened polaron binding energies. This is consistent 
with a crossover from a small-polaron-dominated regime 
at high temperatures to a large-polaron-dominated 
low-temperature regime. The low temperature phase 
also provided evidence for the coexistence of large and 
small polarons, also in agreement with the ARPES 
results. 

Recent work on Bii_a;Caa;Mn03 (x < 0.5) by Yoon 
et al. (2000) is of particular interest in connection with 
polarization studies for the evidence of charge ordering. 
As shown in Figure 31, Raman scattering offers a unique 
means of probing the unconventional spin and/or charge 
dynamics that arise when charge carriers are placed in 
the complex spin environment of a charge-ordered sys- 
tems. Using circularly (LL) in addition to linearly {xx 
and yy) polarized light anti-symmetric components of 
the Raman tensor were isolated. In cubic crystals such 
as Bii_2;Ca2;Mn03 they transform as the Tig irreducible 
representation (equivalent to A2g in tetragonal materi- 
als like the cuprates). Upon entering the charge-ordered 
phase a quasielastic scattering response appears with the 
Tig symmetry of the spin-chirality operator. Thus it was 
conjectured that the chiral excitations were signatures 
of either a chiral spin-liquid state associated with the 
Mn core spins, or of closed- loop charge motion caused 
by the constraining environment of the complex orbital 
and Neel textures. A possible path for charge motion is 
shown in Figure 31, emphasizing the circular nature of 
charge transfer. It is remarkable, that the spectral shape 
and the temperature variation of the characteristic en- 
ergy are quite similar to the low-energy response in the 
cuprates (sec section IV. D. 3) although a state with static 
order is entered in Bii_a;Caa;Mn03. It is interesting and 
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perhaps important how the two types of response are re- 
lated. 

Recently Saitoh et al. (2001) have performed Ra- 
man measurements on detwinned and orbitally ordered 
LaMnOs and have observed multiple peak structures 
which they interpret as orbital excitations or "orbitons" . 
While this has been challenged by Griininger et al. (2002) 
on the basis of selection rules, more recent measurements 
by Kriiger et al. (2004) related the peak features to sec- 
ond order phonon scattering activated via the Franck- 
Condon mechanism (Perebeinos and Allen, 2001). Even 
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FIG. 31 (a) Temperature dependence of the Raman spectra 
in the LL scattering geometry in Bio.igCao.siMnOa (b) Frac- 
tional change in the integrated quasiclastic Raman scattering 
intensity (50 - 350 cm^^) as a function of T/Tco for samples 
with charge-ordering temperatures Tco = 165 K (circles) and 
210 K (triangles). Inset: Fractional change in the quasielastic 
scattering amplitude A and fluctuation rate F as a function 
of T/Tco obtained via a fit of the data with Eq. (49). Lines 
are guides to the eye. (c) Example of a closed-loop path for 
charge motion in the charge-ordered phase (x = 0.5) which 
is not precluded by either the orbital configuration or by the 
spin environment. Filled and empty circles represent Mn^"*" 
and Mn'*"'" sites, respectively. From Yoon et al. (2000). 



if we could only scratch this interesting subject we hope 
we could demonstrate, that Raman measurements con- 
tinue to be of merit to study the interplay of strong cor- 
relations and electron-phonon coupling and the novel ex- 
citations which emerge in orbitally ordered systems. 

In this section wc have shown how the Raman spectra 
evolve as the degree of correlations increases in different 
materials. One common aspect is the non-trivial polar- 
ization and temperature dependence of the spectra which 
emerge in materials with increasing complexity. Finally, 
it was shown that Raman scattering can be applied to 
materials with varying degrees of competition between 
ordered states. Nowhere is this more apparent than in 
the results on the cuprates with high superconducting 
transition temperature. Due to the large amount of work 
devoted to these materials and the complexities of the is- 
sues raised, we split off the discussion of cuprates and 
related compounds into the following separate section. 



IV. HIGH TEMPERATURE SUPERCONDUCTING 
CUPRATES 

The story of the cuprates began with the discovery 
of superconductivity by Bednorz and Miiller (1986) in 
La2-xBaxCu04 with a; ~ 0.1 . . .0.2. Soon after the con- 
firmation by Cava et al. (1987), YBa2Cu307 with a T,. 
above 90 K was synthesized (Wu et al., 1987). These 
unexpected results led to further discoveries of supercon- 
ductivity in materials with a layered crystal structure 
with one or more Cu02 planes per unit cell (Fig. 32). 
Materials synthesis has yielded compounds with increas- 
ingly higher Tc, and major advances have been made in 
single-crystal growth methods, providing a diversity of 
samples. Independent of the material family and its re- 
spective maximal Tc, superconductivity exists for doping 
levels 0.05 <p < 0.27 and 0.12 < n < 0.18 with p and 
n the number of holes and electrons, respectively, per 
plaquette (Cu02 formula unit) (Loram and Tallon, 2001; 
Onose et al., 2001) as shown schematically in Figure 32. 
n = = p denotes half filling. 

It became clear soon thereafter that the cuprates are 
doped Mott insulators with strong electronic correlations 
dominating the entire phase diagram (Fig. 32). The 
emergence of high temperature superconducting phases 
in materials from which strong correlations yield anti- 
ferromagnetism and large departures from a canonical 
Fermi liquid theory has highlighted our limited under- 
standing of electronic correlations. While overdoped sys- 
tems seem to display a behavior of the resistivity close 
to at low temperature and well-defined quasiparticles 
in ARPES studies, strong deviations already occur for 



Recent work can be found in (Ando et al., 2004; Eisaki et al., 
2004; Erb et al., 1996a; Hardy et al., 1993; Liang et al, 2000, 
2002; Onose et al, 2001) 
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FIG. 32 Schematic phase diagram of the cuprates. On the hole-doped (p) side long-range antiferromagnetic (AF) order 
disappears rapidly. The maximal superconducting (SC) transition temperature Tc is strongly material dependent but always 
observed at p ~ 0.16. On the electron-doped (n) side the AF phase is more extended. Tc does not exceed 30 K at n ~ 0.14. 
T* represents the approximate crossover temperature to the pseudogap regime (Timusk and Statt, 1999). On the l.h.s and the 
r.h.s. the structures of prototypical Nd2-xCexCu04 and La2-xSrxCu04, respectively, are shown. The atoms are: Cu - red, O 
- blue, La,Sr and Nd,Ce - yellow. All cuprates are tetragonal or close to tetragonal with small material-dependent deviations. 
Nd2-xCexCu04 crystallizes in T' structure without O octahedra and a slightly shorter c-axis. La2-xSrxCu04 and all other 
hole-doped materials have octahedra which are cut into half for materials with more than one Cu02 plane. 



optimally doped systems and become increasingly pro- 
nounced as the antiferromagnetic phase is approached. 
At optimal doping, where Tc is maximal, the materials 
have already high normal state resistivities and, hence, 
are canonical examples that bad metals make good su- 
perconductors. Yet we still do not have a theoretical 
framework to understand why. 

Upon underdoping, the cuprates develop strong elec- 
tronic anisotropics in the Cu02 planes, which can be 
thought of as a signature of correlations. For this rea- 
son, momentum resolution is crucial for understanding 
the physical properties. ARPES has been very important 
from the beginning, revealing, among many other things, 
a strong k dependence of both the superconducting en- 
ergy gap and the pseudogap in the normal state (Cam- 
puzano tt ai, 2002; Damascelli et ai, 2003). As a more 
subtle effect, the quasiparticle weight and the incoher- 
ent part of the spectral function were observed to have 
a substantial variation with k and doping p, implying 
the importance of correlation effects and the existence of 
strongly momentum-dependent interactions (Damascelli 
et ai, 2003). More recent studies of scanning tunneling 
microscopy (STM) have indicated a presence of nanoscale 
disorder (McElroy et ai, 2005) in addition to the strong 
anisotropics identified from ARPES studies (Damascelli 
et ai, 2003). Both tools have combined to give insight of 
the tendencies these compounds have towards electronic 



ordering before the antiferromagnetic phase is reached. 
This demonstrates directly that the simultaneous under- 
standing of both single- and two-particle response func- 
tions is important. 

In this context, electronic Raman scattering has played 
a major role in characterizing the anisotropic dynamics 
of electrons across the phase diagram. These include the 
intense study of antiferromagnetism, where Raman mea- 
surements on the parent insulating cuprate compounds 
were the first to yield an estimate of the magnetic ex- 
change J from the energy of the two-magnon scattering 
peak in the Big channel (Lyons et ai, 1988). The dis- 
covery of the broad and fiat electronic continuum in the 
normal state of the cuprates close to optimal doping (Bo- 
zovic et ai, 1987) became the signature of the anoma- 
lous and strange metallic phase at high temperatures 
and spawned important ideas by Varma et al. (1989a) 
which honed in on the physics at play in these mate- 
rials. The polarization dependence of this background 
above (Slakey et ai, 1991; Staufer et ai, 1990) as well 
as below (Cooper et ai, 1988b; Hackl et ai, 1988b; 



A selection of references is (Damascelli et ai, 2003; Hanaguri 
et al, 2004; Hoffman et al, 2002; Howald et al, 2003; Vershinin 
et ai, 2004). 



39 



Slakey et ai, 1990b) visualized strong anisotropic inter- 
actions in these materials. Specifically, the observation 
of a polarization-dependent gap opening in the spectra 
below Tg was instrumental in solidifying the symmetry 
and orientation of the superconducting order parameter. 
Finally, more recent work has focused on elucidating the 
behavior of competing orders and complexity which come 
hand-in-hand with strong correlations. 

It is safe to say that no other system has bcicn studied 
so intensely via a bevy of experimental tools in recent 
years as the cuprates. This is certainly true for Raman 
spectroscopy, where the light scattering cross sections 
have been analyzed in a rich number and quality of ma- 
terials. In this section, we discuss Raman results on the 
cuprates and related materials, with an overview of re- 
lating findings to the physics uncovered in other systems. 
We show issues in which consensus has been reached and 
other issues which are controversial and require further 
analysis. 



A. From a Doped Mott Insulator to a Fermi Liquid 

The so-called parent compoimds of the cuprate super- 
conductors are antiferromagnetic Mott insulators. At 
half filling, charges are completely localized and excited 
states are separated by the Coulomb energy U prohibit- 
ing double occupancy in the Hubbard model. The 
ground state is a Heisenberg antiferromagnet which de- 
velops 3D long-range order below a Neel temperature of 
typically 300 K. The excitations are spin waves with a 
very small (in-plane) anisotropy gap at q and a 
nearest-neighbor exchange coupling J. The spin waves 
have been studied right at the beginning by Raman scat- 
tering, and J was determined from the two-magnon peak 
to be of order 100 meV.^''' Very early the existence of 
chiral spin excitations, i.e. excitations where the spins 
are rotated out of the easy (CUO2) plane, was demon- 
strated (Sulewski et ai, 1991). The resonance profile^° 
has been described in terms of the triple resonance theory 
of Chubukov and Frenkel (1995a, b). Very recently, one 
magnon excitations at q = were observed in the Neel 
state of La2-xSrxCu04 (Gozar et ai, 2004) and related 
to Dzyaloshinskii-Moriya and XY optical modes result- 
ing from a spin gap by Silva Neto and Benfatto (2005). 
The experiments arc the q = and hil manifesta- 
tion of canted spins which give also rise to scattering in 
the two-spin-flip channel at large energies in A2g symme- 
try (Sulewski et al., 1991). Here, a lot more information, 



The subject has been reviewed in detailed articles and books 
(Pulde, 1995; Gcbhard, 1997) 

(Blumberg et ai, 1996; Knoll et al, 1990; Lyons et al, 1988, 
1989; Riibhausen et al, 1997; Sugai et al, 1988; Sulewski et al, 
1990, 1991) 

A selection of references is (Blumberg et al, 1996; Knoll et al, 
1996; Lyons et al, 1988; Rubhausen et al, 1997). 



such as the Dzyaloshinskii-Moriya vector, the size of the 
anisotropy gap, and the spin-lattice interaction strength, 
could be derived (Gozar et al., 2005b). 

If the insulator is doped off half filling, antiferromag- 
netic long range order disappears quickly with hole dop- 
ing but survives up to higher electron doping levels (Fig- 
ure 32). Yet short ranged antifcTromagnetic correlations 
are not quenched even at high doing levels, as seen by the 
persistence of the two-magnon peak in Big Raman scat- 
tering experiments on both the hole-doped (Blumberg 
et al., 1994; Reznik et al., 1993; Rubhausen et al., 1999; 
Sugai et al, 2003) and the electron-doped side (Onose 
et al., 2004) of the phase diagram. This is summarized 
in Figure 33 for a number of compounds and shows how 
the two-magnon peak softens and broadens with doping, 
eventually merging into the continuum at higher dop- 
ing levels. Since the two-magnon intensity is dominated 
by a double spin flip of nearest neighbors, it can be ob- 
served even for small magnetic correlation lengths of 
the order of few lattice constants (see Figure 17). The 
two-magnon peak persists at least up to optimal dop- 
ing for hole-doped systems (Fig. 33). The position of 
the maximum and the peak intensity decrease by fac- 
tors of roughly 2 and 20, respectively, for < p < 0.16. 
Consistent with the Raman results, neutron scattering 
experiments in La2_xSrxCu04 revealed magnetic excita- 
tions in the complete superconducting range (Wakimoto 
et al, 2004). 

The broadening of the two-magnon response with dop- 
ing and temperature has been studied theoretically on 
clusters for the t — J model (Prelovsek and Jaklic, 1998) 
and by quantum Monte Carlo (QMC) techniques (Sand- 
vik et al., 1998), respectively. Knoll et al. (1990) analyzed 
their experiments at elevated temperatures in terms of 
spin-lattice coupling. While many features of the data 
at 1/2 filling for undoped cuprates can be captured by 
studies of Help Eq- (46), in general the evolution of 
the magnon lineshape with doping and temperature and 
its full polarization dependence is not very well under- 
stood. In particular, it is an unsettled issue, how, when 
a more metallic state develops for higher doping levels, 
the magnon line merges into the relatively featureless 
continuum (shown in Figures 33 and 34) that extends 
well-beyond all relevant energy scales such as hqvp, the 
superconducting energy gap A or the maximal phonon 
energy ftwD- In strongly overdoped yet superconduct- 
ing samples (0.20 < p < 0.27), the physical properties 
in the normal state are still not those of a conventional 
metal. The continuum itself displays significant polariza- 
tion and doping dependence, which will be addressed in 
section IV. D. 

For p > 0.05 (n > 0.12), superconductivity emerges 
and reaches a maximal Tc at approximately p = 0.16 (n = 
0.14). In the Raman spectra, superconductivity induced 
peaks emerge out of the flat continuum in the normal 
state, accompanied by spectral weight reorganization for 
temperatures below Tc, as shown for Bi2Sr2CaCu208+5 
in Figure 34. 
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In the following section we focus on the polarization 
dependence of the Raman response below Tc in both 
hole- and electron-doped cuprates. We summarize how 
symmetry arguments can be used to obtain the momen- 
tum dependence of electronic properties in general and, 
specifically, that of the energy gap A(k). The doping 
dependence in hole-doped systems is postponed and will 
be discussed in detail in section IV. C. 



B. Superconducting Energy Gap and Symmetry 
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FIG. 34 Raman spectra in the normal state (a) and supercon- 
ducting state (b) for Aig -\- B2g (left panel) and Big {+A2g) 
(right panel) orientations in Bi2Sr2CaCu208+5. From Ya- 
manaka et al. (1988). 



Early Raman results for the electronic continuum 
showed only little difference between the normal and the 
superconducting states (Lyons et al, 1987). We know 
now that the relatively high defect concentration in the 
first samples suppressed the structures related to the 
gap. The synthesis of flux-grown YBa2Cu307 with a 
sufficiently small number of defects (Kaiser et al., 1987; 
Schneemeyer et al., 1987) improved the situation rapidly, 
and a clear indication of the pair-breaking effect was ob- 
tained by Cooper et al. (1988a) on YBa2Cu307, showing 
the emergence of a peak and reorganized spectral weight 
occurring for temperatures below Tc as expected from 
theory. Soon thereafter a strong polarization dependence 
of the Raman spectra was observed (Fig. 34), which can 
be considered the first spectroscopic evidence of a gap 
anisotropy (Cooper et al., 1988b; Hackl et al., 1988b; 
Slakey et al., 1990a; Yamanaka et al., 1988). In contrast 
to conventional materials (see, e.g.. Fig. 24), there is no 
sharp onset of the scattering intensity at a threshold. As 
an explanation of the continuous increase of the scatter- 
ing intensity at small frequencies the possibility of nodes 
was discussed early (Hackl et al., 1988b; Monien and Za- 
wadowski, 1989). Originally shown in Bi2Sr2CaCu208+5 
and YBa2Cu307, the studies were quickly extended to 
other optimally hole doped cuprates. It was shown by 
Kang et al. (1996, 1997) that magnetic fields suppress 
the peaks, independently indicating their relationship to 
superconductivity. 
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I. Symmetry: Big and Big 

The polarization dependence of both the peak fre- 
quencies and the low-energy slopes of the supcTConduct- 
ing spectra had been a vexing problem for several years 
following the first observation of the reorganized spec- 
tral weight and the polarization-dependent response. As 
shown in Figure 34, the peak in the Big response devel- 
oped at roughly 30 percent higher energies than for 
or Aig polarization geometries, and the low frequency 
spectra rose as il^ for Big and linearly with CI for other 
channels. In addition, the temperature dependent deple- 
tion at low frequency shifts was faster in Big symmetry 
than in other channels. 

The satisfactory description emerged as soon as the 
momentum dependences of the Raman vertices of differ- 
ent symmetries /z, 7p(k), and of the energy gap A(k) 
were properly taken into account, as outlined in Section 

II. D. 6. As indicated by a few experiments before (Hardy 
et al, 1993: Slicn et al, 1993) and corroborated by many 
more later (Scalapino, 1995), a gap with (1^2 _y2 symme- 
try is the most compatible with the results in the cuprates 
(Dcvercaux et al., 1994a). 

The predictions of d— wave theory - such as fre- 
quency power-laws, temperature dependences, and rel- 
ative peak positions - were found to be consistent 
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FIG. 35 Fit to the Big and Big data on as-grown (Tc=86 K, 
top panel) and oxygen- annealed (Tc=79 K, bottom panel) 
Bi2Sr2CaCu208+« from Staufor et al. (1992). The inset is 
a log-log plot of the low frequency Big response, showing a 
cross-over from linear to cubic behavior at a characteristic fre- 
quency a;*/Ao=0.38 and 0.45 for the top and bottom panels, 
respectively. From Devereaux (1995b). 



with many optimally hole-doped compounds, such as 
YBa2Cu307 (Chen et ai, 1994a; Devereaux and Einzel, 
1995a), BisSraCaCuaOg+a (Blumberg et al, 1997a; De- 
vereaux et al, 1994a; Devereaux and Einzel, 1995a), 
La2-xSrxCu04 (Chen et al, 1994b), TlaBazCuOe 
(Blumberg et al., 1997b; Devereaux and Einzel, 1995a; 
Kang et al., 1996, 1997), BiaSraCuOe+a (Einzel and 
Hackl, 1996), TlaBaaCaCuaOs (Kang et al, 1997; Mak- 
simov et ai, 1990), TlaBaaCaaCusOiO (Hoffmann et ai, 
1994), and later HgBa2Ca2Cu308+5 (Sacuto et ai, 1998, 
2000), Bi2Sr2Ca2Cu30io-H5 (Limonov et al, 2002a) and 
HgBa2Cu04+5 (Gallais et ai, 2004) The identification of 
twice the maximal gap from the peak in Big channels 
yielded 2Amax — 8 ksTc for all compounds, well above 
the weak coupling values of 4.28 ksTc for d-wave pairing, 
indicating the strong-coupling nature of the pair state^^. 

An additional confirmation of d— wave pairing came 
from the impurity effects on the Raman spectra in the su- 
perconducting state (Devereaux, 1995b; Limonov et al, 
2002b; Misochko et al., 1999). Raman scattering, like 
other types of responses (ARPES, infrared and tunneling 
spectroscopy, NMR, etc.) couples only to the magnitude 
and not to the phase of the order parameter and, hence, 
cannot discriminate between a dx2-y2 and |rfj,2_y2| gap, 
i.e. the sign change of the c?-wave gap is not accessible. 
Impurities however could distinguish between energy gap 
which were conventional yet had accidental nodes from 
a pure d— wave energy gap (Borkowski and Hirschfeld, 
1994). While peaks in the Big channel are generally 
suppressed by impurities, the theory for Raman scatter- 
ing in disordered d— wave superconductors predicted a 
crossover from linear to cubic frequency dependence at 
a characteristic frequency cj* set by the impurity con- 
centration, while a true threshold would develop if the 
nodes were accidental (Devereaux, 1995b). Shown in 
Fig. 35, along with the improvement of sample quality 
and instrumentation, the crossover from a linear to a cu- 
bic frequency dependence was found to directly demon- 
strate the influence of impurities and to further solidify 
the d— wave picture for hole-doped cuprates, as cemented 
by SQUID measurements, reviewed in Refs. (Tsuei and 
Kirtley, 2000; Van Harlingen, 1995)^2. 



Due to this high ratio, we note however that in certain cases 
(Martinho et ai, 2004; Zeyher and Greco, 2002) the Big peak is 
thought not to be related directly to a 2A feature. However we 
note that 2A/A;bTc ~ 10 emerges from strong coupling d— wave 
treatments (Monthoux and Scalapino, 1994), and is consistent 
with the broadening of the Big signal observed near 2A. 
We emphasize again that disorder rapidly suppresses gap features 
in the Raman spectra (see section II. D. 6). Even in a d-wave su- 
perconductor where Tc reacts sensitively to impurities the peaJjs 
disappear long before Tc vanishes. Hence, the characterization 
of the samples is a central issue. Doping generally introduces 
defects along with spins or carriers such as Ni in the Cu02 plane 
or Sr in La2-xSrxCu04 and clusters of oxygen in YBa2Cu306+x 
(Erb et ai, 1996b; Pekker et al, 1991) suppressing the gap struc- 
tures rapidly. 
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Additional attributes have been intensely investigated, 
such as the agreement of the effective mass approxima- 
tion calculated within LDA for d— wave superconductors 
in YBa2Cu307 (Strohm and Cardona, 1997), effects of 
orthorhombicity and mixing of s— wave components al- 
lowed by symmetry (Nemetschek et ai, 1998; Strohm 
and Cardona, 1997), as well as the issue of the Aig peak 
seen in the superconducting state (Krantz and Cardona, 
1994), which is among the complicated though eventually 
crucial problems on the way to a better understanding of 
the cuprates. 



2. The Aig problem - Zn, Ni and Pressure 

While the Big and the B2g spectra are in reasonable 
agreement with the predictions on the basis of d-wave 
pairing various problems arose in Aig symmetry which 
slowed down the acceptance of the model in the begin- 
ning. One of the objections was the strong intensity 
found in the Aig channel. As discussed in Section II. D. 6, 
the Aig response is complicated due to the backflow 
terms which are as singular as the bare terms for frequen- 
cies at the gap edge, leading to cancellations of diverging 
intensities. For both a cylindrical Fermi surface and a 
tight-binding band structure (Devereaux et ai, 1994b; 
Krantz and Cardona, 1994), the effective mass approxi- 
mation predicts a suppression of the Aig intensity in com- 
parison to other channels, in contrast to what is found 
experimentally, as shown in Figures 34 and 38. Since the 
effective mass approximation is of questionable utility in 
highly correlated systems like the cuprates, this objection 
was not as serious, as an overall comparison of spectral 
intensities is not possible without detailed knowledge of 
the Raman scattering matrix elements given in Eq. (36). 
Good fits to the data of the Aig, B2g and Big symmetries 
(with the overall intensities as free parameters) could be 
obtained (Devereaux and Einzel, 1995a). However, a sen- 
sitivity to band structure and higher harmonics of the 
energy gap and Raman vertices found in the calculations 
implied a similar sensitivity to details of the materials 
which was not observed in the data. 

Further, it was argued by Krantz and Cardona (1994) 
that scattering in a multi-band system - such as in the 
Cu02 bilayer - would yield a sharp intensity at twice the 
gap edge, which was again inconsistent with the data on 
single- and double-layer systems available at that time. 

The issue of multiband scattering was partially re- 
solved by Devereaux et al. (1996). It was found that if 
the energy gaps were equal on at least two sheets of the 
Fermi surface split by bi-layer hopping, a diverging in- 
tensity would be possible. However the intensity is only 
proportional to the difference of the individual Raman 
vertices of the bands. This is a qualitative reason why 
the multi-band case would give peaks only under quite 
special conditions, implying that the fits obtained from 
the single band case are still valid. Yet, it does not sat- 
isfactorily explain the sensitivity of the Aig response to 



band and energy gap anisotropy factors. 

A recent undertaking to understand the origin of the 
Aig intensity has involved the response of the peak to 
partial replacement (up to a few percent) of Cu by Zn 
and/or Ni in YBajCusOy-i (Gallais et al, 2002; Le 
Tacon et ai, 2005; Limonov et ai, 2002b; Martinho et ai, 
2004). Martinho et al. (2004) found that the intensity of 
the Aig peak was insensitive to either Ni or Zn doping, 
while the peak in Big is suppressed by Zn, in agreement 
with earlier results (Limonov et ai, 2002b). As shown 
in Figure 36, Gallais et al. (2002) found that the peak 
position in Aig was reduced by Ni impurities and made 
the important observation that it followed that of the 
magnetic spin resonance mode (Sidis et ai, 2004). Le 
Tacon et al. (2005) found that Zn doping increased the 
low frequency spectral weight and suppressed the Big 
peak without changing its position. It is not clear why 
the position of the Big peak does not change while Tc 
decreases. Possible explanations are an inhomogeneous 
distribution of the impurities or an accidental cancella- 
tion effect between the decreasing 2A(Tc) and an increase 
of the peak energy fi^^^^^^ because of defects (Devereaux, 
1995b). The shift of the Aig peak with Zn argued for a 
two-component picture including a contribution from a 
collective mode such as the tt resonance in the spin wave 
spectrum (Hinkov et ai, 2004). 

Two approaches concerning the presence of a collective 
mode in the spectra superimposed upon a well-screened 
background have been put forward recently by Venturini 
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FIG. 36 Raman and neutron-resonance peak ener- 
gies as a function of the critical temperature Tc in 
YBa2(Cui_xZnx)3 06.95 for x = 0, 0.01, and 0.03. The hor- 
izontal line for the Big peak positions is just a guide to the 
eye, while the Aig peak and the neutron resonance are fit- 
ted by a straight line representing SkeTc. From Gallais et al. 
(2002). 
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et al. (2000) and Zeyher and Greco (2002). Vcnturini 
et al. (2000) showed that a collective mode, uniquely ap- 
pearing in Aig and not other channels due to symmetry, 
emerges due to coupling to the 41 meV spin resonance 
mode seen in neutron scattering measurements in a num- 
ber of compounds (Sidis et al., 2004). Zeyher and Greco 
(2002) argued that the peak in the B-^g channel may be 
identified as a collective mode split off from 2A (Blum- 
berg et al., 1997a) as a consequence of the simultaneous 
presence of long-range d-CDW and d-superconducting or- 
der and that the Aig peak may be related to a supercon- 
ducting amplitude mode. In either case the sensitivity to 
parameters characterizing the anisotropics of energy gap, 
band structure, and Raman vertices were not present, yet 
it is still unclear which, if either, are able to explain the 
Aig peak in the superconducting state. 

Yet a remarkable set of experiments by Goncharov and 
Struzhkin (2003) plotted in Figure 37 highlights some 
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FIG. 37 Pressure dependences of the superconducting peaks 
in YBa2Cu306.9B for Big (full circles) and A\g (open circles) 
symmetries. The position of the B\g phonon is also indicated 

(triangles) . The inset shows a comparison to the Raman data 
on Bi2Sr2CaCu208+i samples at ambient pressure with dif- 
ferent doping levels (Kcndziora and Rosenberg, 1995). PYom 
Goncharov and Struzhkin (2003). Reproduced with permis- 
sion from Wiley © 2003. 



shortcomings in the above mentioned scenarios. Upon 
pressure the peak in the Big channel and T,, move down- 
ward in a similar fashion as upon overdoping (discussed 
in detail in Section IV. C), and the Big phonon hardens 
as expected, while the peak position in Aig remains re- 
markably constant. This implies that the Big peak is 
intimately tied to superconductivity and the Aig peak 
has a substantial contribution from a channel which is 
relatively insensitive to pressure and changes in super- 
conducting properties. In the discussions of possible can- 
didates for the mode highlighted above, the pressure sen- 
sitivity of spin, charge, and superconducting order would 
all be expected to be large. Thus the origin of this peak 
remains presently unclear. One can speculate that it may 
be related to possible phonon modes involving mixed Bi- 
O and Ba-0 which have been indicated to be less sensitive 
to pressure, but clearly further work is needed to clarify 
this matter (Goncharov and Struzhkin, 2003). 

In multi-layer compounds with more than 2 adjacent 
Cu02 planes, a strong Aig peak in the superconducting 
state occurs at roughly the peak frequency of the Big 
channel, as shown in Figure 38 for Bi2Sr2Ca2Cu30io+5. 
In addition, strong phonon resonances were found here 
and for multilayer Hg compounds (Hadjiev et al, 1998), 
where some Aig c-axis phonons shift by as much as 20 
cm""'^ at the superconducting transition, more than in 
any other high Tc compound, implying a strong renor- 
malization of the Aig continuum below Tc. 

Recently, Munzar and Cardona (2003) argued that a 
c— axis plasma mode would be expected to be Raman ac- 
tive in cuprate materials with more than 2 Cu02 planes, 
giving rise to an additional contribution in the Aig chan- 
nel only due to mass fluctuations with opposite sign on 
different Cu02 layers. The position of the plasma reso- 
nance was predicted to lie in a frequency range close to 
the continuum peak, as shown in Figure 38. 



3. Resonance Effects 

As another experimental knob to turn, resonance stud- 
ies of the peaks developing in the superconducting state 
provide information on the character and the interactions 
of the quasiparticles forming the condensate. For this 
reason, studies of the gap feature as a function of the 
incoming photon energy were pursued. Effects in the 
energy range of the gap were observed in various com- 
pounds with one (Blumberg et al., 2002; Kang et al., 
1996) more than one Cu02 layer (Budelmann et al., 2005; 
Hadjiev et al., 1998; Limonov et al., 2002a; Riibhausen 
et al, 1999: Sacuto et al, 1998). There are two dis- 
tinct, though not necessarily independent, effects: (i) 
the spectral weight and the shape of some phonons of 
predominantly Aig symmetry change more or less dra- 
matically below Tc (Hadjiev et al., 1998; Limonov et al., 
2002a), and (ii) the intensity of the electronic continuum 
is amplified in the vicinity of 2Amax (Fig. 38) (Blumberg 
et al., 2002; Limonov et al, 2002a; Riibhausen et al.. 
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FIG. 38 Resonance enhancement of the of the Aig Raman 
response in Bi2Sr2Ca2Cu30io+« at 10 K. From Limonov et al. 
(2002a). 



1999). The resonance effects typically occur for excita- 
tion energies of 2 ± 0.2 eV and are particularly strong 
in compounds with 3 or 4 Cu02 layers (Limonov et at, 
2002a). ^■^ Importantly, the normal state is only weakly 
affected. Hence, although the resonance energy is close to 
the 1.6 eV absorption edge (Singley et ai, 2001; Uchida 
et al, 1991), constraints are imposed on explanations 
in terms of interband transitions between the lower and 
the upper Hubbard band (c.f., Eq. (44) and Shastry and 
Shraiman (1990)), since the width of the resonance is 
much larger than 2Amax- At least in the triple-layer com- 
pounds, an additional channel can be opened up below 
Tc due to the c-axis plasma resonance (Munzar and Car- 
dona, 2003) which may interact with low-energy phonons 
and renormalize their shape and intensity substantially. 
The existence of a resonating continuum in double-layer 
compounds is still a matter of debate. Some authors 
find the pair-breaking peaks to resonate (Blumberg et al., 
1997a; Budelmann et ai, 2005; Riibhausen et ai, 1999), 
while others do not (Limonov et at, 2002a; Venturini 
et ai, 2002c). 

Given this background, the strong resonance effects 
in electron-doped Nd2-xCexCu04 were a quite interest- 
ing and somewhat unexpected feature (Blumberg et at, 
2002). While the i?2g response in the superconducting 
state is strongly enhanced toward the red, neither the 
Big spectra below Tc nor the normal state spectra in 
general are particularly sensitive to the energy of the ex- 
citing light. An explanation in terms of Hubbard physics 
is certainly tempting but, for the symmetry and tem- 
perature dependence, not completely exhaustive. In our 
opinion, the subject needs further experimental and theo- 



retical clarification before arriving at a level of predictive 
power. 

After the observation of resonance effects in 
Nd2-xCexCu04 and along with an improved mate- 
rial quality the number of studies in electron-doped 
systems increased continuously and facilitated several 
interesting insights which will be summarized in the 
following paragraph. 



4. Electron versus Hole-Doped Materials 

Some cuprates such as Nd2-xCexCu04 crystallize in 
the T' structure which is characterized by missing oxygen 
octahedra (see Figure 32) and, as a consequence, a short 
c-axis (Tokura et at, 1989). The maximal Tc of approx- 
imately 30 K is obtained in thin films of La2_xCexCu04 
(Naito and Hepp, 2000). None of the electron-doped 
cuprates is completely ordered, and oxygen appearing in 
an apex position is the main defect (Radaelli et ai, 1994) 
even in superconducting samples. 

Phonons (Heyen et ai, 1991) and crystal-field excita- 
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It is interesting to note that the resonance of the two-magnon 
peak is close to the charge transfer gap of 2.5 eV or higher (Blum- 
berg et al, 1997a; Knoll et al, 1996; Rubhausen et al, 1997). 
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FIG. 39 Doping dependence of the low energy electronic Ra- 
man response of Pr2-xCexCu04 single crystals and thin fihns 
for B2g, Big and Aig channels obtained with 647 nm exci- 
tation. The columns are arranged from left to right in order 
of increasing cerium doping. Abbreviations UND, OPT and 
OVD refer to under-doped, optimally doped and over-doped 
samples, respectively. The normal state response (light/red) 
measured just above the respective Tc is decomposed for the 
B2g and Big channels into a Drude-Iike component with a 
constant carrier lifetime (green dotted line) and an extended 
continuum (yellow dotted line). Superconducting spectra 
(dark/blue) are taken at T 4 K. For the OPT crystal a 
low-frequency lj^ power law is shown at Big symmetry. From 
Qazilbash et al. (2005). 
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tions (Jandl et al, 1993) in Nd2-xCcxCu04 were stud- 
ied soon after the discovery (Tokura et ai, 1989) and 
explained thoroughly. An exception was the Aig line at 
590 cm~^, assigned by Heyen et al. (1991) as a locahzed 
mode of interstitial oxygen in apex position (cf. Fig. 32). 
Later on, Onose et al. (1999) could indeed demonstrate 
that the mode is suppressed after annealing the samples 
in reducing atmosphere in order to induce or enhance Tc. 
The first electronic Raman spectra in the superconduct- 
ing state of slightly overdoped Ndi.84Ceo.i6Cu04 showed 
a small gap anisotropy (Stadlober et ai, 1995). The ratio 
4.1 < 2A/kBTc < 4.9 is similar to that of strong-coupling 
conventional superconductors like Pb, Nb or NbsSn (see 
Table 11). The temperature dependence of the gap is 
relatively close to the BCS prediction (Blumberg et ai, 
2002; Stadlober et al., 1995). This motivated an inter- 
pretation in terms of an anisotropic s-wave gap, yielding 
a reasonable description of both the shape and the posi- 
tions of the Big and pair-breaking peaks (Stadlober 
et al., 1995). 

Similar spectral shapes were also found at other dop- 
ing levels for both Nd2-xCexCu04 and Pr2_xCexCu04, 
as shown in Figure 39. The low-energy sides of the peaks 
were found to be almost doping independent and closer 
to those of hole-doped cuprates than to those in conven- 
tional superconductors when samples and instrumenta- 
tion facilitated improved measurements close to f2 = 0. 

Recent ARPES (Armitage et al, 2001; Matsui et al., 
2005) and interferometric experiments (Chesca et at, 
2003; Tsuei and Kirtley, 2000) provided evidence of a 
d-type gap bringing antiferromagnetic spin fluctuations 
back into play as a possible coupling mechanism. Since 
the diameter of the Fermi surface encircling (tt, tt) is 
smaller here than in hole-doped systems, the antiferro- 
magnetic ordering vector Qaf = (t^iTt) connects spots 
close to (7r/2,7r/2) rather than in the vicinity of (7r,0) 
for p doping and the gap is expected to have a maximum 
at the hot spots and to decrease towards (tt, 0). 

On this basis Blumberg et al. (2002) proposed a novel 
type of non-monotonic d-wave gap with maxima only ap- 
proximately 15° away from the diagonal. An explicit 
calculation using the suggested A(k) (Blumberg et al., 
2002) in a one-band model approximately reproduces the 
overall line shapes but reveals discrepancies in the peak 
positions (Venturini et al., 2003) (see also Blumberg et al. 
(2003)). Better agreement between experiment and the- 
ory can be obtained with a monotonic d— wave gap in a 
two band picture (Liu et al., 2005) for which evidence has 
been found by magnetotransport (Fournier et al, 1997). 

Four remarks are to be considered: (i) In ARPES, the 
gap maximum is found approximately at the hot spot 
where a pseudogap is observed above Tc (Matsui et al, 
2005). (ii) No quasiparticle peaks indicating coherence in 
the superconducting state are resolved in ARPES (Ar- 
mftage et al, 2001; Matsui et al, 2005). (iii) Twice 
the gap energy observed by ARPES is ~ 40% smaller 
than that derived from the Raman spectra (Matsui et al, 
2005). (iv) Phase-sensitive experiments (Alff et al, 1999) 



and results on the magnetic penetration depth (Kokales 
et al, 2000; Skinta et al, 2002a) are supportive of an 
s-type gap. 

Hence in spite of mounting evidence of d-wave pair- 
ing, explanations for the symmetry dependence of the 
Raman spectra and several other experiments are still 
missing. A possible explanation could lie in a crossover 
from d to s pairing upon increasing doping (Skinta et al, 
2002b). However, Raman data from differently doped 
Pr2_xCexCu04 (see Figure 39) do not show a variation 
of the lineshape and sufficiently strong symmetry depen- 
dence of the pair-breaking features in the proper doping 
range (Qazilbash et al, 2005). 

In the hole-doped systems, on the other hand, strong 
variations of the lineshapes at Aig and Big symmetry arc 
found although there is little doubt about the persistence 
of the (dominant) d-wave nature of the superconduct- 
ing gap in the entire phase diagram (Tsuei and Kirtley, 
2000). The doping and polarization dependence of the 
Raman spectra in superconducting p-type cuprates will 
be the subject of the next section. 



C. Superconducting Gap: Doping Dependence 

One of the early Raman experiments on doping effects 
in overdoped Bi2Sr2CaCu208-i-5 showed the energy of the 
pair-breaking peak in Big symmetry to decrease much 
faster than those at the other symmetries and Tc (Staufer 
et al, 1992). The first systematic study was performed by 
Kendziora and Rosenberg (1995) on Bi2Sr2CaCu208+5, 
shown in the inset of Figure 37. The doping level in this 
material can be varied continuously and reliably between 
0.15 < p < 0.23 by changing the oxygen concentration 
(Kendziora et al, 1993; Triscone et al, 1991). At lower 
oxygen doping, the structure is possibly metastable or 
unstable, at higher doping the oxygen diffuses out even 
at room temperature. Underdoping is better achieved 
by replacing Ca with Y. As an essential result, the pair- 
breaking peaks at Aig, Big, and B2g symmetry arc found 
to depend in different ways on doping and, consequently, 
on Tc- This doping dependence is shown in Figure 40 
for recent measurements on Bi2Sr2Cai_xYxCu208+5 by 
Sugai et al (2003). Most remarkably, the Big pair- 
breaking peak is proportional to (1 — p) rather than Tc in 
the doping range indicated, but fades in intensity for un- 
derdoped materials. The peak energies in B2g symmetry 
however scale more or less with the transition temper- 
ature. For p < 0.15 the Big peak becomes very weak, 
yet a clear 2A peak survives in the B^g channel for all 
doping levels. 

In the decade to follow there were numer- 
ous studies on differently doped cuprates, in- 
cluding Bi2Sr2CaCu208+6,24 YBa2Cu306+x,2^ 



(Blumberg et al., 1997a; Budelmann et al, 2005; Hackl et al.. 
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La2_xSrxCu04,2** HgBa2Cu04,27 HgBasCazCuaOs,^^ 
TlaBaaCuOe,^^ TlaBaaCaCuzOs,^^ and 

Tl2Ba2Ca2Cu30io A compilation of experimental 
results is shown in Fig. 41. The scatter of the data 
points partially reflects the development of the sample 
quality. However, there are also discrepancies between 
the results which are related to the interpretation of the 
experiments. 
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FIG. 40 Doping dependence of the B\g and B^g Raman spec- 
tra in Bi-2212 at 5 K and 100 K. From Sugai et al. (2003). 
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We first summarize the generally accepted features 
close to and above optimal doping: 

• at p ~ 0.16 the peaks in the three Raman active 
symmetries are in relative positions expected for 
d wave pairing, the ratio 2Aniax/fcs?c is approxi- 
mately 8,^^ 

• in the overdoped range, p > 0.16, the Big 

peak frequency decreases faster than Tc obeying 
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FIG. 41 Compilation of the position Slpeak of the peak in 
the Raman response in the superconducting state normal- 
ized to T^^^. Big (open diamonds) and B2g (squares) ori- 
entations are shown for a variety of compounds. The re- 
spective references are color coded as follows: Bi-2212 [red 
(Venturini et al., 2002c), green (Sugai and Hosokawa, 2000; 
Sugai et al, 2003), cyan (Kendziora and Rosenberg, 1995), 
yellow (Blumberg et ai, 1997a; Liu et al, 1999), pink (Hewitt 
and Irwin, 2002), and gold (Masui et al, 2003)], LSCO [blue 
(Sugai et al, 2003)], Bi-2223 [grey (Masui et al, 2003)], Tl- 
2201 [purple (Gasparov et ai, 1997, 1998; Kang et ai, 1997; 
Nemetschek et al, 1993)], Tl-2223 [olive (Hoffmann et al, 
1994; Stadlober et ai, 1994)], and Hg-1201 [black (Gallais et 
ai, 2005)]. The dashed black line is the interpolation formula 



6[1 - 82.6(p - 0.16)^]. For comparison, results 



for twice the maximal leading edge gap of ARPES (circles) 
(Gampuzano et al, 2002) and for the peak-to-peak energy in 
the tunneling density of states (triangles) (Zasadzinski et at, 
2002) are included. 
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• whenever the i?2g peak can be observed its maxi- 
mum follows Tc, Op^'^k oc Tc,^^ 

• the Aig peak frequency follows cither the magnetic 
(tt, tt) mode^'' or, in the case of resonantly enhanced 
light scattering, ^p^l^.^^ 

The gap close to the node, which is projected out in 
symmetry, is also found by penetration depth, low 
bias STM, and the Ncrnst effect to more or less follow the 
transition temperature for all doping levels (Deutscher, 
1999, 2005; Panagopoulos and Xiang, 1998; Xu et al., 
2000). However it is inconsistent with peak-to-peak 
measurements of tunneling density of states (Zasadzinski 
et al, 2002), the energy of the (tt, 0) peak in the spectral 
function(Campuzano et al, 2002), and thermal conduc- 
tivity (Sutherland et al., 2005), which when interpreted 
in terms of li— wave quasiparticle picture (Durst and Lee, 
2000) indicate that the gap energy continues to increase 
with decreased doping. This highlights one of the major 
issues in the cuprates of whether the superconducting 
pairing energy rises with underdoping or follows T^, and 
this certainly requires further study. 

In many expca'iments, a second energy scale is observed 
which varies approximately as po — p, similar to fipJak- 
The majority of authors (Timusk and Statt, 1999) call it 
a pseudogap A* opening below the crossover temperature 
T* (see Fig. 32) with A* oc T*. The understanding of the 
pseudogap is a matter of intense research at present. A 
possible relation to the Big Raman data will be discussed 
below in section IV. C. 

In the underdoped range, the interpretation is more 
controversial. The main issues are whether or not the 
B2g pair-breaking peak can be observed at all doping 
levels and how the B\g spectra evolve for p < 0.16. 

The i?2g problem seems to converge with the improve- 
ment of the sample qiiality. The cleaner the samples, the 
clearer the i?2g pair-breaking peak, as expected theoreti- 
cally (Devereaux, 1995b). The problem can be visualized 
in YBa2Cu306-i-x where oxygen tends to cluster and to 
form pinning centers if the Cu-0 chains are not com- 
pletely filled (Erb et al., 1996b). However, in addition to 
fully oxygenated YBa2Cu307, partially ordered phases 
exist for x = 0.5 with Tc ~ 60 K and x — 0.353 with 
Tc ^ (Liang et al., 2000, 2002) where every second 
and third chain is filled, respectively. Therefore, the max- 
ima are pronounced in YBa2Cu306.5o and YBa2Cu306.98 
(Opel et al., 2000), while the peak is smeared out at op- 
timal doping, YBa2Cu3 06.93, and practically disappears 
shghtly below (Sugai et al., 2003). The problem of oxy- 
gen clustering exists also in Bi2Sr2CaCu208+5, but the 



resulting potentials are weaker and essentially doping in- 
dependent, as can be seen directly from the comparison of 
different overdoped samples (Opel et al, 2000; Venturini 
et al, 2002c) exhibiting nearly constant intensity of the 
-629 pair-breaking features. In a similar fashion, partial 
replacement of Ca by Y does not create a strong impurity 
potential either. This was directly shown by electron spin 
resonance (ESR) (Janossy et al, 2003) in YBa2Cu306.i, 
where Ca in place of Y is a very weak impurity which 
does not localize carriers even at low temperature and 
doping. We conclude that Y and O doping can be used 
more or less simultaneously in Bi2Sr2CaCu2 08+5 as long 
as Y is distributed statistically. The direct comparison 
of the data presented in Fig. 40 with those of Opel et al. 
(2000) and Venturini et al. (2002c) indeed shows that 
the B2g pair-breaking peaks have doping and dopant (O, 
Y) independent intensities if a high sample quality can 
be maintained. Under these conditions, the i?2g maxima 
seem to exist at all doping levels and apparently scale 
with Tc- Since in YBa2Cu306+x and La2-xSrxCu04 dop- 
ing changes both carrier concentration p and mean free 
path the pair breaking peaks appear and disappear 
depending on the doping-dependent order. 

It is important to remember these considerations 
for the analysis of the B\g data. This is partic- 
ularly important when intensity issues are discussed. 
Therefore, we focus first on Y-under- and 0-overdoped 
Bi2Sr2CaCu208+5 (Fig. 40) to determine the intensity 
evolution of the Big pair-breaking structure. It becomes 
very weak right below optimal doping while no significant 
weakening of the i?2g structures is observed. The gradual 
suppression of the Big coherence peaks in Raman goes 
along with the disappearance of the 2A peaks in scanning 
tunneling microscopy (McElroy et al, 2003). This is cor- 
roborated by results in partially ordered YBa2Cu306.5, 
in Lai.9Sro.ioCu04, HgBa2Cu04, and HgBa2Ca2Cu308, 
where the i?2g structures are well resolved, while in Big 
symmetry, no pair-breaking effect could be found (Gal- 
lais et al, 2005; Naeini et al, 1999; Opel et al, 2000; 
Venturini et al, 2002b). 

When the peak is observed in Big in the underdoped 
region, as shown in Figure 41, its position does not track 

Tc and lies within the scatter of points or slightly below 
the values of 2A generally obtained via ARPES (Cam- 
puzano et al, 2002) and peak-to-peak positions in tun- 
neling measurements (Zasadzinski et al, 2002)."^^ It has 
been argued by Chubukov et al (1999, 2006) and Zeyher 
and Greco (2002) that the Big peak in Raman measure- 
ments in underdoped compounds may be due to a collec- 
tive mode appearing below Tc from either spin-coupling 



(Gallais et at, 2005; Kcndziora and Rosenberg, 1995; Misochko 
et al, 1999; Opel et al, 2000; Venturini et al, 2002b,c) 
(Gallais et al, 2002; Le Tacon et al, 2005) 
(Limonov et al, 2002a) 



We point out however that there is considerable uneertainty in 
determining the gap from ARPES, as both the leading edge and 
the peak of the spectrum have been used. In either ease, the anti- 
nodal spectral function is very broad in underdoped systems and 
identifying an energy scale from a broad feature is not without 
uncertainty. 
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or d— CDW order, respectively, split off from twice the 
gap maximum 2A. Yet this interpretation is inconsis- 
tent with the doping behavior of the pairing gap deter- 
mined from B2g symmetry, as discussed above. In addi- 
tion, the findings of the doping behavior of the anti- nodal 
quasiparticles in the normal state (listed below in Sec- 
tion IV. C) indicate that anti-nodal quasiparticles have 
become incoherent already above optimal doping. There- 
fore it is not intuitively obvious how a propagating mode 
could emerge deep in a region of incoherence. An alter- 
native scenario is that the peak observed in Big channels 
is strongly altered by incoherence and is a measure of the 
binding energy of localized anti-nodal electrons, such as 
via the formation of resonance- valence bond singlets or 
small polarons. Yet why a peak should appear remains 
to be understood. In any case, this issue highlights an- 
other vexing problem in the cuprates: how involved are 
the anti-nodal electrons in superconductivity for under- 
doped systems. 

There were reports about Big pair-breaking peaks at 
low doping p <0.1 (Blumberg et al., 1997a; Slakey et al., 
1990b) appearing already above Tc, at a doping indepen- 
dent position of approximately 600 cm~^. In the dis- 
cussion of the existence of preformed pairs in the pseu- 
dogap state above Tc and of collective modes inside the 
gap (see Section II. D. 7) that are indicative of the pair- 
ing potential, this result has obviously some importance. 
However, in the vast majority of the experiments, the ob- 
servation could not be confirmed. Further complication 
comes from the existence of an oxygen impurity mode in 
close vicinity (Hewitt et al, 1999; Quilty et al, 1998). 
Future studies have to clarify this issue. 

In conclusion, Raman scattering experiments reveal 
two energy scales in the superconducting state which 
have a different dependence on doping: while the peaks 
in B2g symmetry follow the transition temperature, the 
maximum energy of those in Big symmetry decreases as 
(Po ~ p)- In the majority of the experiments, the co- 
herence peaks in Big symmetry are found to fade away 
rapidly for p < 0.16. The origin of the features observed 
at small doping remains controversial. 



D. Normal State: Dichotomy of Nodal and Anti-Nodal 
Electrons 

The studies of the symmetry and doping dependence of 
the Raman spectra in the superconducting state are sug- 
gestive of interactions with a pronounced structure in mo- 
mentum space which, in addition, vary with doping (see, 
e.g.. Fig. 41). One of the key questions is therefore how 
the interactions renormalize the quasiparticles and their 
response functions in the normal state. "^^ Owing to the 



The properties of the normal state have been studied and re- 
viewed extensively in the last two decades. Recent reviews are 



crystal structure of the cuprates, the anisotropy of the in- 
plane to the out-of-plane transport translates into a mo- 
mentum dependence of the in-plane transport properties 
(Devereaux, 2003b; Forro, 1993; Turlakov and Leggett, 
2001). Hence, Raman scattering can substantially sup- 
plement optical conductivity measurements above Tc. An 
anisotropy of the normal-state Raman spectra was actu- 
ally observed soon after the gap anisotropy.^^ 



1. Unconventional Metal-Insulator Transition 

The momentum dependence of the transport as mea- 
sured by Raman scattering is indeed dramatic. Gener- 
ally, the nodal region along the diagonal of the Brillouin 
zone {B2g symmetry) remain essentially unchanged at 
all doping levels, while the (tt, 0) regions (Big symme- 
try) suffer an overall loss of oscillator strength by ap- 
proximately an order of magnitude if the doping level is 
reduced from 0.22 to 0.10 (Fig. 42). When this doping 
dependence was first observed, Katsufuji et al. (1994) re- 
alized that the variation of the Big spectra with p cannot 
be explained by only changing the filling of a rigid band 
where the ratio of the Big to the B2g scattering intensi- 
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FIG. 42 Direct comparison of the low-energy Big and Big Rar 
man response functions measured at temperatures indicated 
(all below Tc) in La2-a;Sra;Cu04 for different values of Sr dop- 
ing. The scale is the same for all frames. From Naeini et al. 
(1999). 



by Timusk and Statt (1999), Loram and Tallon (2001) and Basov 
and Timusk (2005). A thorough study of the dc and Hall trans- 
port properties in various compounds performed on high-quality 
samples of the last generation was presented recently by Ando 
et al. (2004) 

(Blumberg et al., 1999; Hackl et al., 1996; Katsufuji et al., 1994; 
Nacini et al, 1999; Opel et al, 2000; Rcznik et al, 1993; Slakey 
et al., 1991; Staufer et al, 1990; Yamanaka et al., 1996) 
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ties is essentially determined by it/t'Y with t and t' the 
nearest and next-nearest neighbor hopping matrix ele- 
ments, respectively (Einzel and Hackl, 1996). Changing 
the ratio {t/t') enough to account for the observed Ra- 
man intensities leads to unrealistic band structures. The 
suggested importance of correlation effects on the Raman 
intensities was pointed out early on. 

The intensity changes of the Big spectra as a func- 
tion of doping (Fig. 42) are accompanied by a qualitative 
change of the temperature dependence at a given doping 
level. For Bi2Sr2CaCu208+5, the effect is shown in the 
range 0.15 < p < 0.23 (Fig. 43). For B2g symmetry, there 
is little change of both the overall intensity and the ini- 
tial slope which depends on temperature as expected for 
a metal. As opposed to symmetry, the temperature 
dependence of the Big spectra reverses sign (Fig. 43 (a)), 
indicating non-metallic behavior at p = 0.15. As outlined 
above (Eq. (55)), the initial slope of the response function 
is proportional to a transport lifetime r or, equivalently, a 
conductivity. Accordingly, the inverse of it, [dx'/Jdft]~^ 
is a resistivity which, for the Raman selection rules, is 
momentum sensitive as indicated by the index fi. 

It has been shown by Opel et al. (2000) that both the 
energy dependences and the magnitudes of T^{^,T) = 
h/Tf,{n,T) and m*(f],r)/m& = 1 + A^(f7,T) can be de- 
termined from the spectra using a memory function anal- 
ysis (Gotze and Wolfle, 1972) in combination with an 
energy integral over x'/i/^ ("sum rule"). By extrapo- 
lating the results to = 0, very reliable numbers for 
r{n 0,r) can be obtained.'"' 

In Fig. 44, the "Raman resistivities" r^(rj,T) in the 
limit r2 = are plotted for /z = Big, B2g at various tem- 
peratures and doping levels. The B2g results compare 
well to conventional resistivities p{T) using the Drude 
expression r^(f2 = 0,T) = e o {oj pi)^ p{T) with the renor- 
malized plasma frequency a)pi. With the plasma fre- 
quency in eV and the resistivity in /ificm, one finds 
r^(r) ~ 1.08(i?pi)V(r). At p = 0.23, the transport 
is essentially isotropic. Below p = 0.22, anisotropy de- 
velops quite abruptly, and for p < 0.16, the tempera- 
ture dependence of the Tb^ = 0, T) becomes non- 
metallic. This crossover behavior has been interpreted in 



The memory function or extended Drude analysis is particularly 
useful for a single-component response. Further details and lim- 
itations are discussed in the appendix of reference (Opel et al., 
2000). In Raman scattering the integral is as crucial as the 
plasma frequency in IR spectroscopy if magnitudes and energy 
dependences are to be derived in a similar way as in the analysis 
of the reflectivity (Basov and Timusk, 2005). "Sum rule" might 
be somewhat misleading and should be used with care since the 
integral over x'lii^)/^ <^ a' {Q,) (see Eq. (6)) is not a conserved 
quantity like the number of carriers in the famous f-sum rule. 
Without the integral the energy dependence of one of the quan- 
tities must be dropped and/or a fit to model functions for r^(S7) 
and 1 -I- A^(Q) is required (Blumberg et al., 1999; Hackl et al., 
1996; Naeini et al., 1999; Slakey et al., 1991; Yamanaka et al., 
1996). 



terms of an unconventional metal-insulator transition at 
0.20 < p < 0.22, where the antinodal transport is grad- 
ually quenched while the nodal one remains essentially 
unaffected (Venturini et al, 2002a). 

Since the current vertex has a similar k dependence as 
the B2g Raman vertex, IR spectroscopy and, similarly, 
dc transport project out mainly the nodal part of the 
Fermi surface (Devereaux, 2003b), and show therefore 
good qualitative agreement with the B2g Raman results, 
which exhibit little variation with doping for 0.1 < p < 
0.23 beyond the change of ujp\, as shown in Figure 44. 
The suppression of the antinodal transport goes along 
with a reduction of the quasiparticle strength in the 
vicinity of (tt, 0) (Kim et al., 2003), and a collapse of 
the Korringa law {TiT)~^ = const as observed by NMR, 
where Ti is the spin-lattice relaxation time (AUoul et at, 
1989; Billinge et al., 2003). 

Below p ~ 0.16, insulating behavior can also be ob- 
served in conventional transport at low temperature if 
superconductivity is suppressed with high magnetic fields 
(Ando et al, 2004a; Boebinger et al, 1996). The loga- 
rithmic divergence of the resistivity in the limit T — + 
indicates localization of the carriers. The crossover from 
metallic to non-metallic behavior occurs at temperatures 
T < Tc, well below those found in Big Raman scatter- 
ing. This is consistent with the observation that the i?2g 
spectra do not show any anomaly in the normal state. 
The synopsis of NMR, conventional, and Raman trans- 
port leads to the quite consistent conclusion that carriers 
get gradually localized upon decreasing doping. Accord- 
ing to the Raman results, the suppression of free-carrier 
transport starts at p ~ 0.21 at the anti-nodal regions of 
the Fermi surface and gradually moves towards the nodal 
points with decreasing p. The conductivity apparently 
disappears only very close to or at zero doping, making 
connection to recent ARPES results which show quasi- 
particles along the nodal direction even at the lowest (yet 
finite) doping level (Yoshida et al., 2003). 

Although there are no studies of both the doping and 
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FIG. 43 Direct comparison of the low-energy Big and B2g Ra- 
man response functions measured at temperatures indicated 
(all above Tc) in Bi2Sr2CaCu208+a for different doping levels 
p = 0.15, 0.20, and 0.22 (from left to right). From Venturini 
et al. (2002a). 
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temperature dependences in n-doped cuprates, the com- 
parison of the low-temperature data of Pr2_xCexCu204 
(Fig. 39) to those in BiaSraCaCuzOg-HS (Fig. 40) or 
La2-xSrxCu204 (Fig. 42) is worthwhile. In fact, a sig- 
nificant decrease of the overall Big intensity for small n 
is revealed whereas the continuum at B2g symmetry is 
not as doping independent as on the p-doped side, but 
rather becomes weaker towards lower doping. Raman 
relaxation rates have been determined explicitly only 
for Ndi.85Ceo.i5Cu204 (Koitzsch et ai, 2003). Both in 
Big and i?2g symmetry, Raman "resistivities" are close 
to those found by conventional transport, as shown in 
Figure 45. This is certainly at variance with the re- 
sults in hole-doped systems at comparable carrier con- 
centrations. From the viewpoint of Raman scattering, 
n-doped cuprates at optimal doping in the normal state 
look more like over-doped p-type samples with isotropic 
transport properties. The interpretation is not clear at 
the moment, but this similarity may be perhaps related 
to the rather involved band structure on the rt-doped side 
(Fournier et ai, 1997; Onose et ai, 2004). 



2. Quantum Critical Point(s) 

The quite unusual transport properties described 
above are accompanied by various crossover phenomena, 
such as the opening of a pseudogap (Timusk and Statt, 
1999), the collapse of the Korringa law (Billinge et ai, 
2003), and several other intriguing observations (Loram 
and Tallon, 2001). Recent studies show pseudogap phe- 
nomena also in electron doped cuprates (Alff et ai, 2003; 
Onose et ai, 2001, 2004). All those anomalies fit into 
the broader context of a quantum critical point (QCP) 
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FIG. 44 Raman relaxation rates ("resistivities") at B\g and 
Big symmetries of Bi2Sr2CaCu208+i . -Big and -829 is sensi- 
tive to anti-nodal and nodal regions of the Brillouin zone as 
indicated. The dashed lines represent transport data. Using 
plasma frequencies from IR spectroscopies, the Drude formula 
is used for the conversion from resistivities to relaxation rates. 
From Hackl et al. (2005). 



buried below the superconducting phase. A QCP oc- 
curs in the phase diagram at T = at a critical value Xc 
of a control parameter, such as doping x and/or pressure. 
The most general property of a QCP is the existence of 
thermal and quantum fluctuations up to high tempera- 
tures for a; ~ which prevent the transition into an 
ordered phase at finite T. The end point of the Neel 
phase at p ~ 0.02 as well as p = 0.05 and p = 0.27 delim- 
iting superconductivity are examples on the hole-doped 
side. However, an interesting putative QCP may be hid- 
den below Tc at the zero temperature extrapolated value 
of r*(p, n) (see Fig. 32), indicating competition between 
different types of order and superconductivity (Ander- 
gassen et ai, 2001; Castellani et ai, 1997; Chakravarty 
et ai, 2001; Kivelson et at, 2003; Sachdev, 1999). Below 
T* partial or even long range order can be established. 

In a normal metal, the kinetic energy of the electrons 
i?kin is much larger than the Coulomb energy U because 
of screening. If the ratio U /-Bkin increases upon decreas- 
ing carrier density and approaches 1 various types of in- 
stabilities can arise which usually induce a transition to 
an insulator. Classical examples are the Mott transi- 
tion or the Wigner crystal. The cuprates are close to 
this limit, and several additional possibilities have been 
discussed. Prominent examples are spontaneous orbital 
currents (Varma, 1997), spin ordering (Kivelson et al, 
2003; Machida, 1989; Tranquada et ai, 1995, 2004; Za- 
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tries of Ndi.85Ceo.i5Cu204 (Koitzsch et al, 2003). Big and 
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For general references see Sachdev (1999); Vojta (2003). 
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anen and Gunnarsson, 1989), charge ordering such as 
stripes or density waves (Andergassen et ai, 2001; Castel- 
lani et ai, 1995; Chakravarty et ai, 2001; Kivelson et ai, 
2003) and Fermi surface deformation fluctuations (Met- 
zncr et ai, 2003). In all cases, static order is found only, 
if at all, at very low doping and/or in specifically modi- 
fied structures such as La2_y_x^eySrxCu04 for y ~ 0.4 
and Re = Nd, Eu (Klauss et ai, 2000; Tranquada et ai, 
1995) or the nickelates. In Lai.775Sro.225Ni04 (Pashke- 
vich et at, 2000) and Lai.67Sro.33Ni04 (Blumberg et ai, 
1998; Yamamoto et ai, 1998) the formation of static dis- 
tortions due to the formation of stripes are big enough 
to split phonon and magnon lines due to zone doubling 
and to suppress the Big continuum at low energies. 

Of course, static order is much easier to verify than 
fiuctuating order. Nevertheless, fluctuating incommen- 
surate spin order has been observed in La2-xSrxCu04 
for X > 0.055 (Fujita et ai, 2002), YBaaCuaOe.ss, and 
YBa2Cu306.6 (Hinkov et ai, 2004), while commensu- 
rate spin order is seen at low energies in YBa2Cu3 06.353 
(Stock et ai, 2006). This means that in a temperature 
dependent volume characterized by a coherence length 
^s(r) a superstructure of the antiferromagnetically or- 
dered spins is established for higher doping levels. At 
the same time the charges which are expelled from the 
antiferromagnetic regions becoming spatially organized 
as well. This type of order has similarity with a liquid 
crystal and is sometimes referred to as nematic (Kivel- 
son et ai, 1998, 2003). From these experiments it can- 
not be decided whether the fluctuating magnetic super- 
structure is a property of the spins or of the charges. 
For statically ordered La2--y_xNdySrxCu04 with doping 
X = 1/8, Tranquada et al. (1995) showed that charge or- 
dering precedes spin ordering when the temperature is 
reduced. It has indeed been shown that charge ordering 
fluctuations are quite common phenomena in correlated 
systems (Castellani et ai, 1995; Metzner et ai, 2003). 

As long as there is no long range or static order which 
competes with superconductivity, the fluctuations can es- 
tablish Cooper pairing (Chakravarty et ai, 2001; Perali 
et at, 1996). For this reason, the understanding of the 
dynamics of incipient charge and spin order is a very in- 
teresting subject in the cuprates. 



3. Role of Fluctuations and Incipient Ordering Phenomena at 
small doping 

The study of dynamical order requires inelastic probes 
typically sensitive at finite q. In the case of charge or- 
dering, resonant X-ray scattering is the most promising 
as it can provide direct evidence for the charge order- 
ing at a particular wavevector (Abbamonte et ai, 2004), 
while neutrons can be used only if the charges modulate 
the lattice or the spin structure. While optical meth- 
ods are confined to q = two excitations with oppo- 
site momenta can be created, such as multi-phonon or 
two magnon scattering. Therefore light can be scattered 



from Cooper pairs or if an electron-hole pair couples to 
two excitations. It depends on the individual context 
how fluctuation phenomena are best studied. 

In the copper-oxygen systems, an unexpected addi- 
tional component in the Big spectra was observed in 
Lai.9oSro.ioCu04: as shown in Figure 46, a peak is found 
in the 100 cm~^ range which gains intensity and moves to 
very low but flnite energy with decreasing temperature 
(Tassini et al, 2005; Venturini et ai, 2002b) ."^ At high 
energy for all temperatures, and at high temperature for 
all energies, the spectra of the cuprates with p ~ 0.10 
in general are similar and show a strong anisotropy 
between nodal and anti-nodal regions in the Brillouin 
zone (see Fig. 42, 43, and 44). At low temperature 
the low energy feature is unique to Lai.9oSro.ioCu04, 
while the spectra become flatter in YBa2Cu306.5 and 
Bi2Sr2(Yo.38Cao.62)Cu208-i-5 (see Fig. 43). To appreci- 
ate the differences, one has to recall that the conventional 
transport properties of all these compounds are very sim- 
ilar (Ando et al., 2004). Before we discuss possible inter- 
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FIG. 46 Low energy response of underdoped La2-xSrxCu04. 
A Drude-like peak (Zawadowski and Cardona, 1990) with a 
characteristic energy i}c{x,T) is revealed after subtraction of 
the 2D response of the Cu02 planes. At a; = 0.02 (a) and 
0.10 (b) the additional response is observed in B2g and Big 
symmetry, respectively. The styles of the lines (colours) do 
not correspond to similar temperatures but rather highlight 
the scaling of the response with temperature: similar spectra 
are obtained if the temperatures differ by approximately a 
factor of 2. From Tassini et al. (2005). 



We note that the the lines observed recently by Gozar et al. 

(2004) in La2-xSrxCu04 {x = 0, 0.01) have nothing in com- 
mon with the response discussed here and were clearly identi- 
fied as one-magnon excitations (see also Silva Neto and Benfatto 

(2005) ). 
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pretations we wish to compare results for different doping 
levels and materials. 

From the comparison of the low energy peak to in- 
frared results, a relationship to fluctuating stripe or- 
der was conjectured (Vcnturini et ai, 2002b). For fur- 
ther support, it seems worthwhile to exploit the selection 
rules Raman scattering offers. The effect in Big symme- 
try is indeed compatible with the orientation of stripes 
along the Cu-0 bonds at p — x — 0.10 (Fujita et ai, 
2002). This is simply because order along the principal 
axes of an essentially tetragonal lattice corresponds to 
an orthorhombic distortion. Fluctuations correspond to 
micro- twinning, meaning that x and y cannot be accessed 
individually. Hence, Big symmetry measuring only the 
difference xx — yy (on a microscopic scale) projects out 
exactly this distortion. As a familiar example, we re- 
call that fully oxygenated YBa2Cu3 07 has Cu-0 chains 
along the crystallographic b axis making b > a. Even if 
the sample is twinned, the chain contributions are always 
superimposed on the Big spectra, while the xx and yy 
spectra are equal (as opposed to a single-domain crystal) . 
For X < 0.055, a reorientation by 45° is observed by neu- 
tron scattering (Fujita et ai, 2002). Since Big and B2g 
are equivalent modulo a 7r/4 rotation in the basal plane of 
a tetragonal lattice, the related peak should now appear 
in B2g rather than in Big symmetry. In fact, this has 
be observed in Lai.98Sro.o2Cu04, as shown in Figure 46 
(Tassini et ai, 2005). 




FIG. 47 Raman response x'/i.i'^yT) of 

(Yo.97Cao.o3)Ba2Cu306.i in Big (a) and B2g (b) sym- 
metry. The doping level is close to p = 0.02. From Hackl 
et al. (2005). 



Spectra similar to those in Lai.98Sro.o2Cu04 (Fig. 46) 
are found in (Yo.97Cao.o3)Ba2Cu306.i with p ~ 0.02 
(Fig. 47 (b)). Owing to the homogeneity of the sam- 
ple, the peak is very narrow and well defined. As the low 
energy peak emerges, the intensity of the continuum is 
suppressed by roughly 30 % over an energy range of ap- 
proximately 550 cm~^ in a similar though much stronger 
way than at higher doping (Nemetschek et ai, 1997; Opel 
et ai, 2000). In Big symmetry, there is an overall loss 
of spectral weight up to much higher energies, similar 
to YBa2Cu306+x and Bi2Sr2CaCu208+5 at higher dop- 
ings, demonstrating the transition to a correlated insula- 
tor (Freericks and Devereaux, 2001). 

Simultaneously with the studies in Lai.9oSro.ioCu04, 
the quasi ID ladder compound Sri4_a;Caa;Cu2404i was 
investigated (Blumberg et ai, 2002; Gozar et ai, 2003) 
which is considered a model system for the high- Tc 
cuprates (Dagotto and Rice, 1996; Dagotto, 1999; Sigrist 
et ai, 1994). At temperatures above approximately 
450 K low energy spectral peaks similar to those in 
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FIG. 48 Temperature dependent Raman response for cc po- 
larization in ladder compounds. Upper inset: fit of the data 
with Eq. 49 plus a small background. Lower inset: Two-leg 
ladder structure. From Gozar et al. (2003). 
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Lai.goSro.ioCu04 arc found (Fig. 48), although their en- 
ergies and widths are an order of magnitude smaller. The 
width obeys an Arrhenius law with a doping independent 
activation energy A of approximately 2100 K. In the in- 
sulator (a; = 0), the conductivity (for T < 300 K) reveals 
a similar A while for x = 12 the conductivity is metallic 
above 70 K (Gozar et al., 2003; Gozar and Blumberg, 
2005a) . The position of the maximum and the width de- 
crease upon cooling. In both cases, for T < 450 K the 
peaks move below the detection limit of 2 cm~^. Be- 
low 250 K, no indication of the peaks can be detected 
any more. It is interesting to note that at least in un- 
doped Sri4Cu2404i a charge-ordered state develops be- 
low Tco « 250 K (Abbamonte et al, 2004). 

In the ordered state Sri4Cu2404i has an optical re- 
sponse which is well described by that of a pinned charge 
density wave (Blumberg et al., 2002; Littlewood, 1987). 
For T > 250 K, the results are interpreted in terms of a 
damped plasma oscillation above Tco (Blumberg et al, 
2002; Gozar et al, 2003; Gozar and Blumberg, 2005a). 
Since free carriers damp the mode, the width is expected 
to increase proportional the conductivity just opposite to 
what one would expect from free carrier response. This 
complicates the interpretation of the results in metal- 
lic Sr2Cai2Cu2404i where the width still increases with 
temperature while the conductivity decreases. Gozar 
et al. (2003) argue that different regions of the Fermi 
surface contribute to the transport and the damping in 
essentially 2D Sr2Cai2Cu2404i. 

A comparison of panels (a) and (b) of Figure 48 shows 
that the continuum surviving at low temperature is sig- 
nificantly stronger in the metal (Fig. 48 (b)). We there- 
fore think that the low-energy mode found at T > 250 K 
could also originate from incipient charge order. In 
Sri4_a,Caa,Cu2404i the mode is then superimposed on a 
free electron response in a similar way as in the high- Tc 
cuprates. 

It is interesting that similar types of spectra are also 
found in Bii-ajCaj^MnOa (see section III.E), exhibiting 
comparable temperature dependences. While the results 
in ladders are interpreted in terms of an overdamped 
plasma mode in a CDW system above the ordering tem- 
perature (Blumberg et ai, 2002; Gozar et al, 2003; Gozar 
and Blumberg, 2005a), the response in the cuprates was 
proposed to originate from fluctuations of the charge 
density in the vicinity of a charge-ordering instability 
(Caprara et ai, 2005). Since the charge modulation ob- 
served e.g. by tunneling microscopy (Hoffman et al., 
2002; Howald et al., 2003; Vershinin et al., 2004) is at 
finite q, two fluctuations have to be exchanged to fulfill 
the q ~ selection rule in Raman scattering (Caprara et 
ai, 2005; Venturini et ai, 2000). Caprara et al. (2005) 
have shown that charge ordering fluctuations at incom- 
mensurate wavevcctors determined via neutron scatter- 
ing (Fujita et al., 2002) yield proper lineshapes and se- 
lection rules in Lai.9oSro.ioCu04 and Lai.98Sro.o2Cu04. 
To which extent the spin channel is involved and whether 
similar considerations apply to other dopings and cuprate 



families have not been explored yet and remain impor- 
tant future topics. 

In contrast to Sri4_a;CaxCu2404i, the low energy peak 
is strongly doping dependent in La2-xSrxCu04. The 
temperature scale is different by a factor of 2 for the two 
doping levels studied (Fig. 46), implying that r*(0.02) « 
2T*(0.10). The doping level determined from extrap- 
olating T* to zero, T*{pc) = 0, yields 0.15 < Pc < 0.20 
in close vicinity of the QCP inferred from other exper- 
iments. It is therefore possible that the low energy re- 
sponse is related to the thermal and quantum fluctua- 
tions above a hidden critical point. 

In summary, the comparison of the cuprates studied 
suggests that there is a superposition of two anoma- 
lies in Big symmetry: (i) below p ~ 0.22 antinodal 
quasiparticles become localized in all compounds as can 
be observed at sufficiently high temperature; (ii) with 
decreasing temperature a well-defined peak develops in 
Lai.9oSro.ioCu04 at low energy. Although nothing com- 
parable can be resolved in Y- and Bi-based compounds at 
the same doping level, the additional response may not 
necessarily be absent. It rather can mean that the low en- 
ergy features are broader and shifted to higher energies, 
hence becoming very similar in shape to the response 
from the 2D planes. Since the relationship between su- 
perconductivity and spin and/or charge fiuctuations is a 
key issue in the cuprates further work seems worthwhile 
here. 



V. CONCLUSIONS AND OPEN QUESTIONS 

The results discussed in the course of this review 
demonstrate that Raman spectroscopy is and has been 
an invaluable tool to investigate the dynamics of strongly 
correlated electrons. Improvements in the experimental 
technique have opened up the field to a variety of new 
systems. Light scattering has offered unique insights into 
dynamics in different regions of the Brillouin zone, show- 
ing the development of ordering and the competition be- 
tween various phases as the role of correlations increases. 
Materials such as superconductors with charge-density 
wave order, correlated insulators, ruthenates, mangan- 
ites, and finally the superconducting cuprates show the 
rich variety of phenomena which have been unveiled via 
Raman investigations. Along the way, light scattering 
has considerably deepened our knowledge of dynamics 
and correlation effects, and has provided several key in- 
gredients towards the development of a comprehensive 
theoretical description of these materials. 

Summarizing the findings on correlated systems over 
the past several years, the key ingredients stemming from 
Raman investigations include: (i) The development of 
anisotropics as correlations arc increased. Symmetry- 
selective measurements provide a tool to zoom in on the 
dynamics in different regions of the Brillouin zone. This 
is evidenced by Raman scattering studies on MgB2 as 
well as the cuprates in the normal and superconducting 
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phases, and the ordered phases in the ruthenates and 
the manganites. (ii) The existence of collective modes. 
Using Raman techniques the prevalence of certain types 
of order as a consequence were identified. Evidence for 
the importance of collective modes comes from CDW- 
superconductors and magnon scattering in antiferromag- 
netic insulators, and implications for a number of possible 
fluctuation or ordering modes exist in the cupratcs. (iii) 
A qualitative understanding of quantum critical behav- 
ior. The symmetry projection of parts of the Brillouin 
zone completes the picture of the battleground between 
competing orders concomitant with an underlying quan- 
tum phase transition. Here, spectral weight transfers as a 
function of temperature, doping, and pressure on a num- 
ber of materials combined with polarization dependent 
studies have opened a new door in the area of quantum 
criticality. 

It is clear that a number of issues are still of primary 
interest in correlated systems in general. These include 

(i) the origin of the electronic continuum in a number of 
compounds which look surprisingly similar at first pass. 

(ii) The origin of the polarization dependence in mate- 
rials which exhibit instabilities towards ordered phases. 

(iii) The mapping of Brillouin zone-projected electron dy- 
namics close to a quantum critical point. 

The cuprates continue to provide a wealth of infor- 
mation and puzzles in the area of superconductivity and 
strong electronic correlations. Issues in which consen- 
sus has been reached include (i) the presence antiferro- 
magnetic correlations over a wide range of doping levels 
evidenced from the two-magnon peak, (ii) the broad con- 
tinuum in the normal state as a common feature of all 
the cuprates, (iii) the d— wave nature of the pair state 
below Tc in a number of hole-doped materials derived 
from the low frequency power laws and polarization de- 
pendences, (iv) the dichotomy between the dynamics of 
Big and i?2g quasiparticles at low dopings, where the low 
frequency anti-nodal Big behavior is governed by inco- 
herence at the same time as the nodal B2g quasiparti- 
cles show relatively doping-independent metallic charac- 
ter, similar to the findings of transport quantities, and 
(v) the disappearance of this dichotomy for appreciably 
doped samples. 

In the cuprates there are still several unsettled ques- 
tions, (i) What is the origin of the Aig and/or Big peaks 
in the superconducting state? Do they originate from 
the redistribution of superconducting quasiparticles or 
are they collective modes or of lattice origin? (ii) How 
does the nodal/anti- nodal dichotomy picture of coher- 
ence/incoherence merge into metallic description at high 
doping, BCS superconductivity for temperatures below 
Tc, and antiferromagnetism near half filling? (iii) What 
is the microscopic origin for low energy peaks at low dop- 
ing, and how from this can information be obtained on 
the competition between ordered phases? 

These are issues which we believe will form the plan of 
development of Raman scattering in the cuprates as well 
as other materials in the years to come. The continua- 



tion of investigations on many new and well-characterized 
samples will further our knowledge of materials and cor- 
relations, and applications of the use of pressure and 
magnetic field will allow an exploration of quantum crit- 
icality and evolution of anisotropic electron dynamics in 
a variety of systems. The rapid development of Raman 
scattering as we have outlined indicate that the study of 
the electronic dynamics of complex materials will remain 
a vibrant and promising field of research. 
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